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ABSTRACT 


This  dissertation  presents  both  analytic  and 
processing  techniques  for  various  radar  imaging  systems. 

A  two  dimensional  system  classif ication  method,  which 
is  very  general  and  hence  applies  to  the  special  case  of 
radar  imaging  systems  as  well,  is  proposed  to  assist  in 
understanding  the  structure  and  describing  the  limitations 
of  2-D  systems.  Once  a  given  system  is  identified  with  the 
simplest  possible  class,  the  specific  techniques  can  be 
directly  utilized  to  process  the  data  or  reconstruct  the 
images . 

Following  a  review  of  radar  imaging  principles, 
several  coherent  radar  systems  are  analyzed  and  experimented 
upon.  They  include  synthetic  aperture  radar  (SAR)  ground 
mapping,  imaging  of  an  aircraft  target  from  turntable  data, 
and  imaging  of  a  flying  aircraft  target.  In  each  case 
the  point  spread  function  (PSF)  of  the  imaging  system  is 
derived  or  estimated.  Physical  considerations  are  then 
incorporated  in  mathematical  PSF's  to  categorize  the  imag¬ 
ing  systems  according  to  the  aforementioned  system  classi¬ 
fication  principle  proposed.  Degrees  of  Freedom  (DOF)  under 
different  imaging  geometries  are  analyzed  as  a  means  to 
determine  the  amount  of  information  present  in  the  usually 
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huge  amount  of  raw  radar  data  for  the  purpose  of  efficient 
computation  and  minimal  storage  requirement.  Motion 
compensation,  range  curvature,  range  alignment,  de-chirping, 
FFT,  registration  and  side  lobe  reduction  problems  are  all 
addressed  and  experiments  are  performed  using  data  from 
RAT-SCAT  (for  turntable  imaging)  and  other  facilities. 

The  results  shown  suggest  the  versatility  of  coherent  radar 
imaging . 

Possible  extentions  of  the  current  work  are  discussed. 
The  understanding  of  the  system  cha racteristics ,  in  parti¬ 
cular  the  formation  of  the  radar  image  will  aid  in  the 
advancement  of  techniques  for  radar  image  enhancement, 
encoding,  quantization,  and  restoration. 
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Chapter  1 
INTRODUCTION 

This  dissertation  deals  with  the  radar  imaging  problems. 
The  objective  of  a  radar  imaging  system  is  to  reconstruct 
the  target  image  from  its  radar  returns.  Radar  imaging 
geometry  is  one  of  the  most  important  factors  which  deter¬ 
mines  the  complexity  of  reconstruction  procedures  and  the 
quality  of  the  images  obtainable.  Many  radar  parameters 
function  interactively  making  the  analysis  and  the  optimiza¬ 
tion  of  the  system  very  complicated.  Physical  considera¬ 
tions  as  well  as  mathematical  assumptions  must  be  taken  and 
made  before  a  practical  reconstruction  algorithm  can  be 
realized. 

This  dissertation  presents  both  analytic  and  processing 
techniques  for  various  radar  imaging  systems  from  a  system 
noint  of  view.  Degrees  of  freedom  (DOF)  under  different 
imaging  geometries  are  sought  as  a  means  to  determine  the 
amount  of  information  available  in  the  usually  huge  amount 
of  raw  radar  data  for  the  purpose  of  efficient  image  recon¬ 
struction  and  minimal  storage  requirement. 

A  two-dimensional  system  classification  method,  which 
is  very  general  and  hence  applies  to  the  special  case  of 
radar  imaging  systems  as  well,  is  proposed  in  Chapter  2  as 


1 


a  preliminary  study  to  assist  in  understanding  the  structure 
and  describing  the  limitations  of  2-D  systems,  with  the 
ultimate  goal  of  providing  guidelines  for  any  given  system 
to  reduce  its  reconstruction  or  orocessing  effort.  Once  a 
system  is  identified  with  the  simplest  possible  class,  the 
specific  techniques  can  be  directly  applied  to  2-D  process¬ 
ing.  Mathematical  examples  are  given  and  physical  examples 
are  shown  to  fit  nicely  into  such  categori zations .  In 
addition  to  the  system  classification,  system  decomposition 
is  proposed  which  would  expose  more  system  structure  avail¬ 
able  for  current  fast  processing  techniques. 

An  overview  of  the  radar  signalling  principles  and 
processing  techniques  in  both  the  range  and  azimuth  direc¬ 
tions,  assuming  system  separabi 1 ity ,  is  presented  in 
Chapter  3.  This  chapter  is  devoted  to  reviewing  the  one¬ 
dimensional  radar  signal  processing  methods  which  form  the 
foundation  of  general  two-dimensional  radar  image  process¬ 
ing  discussed  in  later  chapters.  Advantages  as  well  as 
limitations  for  each  simple  case  are  given  and  compared. 

Following  the  chapter  on  radar  signalling  principles, 
several  coherent  radar  systems  are  analyzed  and  experiments 
performed  to  support  the  theoretical  work.  They  include 
synthetic  aperture  radar  ground  mapping  of  the  stripping 
mode,  coherent  imaging  of  aircraft  targets  from  turntable 
data,  and  imaging  of  an  actual  aircraft  target  in  flight. 

In  each  case  the  point  spread  function  of  the  imaging  system 
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is  derived  and  categorized  according  to  the  classification 
method  proposed  in  Chapter  2.  Physical  considerations  are 
then  incorporated  in  mathematical  PSF's  to  determine  the 
actual  reconstruction  algorithms  to  be  used.  Motion  compen¬ 
sation,  de-chirping,  range  alignment,  range  curvature,  image 
registration,  sidelobe  reduction  and  Fourier  transforming 
problems  are  all  discussed  and  numerous  experiments  per¬ 
formed  using  data  from  RAT-SCAT  (for  turntable  imaging)  and 
other  facilities.  The  results  shown  suggest  the  feasibility 
and  versatility  of  coherent  radar  imaging. 

Finally,  possible  extensions  of  the  current  work  are 
discussed.  Postprocessing  upon  the  images  reconstructed 
from  this  work  can  be  helpful  in  extracting  the  most  rele  - 
vant  information  once  the  characteristics  of  various  radar 
imaging  systems  are  fully  understood.  Among  them,  the  radar 
image  enhancement,  image  coding  and  quantization  are  valuable 
examples . 

An  analogy  between  radar  imaging  systems  and  sonar 
imaging  systems  from  the  mathematical  point  of  view  is 
pointed  out  which  could  be  a  significant  extension  of  this 
work.  However,  because  of  the  tremendous  disparity  between 
the  speeds  of  an  EM  wave  and  a  sound  wave,  many  different 
physical  considerations  have  to  be  investigated  before  any 
feasibility  of  the  sonar  imaging  system  could  be  concluded. 
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Chapter  2 


STRUCTURES  OF  TWO-DIMENSIONAL  TRANSFORMATIONS 

2. 1  Introduction 

Modern  day  imaging  systems  often  use  digital  computers 
for  completion  of  the  imaging  process  in  the  form  of  "recon¬ 
struction"  algorithms.  Such  algorithms  may  be  as  simple  as 
edge  enhancement  for  modulation  transfer  function  compen¬ 
sation  to  as  complex  as  polar  coordinate  tomographic  recon¬ 
struction  methods  for  three  dimensional  imaging.  However, 
one  common  thread  throughout  these  methodologies  is  the 
theme  of  avoiding  singularity  (ill-conditioning)  and 
consequent  noise  amplification  while  simultaneously  extrac¬ 
ting  as  much  image  information  as  is  possible  from  the 
original  scene  of  interest.  In  numerical  analysis  terms,  we 
wish  to  reconstruct  the  object  up  to  the  degrees  of  freedom 
(rank)  of  the  imaging  system  without  exceeding  the  computa¬ 
tional  complexity  indigeneous  to  a  particular  imaging 
configuration.  Previous  analysis  for  the  degrees  of  freedom 
of  imaging  systems  includes  the  earlier  work  of  Twomey 
(2-1, 2-2]  followed  by  Gori  and  Guattari  [2-3].  However,  the 
above  degrees  of  freedom  (DOF)  are  to  be  differentiated  from 
the  computational  degrees  of  freedom  (DOF^)  to  be  proposed 
in  this  chapter. 
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As  a  preliminary  study  for  our  research  work  of  radar 
imaging  systems,  this  chapter  is  devoted  to  a  generalization 
of  the  concept  of  the  point  spread  function  (PSF)  matrix 
and  its  role  in  general  two  variable  to  two  variable  dis¬ 
crete  linear  mappings  (the  most  common  of  which  are,  of 
course,  imaging  systems) .  The  importance  of  understanding 
the  system  PFS  lies  in  the  fact  that  the  image  reconstruction 
requires,  in  a  general  sense,  the  inversion  of  the  point 
spread  function  or  its  gramian  [2-1 ,2, 3,4 , 5] .  It  is  obvious 
that  one  would  like  to  make  full  use  of  the  structure  of  the 
PSF  to  ease  this  process.  By  categorizing  the  two- 
dimensional  systems  according  to  their  structure  with 
associated  computational  and  storage  savings  one  could  then 
apply  fast  and  efficient  reconstruction  algorithms  for  each 
class.  Thus,  whenever  a  given  category  is  identified  with 
the  simplest  possible  class,  the  underlying  reconstruction 
techniques  can  be  directly  utilized. 

The  PSF  matrix  is  generically  related  to  optical 
analysis  of  linear  space  invariant  imaging  systems,  but 
today  has  come  to  represent  a  much  more  general  mapping, 
especially  if  one  takes  a  numerical  analysis  view.  Because 
of  the  relatively  advanced  state  of  the  numeric  processing 
art  provided  by  numerical  procedures  in  problems  of 
inversion  of  ill-conditioned  systems  and  because  the  domain 
of  reconstruction  for  modern  day  imaging  systems  is  numeric 
and  discrete,  linear  algebraic  analysis  techniques  are 
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becoming  increasingly  popular  [2-5,6],  In  order  that  we 
take  advantage  of  such  techniques,  the  following  section 
discusses  the  implications  of  formulating  our  two  variable 
linear  mapping  systems  into  one  variable  lexicographic 
ordering  and  the  tensor  to  matrix  simplifications  thus 
introduced  (at  least  in  a  notational  sense) . 

Following  the  lexicographic  discussion,  nine  categories 
of  PSF  matrices  are  discussed,  along  with  their  computa¬ 
tional  degrees  of  freedom  (DOF  )  and  inherent  structural 
relationships.  It  is  shown  that  these  nine  models  cover  all 
two  variable  to  two  variable  linear  mappings  and  that  the 
reconstruction  and  analysis  algorithms  which  relate  to  the 
respective  models  allow  for  generalized  processing  results 
given  the  users'  knowledge  of  what  processing  model  his 
imaging  system  falls  into.  The  models  so  discussed  are  then 
exemplified  by  illustrations  from  radar  imaging,  tomography, 
optical  astigmatism,  curvature  of  field,  rotational  motion 
blur,  and  standard  space  invariant  optical  transfer  function 
distortions . 

2.2  Lexicographic  Ordering 

General  two  variable  to  two  variable  imaging  systems 
can  be  analyzed  from  three  typical  models;  a  continuous- 
continuous  representation,  a  continuous-discrete  represen¬ 
tation,  or  a  discrete-discrete  representation.  Such  models 
result  in  the  following  three  versions  of  the  integral 
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equations  used  for  imaging  (see  Fiaure  2-1). 

g  (x,y)  =  fft  (£,n)Mx,y;£,n)df;dn 
*-//  f  (£,n)h(£,n)d£dn 
g  =  [h] f 


(2-la) 

(2-lb) 

(2-lc) 


In  the  above  f  represents  the  original  object,  g  represents 
the  resulting  image  or  observations  to  be  reconstructed  and 
h  represents  the  linear  mapping  of  object  to  image.  The 
first  equation  represents  a  Fredholm  integral  with  a  four 
variable  kernel.  The  second  represents  a  linear  mapping 
with  a  continuous  -discrete  kernel,  and  the  third  repre¬ 
sents  the  mapping  in  linear  algebraic  notation.  For  the 
discrete  image  cases  we  will  index  the  (x,y)  continuous 
plane  with  a  set  of  discrete  indices  (i,k).  However,  for 
notaticnal  convenience,  these  indices  will  be  lexicogra¬ 
phically  ordered  to  form  a  "row  scanned"  vector  image  g. 

Thus  the  M*N  image  becomes  a  vector  of  dimension  MN*1. 

In  the  discrete-discrete  model  we  assume  the  original 
object  was  N*M  indexed  by  (p,q).  Due  to  our  lexicographic 
"row  scanning"  notation,  the  point  spread  function  becomes 
an  MNxMN  matrix  with  a  "block"  notational  structure. 

(Actually  imaging  a  systems  with  unequal  data  points  in  the 
input  and  output  planes  are  easily  included  in  this  analysis, 
but  for  notational  convenience,  we  will  restrict  our  nota¬ 
tion  to  non-square,  N*M,  but  equal  input  and  output  planes). 
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This  lexicographic  "block"  structure  of  the  point  spread 
function  matrix  can  be  represented  as 


[H] 


[H]  (1,1) 

[h](1'2) 

...  [h](1,m) 

[H]  {2,l) 

[h]<2'2> 

[H] (M'15 


[H] 


(M,  M) 


where  the  individual  matrices  are  N*N  given  as: 

[H]  (l,p)  =  [h.  ]  (1,P)  (2-3) 

K  t  CT 

See  Figure  2-1.  Thus  the  most  general  point  spread  function 
matrix  represents  an  arbitrary  kernel  h(i,k;p,q)  in  a 
discrete  notational  system.  We  will  see  that  the  block 
structure  of  equations  (2-2)  and  (2-3)  will  be  particularly 
powerful  as  an  analysis  tool  because,  in  matrix  algebra, 
such  block  structure  is  invariant  to  matrix  addition,  multi¬ 
plication,  transposition,  inversion,  etc.  Consequently, 
"block  diagonal",  "block  Toeplitz",  "block  circulant",  etc. 
matrices  will  retain  their  lexicographic  structure  through¬ 
out  manipulation  with  matrices  within  their  own  class. 


2 . 3  Point  Spread  Function  Matrices  and  Computational 
Requirements 

In  this  section  we  discuss  various  models  of  point 
spread  function  matrices,  define  the  notations  used  and 
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Fig.  2-1.  Imaging  models 


exploit  ways  to  utilize  the  structure  of  the  system  to 
speed  up  the  reconstruct  ions .  In  its  most  general  form 
[H]  is  that  of  a  nonseparable  space  variant  point  spread 
function  and  there  is  no  structure  without  further  analytic 
knowledge. 

For  our  discrete  model,  three  classes  of  structure 
could  be  used  to  simplify  the  system  description.  These 
structures  are  labeled  i)  space  invariant,  ii)  separate, 
and  iii)  separable.  Let  h(i,k;p,q)  be  the  point  spread 
function  with  input  variable  index  p,q  (as  in  Figure  2-1) 
and  output  variable  index  i,k.  Then  the  definitions  of 
the  three  concepts  mentioned  above  become: 

i )  Space  Invariant: 

invariant  in  one  dimension  only 
invariant  in  both  dimensions 

ii)  Separate : 

h(i,k;p,q)  =  hN (k;q/i)hM(i;p/q) 

Here  h^(i;p/q)  means  processing  along  the  i  or  p  direction 

with  q  as  a  parameter  fixed  for  i  and  p  but  variable  in  the 

k  indexed  dimension.  Note  that  in  h„  and  h„  above  the 

N  M 

conditioning  variables  i  and  q  are  output  and  input  varia¬ 
bles  respectively.  This  means  h^  occurs  before  h^.  Of 
course  it  is  also  possible  to  have  separate  processing  in 
the  other  order,  i.e. 

h(i,k;p,q)  =  hM ( i ; p/k ) hN (k ; q/p) . 
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The  order  of  processing  in  the  two  dimensions  is  not 
commutative,  in  contrast  to  the  separable  class  below. 

In  fact  the  separable  class  will  always  be  a  subset  of  the 
separate  class  . 
iii)  Separable : 

h  (i  ,k;p,q)  =  hM  (  i  ;  p)  h?J  (k  ; q) 

Combinations  of  the  above  three  structures  lead  to  nine 

cases  which  can  be  diagrammatically  depicted  as  in  Figure 

2-2  (a).  Thus  there  are  nine  branches  of  the  system 

classification  tree,  which  have  been  arbitrarily  numbered 

from  1  to  9  in  the  figure. 

If  we  define  the  computational  degrees  of  freedom 

(DOFc)  of  a  class  of  linear  system  matrices  to  be  the 

maximum  number  of  independently-determined  entries  of  the 

matrix  which  belongs  to  this  class,  then  intuitively,  this 

DOFc  of  the  matrix  is  an  indicator  of  the  effort  required 

to  diagonalize  the  matrix  or  to  reconstruct  the  image  in 

the  ideal  case.  Accordingly,  if  a  structure  is  a  subset  of 

another,  its  DOF  must  be  less  than  the  other.  In  the 
c 

following  we  exploit  the  structure  of  the  9  cases  of  Figure 
2-2  (a)  for  two  dimensional  discrete  linear  systems  in  the 
approximate  order  of  decreasing  DOFc  for  given  nontrivial 
M  and  N,  i.e.  Ni>2  and  N>2 .  We  will  assume  that  the  data  are 
stacked  in  appropriate  lexicographic  order  as  defined  in 
the  previous  section. 
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Fig. 2-2.  PSF  matrix  classification  and  decomposition 


Case  1;  Nonseparate  space  variant 


h(i,k;p,q)  =  h(i,k;p,q) 


DOFc((H])  = 


(1,1) 


...  [H] 


(1 ,  M) 


(H]  = 


(2-4) 


IH]'"'1’  ... 


where  [H]^'^  is  N  xN . 


Case  2:  Nonseparate  space  invariant  in  one  dimension 
(assume  invariant  in  (i,p)  coordinate) 

h ( i , k ; p , q )  =  h ( i-p , k ; q ) 

DOFc([Hj)  =  (2M-1 ) N2 


[H]  = 


(2-5a ) 


(M-l) 


Here  matrix  (2-5a)  is  Toeplitz  in  the  (i,p)  superscript 

2 

indices  but  each  submatrix  is  still  arbitrary  (thus  the  N 
factor  in  the  degrees  of  freedom) . 

Had  we  chosen  or  been  faced  with  a  linear  system  that 
was  nonseparate  but  invariant  in  the  (k,q)  dimension,  then 
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H (i , k; p,q)  =  h(i,k-q;p) 
DOFc[H]  =  (2N-1)M2 


( 2-5b) 


(i/P) 


and  aside  from  this  fact,  there  are  no  constraints  among 
different  [H] (l,p) '  s. 


Case  3:  Nonseparate  space  invariant  in  both  dimensions 

h { i , k ; p , q )  =  h(i-p,k-q) 

DOFc (  [H] )  =  (2M-1)  (2N-1) 

_[H](0)  [H]  (-X)  .  .  . 

[H]  (1)  [H]  (0) 

lH]  -  _  .  (2-6) 

*  • 

[H] (M~1)  ...  [H] (0) 
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where 


(HI 


(j)  _ 


1-N 


N-l 


0 


(j) 


The  above  matrix  for  [H]  is  known  as  block  Toeplitz  where 
the  superscript  indices  are  Toeplitz  and  the  submatrices 
are  also  Toeplitz.  This  matrix  structure  is  often  appro¬ 
ximated  by  block  circulants  which  are  then  diagonalizable 
by  the  2-D  Fourier  transform  process  [2-7]. 


Case  4:  Separate  space  variant 
(assume  column  operations  first) 


where 


h(i,k;p,q)  =  hfJ  (k  ;  q/i )  hM  ( i  ;  p/q) 
DOF ( [H] )  =  N2M+M2N 


[h]  =  ny  icM] 
(i) 


IV  = 


[R] 


[R] 


(2) 


[R] 


(M) 


tCM]  " 


[C](1,1)  ...  [c](1,M) 


[C](M,1)  _  [C](M,M) 


( 2—  7  a ) 
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where 


Here  [ ]  is  a  block  diagonal  and  [C^]  is  made  up  of  block 
submatrices  which  themselves  are  diagonal.  Had  we  reversed 
the  order  of  operations  for  case  4,  i.e. 

h  ( i  ,  k ;  p ,  q )  =  hM  ( i  ;  p/k)  hN  (k;  q/p) 
then  we  would  have 

[H]  =  [CM]  fRj,]  (2-7b) 

where  again  (R^l  is  block  diagonal  and  [CM]  consists  of 
diagonal  submatrices.  We  note  that  the  sparse  matrices  of 
the  form  [R^]  and  [CM ]  are  row  operators  and  column  opera¬ 
tors,  respectively,  on  lexicographically  ordered  two 
dimensional  data.  Decomposition  into  forms  like  (2-7a)  or 
(2-7b)  identifies  a  separate  system  with  sequential  row  and 
column  (or  vice  versa)  operations  which  are  not  necessarily 
commutable. 

Case  5:  Separate  space  invariant  in  one  dimension 
(assume  invariant  in  (i,p)  coordinate  and  column  operations 
first) 

h  (i,k;p,q)  =  hN (k ; q/i ) hM (i-p/q ) 

DOF  (  [H] )  =  N2M+  (2M-1)N 
c 
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[H]  =  [Rn]  [Cm] 


(2-8) 


where  [ ]  is  block  diagonal  as  before  but  [ C„ )  is  Toeplitz 


n  index,  i.e. 

[C]  (0)  [C]  (_1) 

( C]  ( 1  M) 

[C]  (1)  [C]  (0) 

• 

[C] (M_1)  .  .  . 

[C)(0> 

fSii  = 


and  the  [C]  ^  are  diagonal.  Again  if  the  role  of  the 
invariance  in  Case  5  is  reversed,  i.e.,  invariant  in  the 
(k,q)  coordinate 

h(i,k;p,q)  =  hN  (k-q/i ) hM (i ; p/q) 


then 


and 


DOFc[H]  =  (2N-DM+M  N 


[H]  =  CRnI  tCM] 

where  [Rj^]  is  block  diagonal  with  block  matrices  beinq 
Toeplitz  and  [C^]  is  subpartitioned  into  arbitrary  diagonal 
matrices  as  in  Case  4.  It  should  be  clear  that  Case  5  is  a 
subset  of  Case  4. 


Case  6;  Separate  space  invariant  in  both  dimensions 
(assume  column  operations  first) 


H(i,k;p,q)  =  hN(k-q/i)hM(i-p/q) 
DOF  (  [H] )  =  (2N-1)M+ (2M-1)N 
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[H]  =  [Rjj]  [CM] 


(2-9) 


Here  [ ]  is  block  diagonal  with  the  matrices  on  the 
diagonal  being  Toeplitz  and  [C^]  is  subpartitioned  into 
Toeplitz  indices  each  submatrix  of  which  is  diagonal. 
Therefore  Case  6  is  a  subset  of  both  Case  5  and  Case  4. 

Case  7:  Separable  space  variant 

h  ( i  ,  k ;  p ,  q )  =  h_M  ( i  ,  p)  (k ,  q ) 

Unlike  the  separate  cases  separable  operations  can  be 
performed  regardless  of  their  order.  This  is  because  of 
the  independence  of  one  dimensions 's  operations  upon  the 
other's.  Thus 

DOF  (  [H] )  =  M2+N2 
c 

[H]  =  feN]  0  [SM]  (2-10) 

where  0  is  the  kronecker  or  direct  product  matrix  operation 

[2-8].  Here  ]  is  N*N  and  [<3^]  is  M*M  resulting  in  the 
2  2 

M  +N  degrees  of  freedom.  Again  [r^]  and  [ *3^ ]  correspond  to 
row  and  column  operations  respectively.  An  equivalent 
imaging  system  in  non- lexicographic  order  for  this  and  all 
separable  cases  would  be  to  take  the  NxM  object  matrix  [f] 
and  form  the  NxM  image  matrix  [g]  as 

[g]  =  [C-Ml  If)  [RnJ 

or  in  lexicographically  stacked  notation 

g  =  [  [fiNJ  ®  [eMl  If 
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Case  8:  Separable  space  invariant  in  one  dimension 
(assume  invariant  in  (i,p)  coordinate) 

h(i,k;p,q)  =  hM (i-p) hN (k;q) 

DOF  (  [H] )  =  (2M-1)+N2 
c 

[H]  =  [«N]  ®  (CMJ  (2-11) 

where  [^  ]  is  now  a  Toeplitz  matrix.  If  the  space  invari¬ 
ance  were  in  the  (k,q)  coordinate,  then  we  would  have 

h (i ,k;p,q)  =  hM (i ;p) h^ (k-q) 

DOF  ( [H] )  =  M2+2N-1 
c 

and  [&N]  would  be  Toeplitz  rather  than  [C^]  . 

Case  9:  Separable  space  invariant  in  both  dimensions 

h  (i  ,  k;  p, q)  =  hM  (i-p)  hN  (k-q) 

DOF  (  [H] )  =  (2M-1)  +  (2N-1) 
c 

[H]  =  [ftN]  ®  IC-MJ  (2-12) 

where  both  and  are  Toeplitz.  This  particular 

case  is  the  well-known  model  for  separable  linear  space 

invariant  imaging  systems  in  which  approximations  to  [R^] 

and  [^-,.1  with  circulants  results  in  the  familiar  2-D  Fourier 
n 

transform  computation. 

2  2 

Note  that  DOFc([H])£M  N  for  all  nine  cases  above  for 
all  non-trivial  2-dimensional  systems  with  M^2  and  N>2. 

If  Cases  1-9  above  are  considered  as  sets,  then  there 
are  covering  relations  among  them  as  depicted  in  Table  2-1. 
From  the  table  it  is  clear  that  Case  1  covers  all  other 
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(read  from  left  to  top) 


=  :  equality 

Z>  :  strictly  covering 

*  :  partially  covering 

Examples:  Case  4  strictly  covers  Cases  5, 6, 7, 8, 9. 


Table  2-1.  Covering  relations  among  9  cases  proposed 
for  two-dimensional  systems 


cases.  Case  4  covers  all  higher  indexed  cases  (i.e.  5,6,7, 
8,9)  as  does  Case  7  (i.e.  8,9).  In  the  cases  of  space 
variance  in  both  dimensions  (Cases  1,4,7),  nonseparate 
case  (Case  1)  covers  separate  case  (Case  4)  which  in  turn 
covers  separable  case  (Case  7).  Similarly,  within  any 
group  of  the  nonseparate  or  separate  or  separable  cases  the 
space  variant  case  covers  1-D  space  invariant  case  which  in 
turn  covers  2-D  space  invariant  case.  They  are  straight 
covering  relations.  However,  nonseparate  space  invariant 
(1-D  or  2-D)  cases  do  not  cover  the  corresponding  separate 
cases,  e.g.  Case  2  ^6  Case  5.  This  is  because  in  the 
separate  cases  the  sequential  processing  property,  i.e., 
row  operations  followed  by  column  operations,  etc.,  allows 
a  more  general  definition  of  space  invariance  in  which  the 
invariance  is  only  in  a  column  by  column  or  row  by  row 
sense,  in  contrast  to  the  traditional  space  invariance  in 
which  the  processing  or  "blurring",  if  any,  is  uniform 
throughout  the  whole  image  plane. 

To  classify  a  given  system  into  the  simplest  possible 
case,  one  could  use  Fig.  2-2 (b)  for  sequential  tests.  The 
terminal  nodes  are  the  corresponding  simplest  cases  classi¬ 
fied.  As  an  alternative,  Fig.  2-2  (a)  can  easily  be  modi¬ 
fied  to  become  a  flow  chart  similar  to  Fig.  2-2 (b) .  The 
difference  is  that  in  Fig.  2-2  (a)  the  separability  property 
is  tested  first  while  in  Fig.  2-2  (b)  the  invariance  proper¬ 
ty  is  tested  first.  It  is  noted,  however,  that  to  identify 
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the  separate  (but  not  separable)  space  invariant  cases 
(Cases  4  and  5),  the  separability  has  to  be  determined 
first  due  to  the  reason  explained  above. 


The  knowledge  of  the  covering  relations  among  the  nine 
cases  is  important  in  the  application  of  system  decomposi¬ 
tion.  In  discrete  linear  image  reconstruction  or  object 
estimation  from  the  sensor  data  or  observations,  the 
relation  between  the  data  and  object,  i.e.  point  spread 
function  matrix  [ H ]  has  to  be  determined  from  the  data 
sensing  geometry  first  and  then  its  inverse  estimated  and 
multiplied  to  the  data  vector  g  to  come  up  with  an  estimate 
of  object  function  f_  as  in  Fig.  2-1.  This  reconstruction 
scheme,  although  straightforward  in  theory,  usually  involves 
huge  amounts  of  computation  because  of  the  complexity  of 
[H] .  In  the  worst  case,  one  would  have  to  resort  to  a  full 
singular  value  decomposition  (SVD)  to  find  a  pseudoinverse 
[HJ  [2-5].  Thus,  it  is  usually  desirable  to  decompose  [H] 
in  such  a  way  that  it  becomes  simpler  to  apply  consecutive 
processings  to  g  than  [H]  ^  itself.  One  possible  way  is  to 
do  some  preprocessing  [ H ^ ^  upon  g,  as  in  Fig.  2-2  (d)  such 
that  [ H ^  3  ^  and  the  remaining  processing,  [H^]  ^ ,  would  be 

much  simpler  than  the  [H]  ^  itself.  This  means  that,  if 

possible,  the  decomposed  component  matrices  should  be  in 
cases  covered  by  the  case  of  [H]  and  the  combined  effort  is 
still  less  than  ! H ]  itself.  Examples  of  pre-processing  in 
the  radar  imagirc  are  motion  compensation,  range  walking 


compensation,  data  presuming,  etc.,  as  in  [2-1]. 

Another  existent  technique  of  system  decomposition  is 
to  process  the  data  in  the  transformed  domain  as  in 
Fig.  2-2  (e) .  Sawchuk  [2-14]  proposed  a  geometrical 
transformation  of  [G]  in  the  reconstruction  of  rotationally 
blurred  images  with  successful  results.  Another  simple 
example  could  be  the  processing  of  polar  coordinate  Fourier 
transformation  which  originally  belongs  to  Case  2.  An 
interpolation  transformation  [G]  [2-17]  from  polar  coor¬ 

dinates  to  rectangular  coordinates  could  make  the  system 
separable. 

It  is  pointed  out  that  separability  and  invariance  are 
not  the  only  structure  that  can  be  utilized.  In  some 
situations,  the  structure  within  our  categorization  could  be 
further  scrutinized.  For  example  although  a  2-D  FFT  belongs 
to  Case  7  which  is  space  variant,  the  availability  of  fast 
algorithms  make  it  competitive  with  the  2-D  separable 
convolution  which  is  Case  9.  In  fact,  in  the  case  of 
space-invariant  processing,  sometimes  it  is  preferable  to 
apply  a  transforming  technique  similar  to  Fig.  2-2  (e)  to  do 
simple  point  by  point  operation  in  the  transformed  eigen- 
space,  i.e.,  the  Fourier  domain.  This  is  one  of  a  few 
examples  where  a  specific  structure  on  top  of  the  case 
structure  could  even  be  utilized  to  further  speed  up  the 
processing  within  that  case.  It  is  the  purpose  of  the 
system  classification  and  decomposition  to  extract  as  much 
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[H]  =  I G ]  ~ 1  [H3  ]  [G] 

IH ] ~ 1  =  [G] “ 1  [H3]_1  [G] 

f  =  [H]"1  g 

=  [G]"1  f  H  3  J _ 1  [G ]  g 

s.t.  [G]  and  [ H 3 ]  are  computationally  simpler 
than  [ H ] 


Figure  2-2 (e) .  Processing  in  transformed  domain 


Fig.  2-2.  PSF  matrix  classification  and 
decomoosi tion 
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system  structure  as  possible  for  the  purpose  of  fast 
processing.  However,  although  we  have  proposed  a  classifi¬ 
cation  method,  the  general  decomposition  is  the  state  of  art 
which  demands  ingenuity.  It  suffices  to  say  that,  in 
general,  the  principle  of  the  decomposition  is  such  that 
more  fast  processing  techniques  can  be  applied  upon  the 
transformed  data.  However,  some  disadvantages  might 
accompany  the  decomposition  of  the  system  in  that  one  might 
have  to  sacrifice  some  degrees  of  freedom  (in  the  informa¬ 
tion  sense)  for  fast  processing.  This  is  because  the 
optimal  mean-square  reconstruction  is  guaranteed  only  by 
the  availability  of  singular  values  of  [H]  itself  [2-6] , 
which  is  usually  very  difficult  to  infer  from  singular 
values  of  product  component  matrices. 

To  summarize  these  9  cases  for  linear  two  dimensional 
systems,  Table  2-2  has  been  prepared.  From  the  table  it  is 
clear  that  a  large  variety  of  systems  and  their  computa¬ 
tional  degrees  of  freedom  can  be  formulated  depending  on 
specific  point  spread  function  circumstances.  To  make 
these  9  cases  more  meaningful  and  relevant  to  actual  imaging 
conditions,  a  group  of  examples  are  developed  in  the 
following  section  as  illustration  of  some  of  the  cases 
discussed  here.  However,  before  proceeding  to  the  follow¬ 
ing  section,  the  reader's  attention  is  drawn  to  the  mathe¬ 
matical  examples  of  PSF's  of  Table  2-3  associated  with  the 
individual  case  numbers  of  Table  2-2.  While  these  are 
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Case  No. _ Description  I  PSF  *>DOF([})])  Equation 


Table  2-2.  P rope  rt  i es  of  the  S)  cases 


Point  Spread  Function 


h(i ,k;p,q)  =  exp{ - ( i +k+p+q) 2) 

h(i-p,k ;q)  =  exp{- ( ( i- p)  +k+q) 2} 
h(k-q,i;p)  =  exp{- ( ( k-q) +i+p) 2} 

h(i-p,k-q)  =  exp{-((i-p)  +  (k-q))2} 

hN(k;q/i)hM(i;p/q)  =  exp{- (k+q+i ) 2} exp{ - (i+p+q) 2} 

2  2 

hN(k;q/p)hM(i  ;p/q)  =  exp{-  (k+q+p) “} exp { -  (i+p+k)  } 

2  2 

hN(k;q/i)hM(i-p/q)  =  exp{-(k+q+i)  }exp{ - ( (i-p)+q)  } 
hN(k;q/p)hM(i-p/k)  =  exp{- (k+q+p) 2}exp{- (( i -q) +k) 2} 
hN(k-q/i)hM(i ;p/q)  =  exp{ - ( (k-q) +i ) 2}exp{ - ( i+p+q) 2} 
hN(k-q/p)hM(i;pA)  =  exp{-((k-q)+p)‘'}exp{-(i+p+k)2} 

hNA-q/i)hM(i-p/q)  =  exp{- ( (k-q)  +i ) 2)  exp{- ( ( i -p)  +q) 2} 
hN(k-q/p)hM(i-p/k)  =  exp{- ( (k-q) +p ) 2) exp{ - ( (i -p) +k) 2] 

2  2 

hN(k;q)hM(i;p)  =  exp{-(k+q)  }exp{- (i+p)"} 

2  ^ 

hN(k;q)hM(i-p)  =  exp{-(k+q)  }exp{- (i -p) "} 
hN(k-q)hM(i ;p)  =  exp{-(k+q)2}exp(-(i-p)2} 

2  2 

hN(k-q)hM(i-p)  =  exp{-(k-q)  }exp{- (i-p)  } 


Table  2-3.  Mathematical  examples  of  the  9  cases 


mathematically  correct,  they  may  not  be  meaningful  in  a 
physical  sense.  For  physically  meaningful  examples,  the 
reader  is  referred  to  the  following  section. 

2 . 4  Physical  Examples 

This  section  is  designed  to  put  the  numeric  structure 
of  the  point  spread  function  matrices  described  above 
into  physically  meaningful  perspective  with  respect  to  some 
typical  realistic  imaging  configurations.  While  naturally 
not  all  imaging  systems  can  be  investigated,  it  is  hoped 
that  enough  illustrations  are  presented  to  provide  a  firm 
grasp  of  the  relationships  of  the  real  world  physical 
phenomena  with  the  numerical  structure  of  the  computer  world 
computational  phenomena  in  implementing  restoration  and/or 
reconstruction  algorithms.  F,xamples  are  taken  from  optical 
imaging,  x-ray  imaging,  and  radar  imaging  and  range  in 
complexity  from  Case  1  (i.e.,  arbitrary  imaging)  to  Case  9 

(that  of  separable  space  invariant  imaging).  The  illust¬ 
rations  are  presented  below. 

Case  1:  Nonseparate  space  variant 

This  imaging  configuration  represents  an  arbitrarily 
comple-  system  in  which  the  object  field  f(p,g)  is  mapped 
into  the  image  field  g(i,k)  by  a  point  spread  function 
matrix  [H]  of  arbitrary  entries.  In  the  extremely  precise 
higher  order  modelling  of  any  linear  imaging  system,  the 
limiting  situation  would  result  in  this  case.  Naturally, 
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such  precision  is  not  warranted  in  many  situations,  thereby 
leading  to  the  simplification  of  succeeding  examples. 


Case  2:  Nonseparate  space  invariant  in  one  dimension 

Three  examples  of  imaging  systems  that  fit  into  this 
category  come  from  radar  imaging  and  x-ray  imaging.  In  the 
former  situation  consider  the  geometry  of  a  synthetic 
aperture  radar  (SAR)  operating  in  the  "stripping"  mode. 

This  configuration  is  illustrated  in  Figure  2-3.  After  some 
analysis  the  point  spread  function  can  be  shown  to  be  [2-9]: 


2,2.2,  fe 


/  2[(x1-x2)  +y1+z2] 

h(x2,t;x1,y1)  =  exp(joiot)A(x1-x2,y1,z2)fs(t - 1 


exp 


2[(x1-x2)2+y2+z2]% 
'3wc - : - 


(2-13) 


where  the  parameters  are  defined  as: 

ojq  :  range  offset  angular  frequency 
uic:  carrier  angular  frequency 

A:  antenna  beam  pattern 

f s :  modulating  function  of  signal  pulse 
z2:  flight  height 

c:  speed  of  light 

x2 :  takes  on  discrete  values  (nvTg) 

v:  aircraft  velocity 

Tg:  interpulse  period 

In  Chapter  4  Eq .  (2-13)  will  be  rederived  and  becomes 
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Eq.  (4-10).  With  (x^,y^)  being  the  input  or  object  plane 
and  (x2,t)  being  the  output  or  image  plane  we  see  that 

h (X2 , t ; x^ , y^ )  =  h (x2~x1 , t ; y1 )  (2-14a) 

or  in  our  discrete  numeric  notation  of  the  previous  section 

h ( i , k ; p , q )  =  h(i-p,k;q)  (2-14b) 

The  above  equation  is  referred  to  as  the  SAR  imaging 
equation  with  range  curvature.  The  range  curvature  forces 
the  point  spread  function  to  be  nonseparate.  However,  we 
will  see  in  subsequent  cases  that,  by  ignoring  such  curva¬ 
ture,  great  simplification  in  processing  can  be  achieved. 

The  second  example  for  this  case  is  obtained  from 
computer  aided  tomography  (CAT)  scanners  in  which  x-ray 
cross-section  slices  are  obtained  by  computer  reconstruction 
of  radiation  projections.  Figure  2-4  illustrates  the  geome¬ 
try  of  such  systems.  Essentially  the  imaging  equation  is 
given  by 

P  (r  ,  9  )  =  fj  f (x,y)h(r,9;x,y)dxdy  (2-15) 

where  the  integration  is  taken  over  the  unit  circle  as  in  [2-4]  . 
Assuming  the  blur  is  space  invariant  along  r  and  independent 
of  9,  McCaughey  and  Andrews  [2-4]  have  shown  that  once  r  and 
9  are  sampled  the  PSF  becomes 

h(r,9;x,y)  =  h  (xcos9+ysin0-r ) 
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f(x,y) 


r 


Fig.  2-4.  CAT  scanners 


/  \\ 


or 


h  ( r ,  0  ;  p ,  <// )  =  h  (pcos  (6-iJO  -r) 

where 

x  =  pcosiJj 
y  =  psinii) 


where  fixed  azimuthal  steps  are  assumed  and  the  super¬ 
scripts  denote  the  azimuth  variables.  E.g.  (i,m)  means 
(®i'^m)’  ^H]  is  c^rculant  the  index  (i-m)  because  of 

the  (6-ip)  dependence  in  the  PSF.  Note  that 

cos  ( 9  .  - iJj  )  =  cos(0  -ij>.).  Because  a  circulant  is  a  subset 
1  m  m  i 

of  a  Toeplitz  process  we  see  that  this  PSF  is  a  member  of 
Case  2,  equation  (2-5a). 

As  for  the  third  example,  consider  the  astigmatism  and 
curvature  of  field  studied  by  Sawchuk  [2-14] .  The  point 
spread  functions  behave  in  a  rather  complex  fashion. 

Figure  2-6  illustrates  the  imaging  geometry.  It  can  be 
shown  that  the  point  spread  function  assumes  the  form 
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h<xr'x0;ur'V 


D (2C+D) u 


4  , 


2  2 
[xrsin (x0-uQ ) ]  [xrcos  (xQ-u0 ) -ur ] 


2  2  4 
D  R  u 

r 


(2-17a) 


2  4  2 
(2C+D)u*R 


h (xr,xQ ;ur,u0 )  =  0,  otherwise 

where 

C:  degree  of  astigmatism 

D:  degree  of  field  curvature 

R:  radius  of  a  circular  exit  pupil. 

A  close  examination  of  Eq.  (2-17a)  shows  that  it  can  be 

rewritten  as 

h(x0 ,xr;u0 ,ur)  =  h(x0-u0,xr;ur)  (2-17b) 

or  in  our  discrete  numeric  notation 


h  ( i , k ; p , q )  =  h(i-p,k;q) 


In  the  above  formulation  it  should  be  noted  that  for  no 
field  curvature  (D=0)  and  therefore  only  astigmatism,  the 
point  spread  function  becomes 


ha(x0'xr;u9'ur1  = 


[4C2RV-(xr~ur)2]^ 

2C2u4 

r 


2  2 
u  -2CRu  <  x  <  u  +2CRu 
r  r  —  r  —  r  r 


and  zero  elsewhere,  which  is  separable  and  a  one  dimensional 
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blur  along  (xr»u  ). 


Case  3;  Monseparate  space  invariant 

This  situation  occurs  most  frequently  in  traditional 
linear  space  invariant  optical  imaging  systems.  The 
concepts  of  coherent  transfer  functions  (CTF) ,  optical 
transfer  functions  (OTF) ,  and  modulation  transfer  functions 
(MTF)  all  exist  in  this  category  of  imaging  situation.  As 
an  example  consider  the  OTF  for  a  diffraction  limited  system 
with  circular  exit  pupils  given  in  Figure  2-5.  Goodman 
[2-10]  has  shown  that  the  imaging  equation  are  given  by 


g (x,y) 


f (C,n)h(x-5,y-n)d^dn 


(2-18) 


for  a  circular  pupil  function  with  diameter  Z, 

T 


(y-o) ' 

17? 


h(x-£,y-n)  =  circ 
and  in  the  notation  of  the  previous  section 
h(i,k;p,q)  =  circ 


( 2-1 9a ) 


jVH 


~2  2 
-Pr  + (k-q) 


1/2 


(2-19b) 


Such  an  imaging  system  has  a  PSF  matrix  with  (2M-1)  (2N-1) 
degrees  of  freedom  arranged  in  a  block  Toeplitz  fashion. 

As  mentioned  earlier,  approximations  to  such  block  Toeplitz 
forms  with  block  circulant  forms  results  in  the  use  of  fast 
Fourier  transform  processing  [2-7]. 


Case  4:  Separate  space  variant 
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Fig.  2-6.  Astigmatism  and  curvature  of  field 


No  meaningful  examples  come  to  the  minds  of  the  author. 
Interested  readers  are  invited  to  contribute  to  this 
illustration  with  examples  from  their  own  experience. 

Case.  5:  Separate  space  invariant  in  one  dimension 

This  class  of  imaging  systems  has  examples  from  push- 
broom  sweep  arrays  (as  in  ERTS  and  LANDSAT  data)  and  SAR 
radar  imaging  configurations.  Proceeding  with  the  radar 
example,  referring  back  to  Figure  2-3  we  see  that  the 
effect  of  range  curvature  caused  the  SAR  model  to  be  non¬ 
separate.  However,  if  our  geometry  and  requirement  is  such 
that  the  range  curvature  effect  is  negligible,  then  our 
point  spread  function  becomes  decoupled  [2-9,2-11,2-12]  and 


h(x2,t;x1,y1)  =  hN  ( t ;  yL )  hM  (Xj-Xj/yj. ) 

=  exp(juot)A(x1-x2,y1,z2)fs 


(2-20a) 


exp  <-  ju 


{-jt°j((xrx2)2+yi+z2)2| 


(2-20b) 


which  will  be  repeated  as  Eq.  (4-12)  in  chapter  4.  In  this 

situation  with  N  samples  in  time  (range)  and  M  samples  in 

2 

azimuth,  we  obtain  (2M-1)N+N  degrees  of  freedom  (see 
Table  2-2).  While  this  system  is  "separate"  it  is  nonsepa- 
rable.  This  means  that  the  order  of  the  processing  is  still 
important.  Specifically,  one  should  process  the  range  data 
first  and  then  process  the  azimuth  data  in  the  image  recon¬ 
struction  effort.  The  azimuth  processing  is  a  function  of 
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range  which  explains  the  conditioning  on  in  equation 
(2-20)  and  the  need  for  the  conical  lens  in  optical  recon¬ 
struction  of  SAR  imagery.  (See  for  instance  Goodman  [2-10] , 
Harger  [2-11] ,  and  Rihaczek  [2-13] )  . 


Case  6:  Separate  space  invariant  in  both  dimension 

Two  examples  are  presented  for  this  case,  an  astigmatic 
processor  from  Goodman  [2-10]  and  rotational  motion  blur 
from  Sawchuk  [2-15] .  Figure  2-7  presents  the  optical  chain 
for  the  astigmatic  processor.  The  imaging  equation  is 

g(x3,y3)  -j(/’f(xry1)h(x3-x1)«(y3-y  )dx  dyx  (2-21) 


This  is  an  astigmatic  process  in  one  dimension  only.  Like¬ 
wise,  a  two-dimensional  (sequential)  astigmatic  process  is 
a  straightforward  extension. 

The  second  example  for  this  case  is  that  of  rotational 
motion  blur  as  illustrated  in  Figure  2-8.  Here  the  blur  is 
greater  the  further  away  from  the  origin.  The  blur  is  an 
angular  spreading  in  0  resulting  in  an  imaging  system  of 


g  (r ,  9  )  =  If  f  (p,<J>)h(p,<J>;r,9)pdpd<j>  (2-2  2a] 


where 


h(p,4>,r,0)  =  =-ha  ( 9  —  4> )  6  (p-r) 

P  w 


[  2-22b) 


and  the  factor  1/p  is  due  to  the  nonuniform  rotational 
speed.  In  terms  of  the  polar  coordinate  system  (p,<J>)  into 
(r,0),  the  imaging  system  is  a  blurring  along  9  only  and  is 
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two  dimensionally  invariant.  However,  in  terms  of  coor¬ 
dinate  system  (p,pi}>)  into  (r,r0)  which  has  a  Jacobian  of 
unity  with  respect  to  the  rectangular  coordinate  system, 
the  system  is  invariant  along  the  mapping  from  p  <j>  to  r0,  but 
changes  its  form  as  the  other  coordinate  p  or  r  changes. 
Intuitively  the  points  farther  away  from  the  origin  are 
subject  to  more  blur  as  a  result  of  nonuniform  rotational 
velocity  which  results  in  a  nonseparable  blurring  in  accor¬ 
dance  with  our  definition. 


Case  7:  Separable  space  variant 

Separable  space  variant  processes  occur  quite  frequent¬ 
ly  in  digital  computation  simply  as  a  convenience  for 
computational  simplicity.  Thus  two  dimensional  unitary 
transforms  for  image  coding  fall  into  this  example.  However 
more  relevant  imaging  situations  in  this  category  include 
imaging  lenses  with  square  apertures  and  the  Fourier  trans¬ 
forming  properties  of  lenses  themselves. 

Consider  a  lens  imaging  system  with  square  aperture 

such  that  p(x,y)  =  p  (x)p  (y)  is  separable  into  x  and  y 

x  y 

coordinate  and  magnification  factor  m  f  -1.  Goodman  12-10] 

has  shown  this  system  to  have  a  PSF  (see  Figure  2-9)  in 

coherent  light  of 

h(xi>y  ,xQ,y0)  =-lffp(x,y)exp 
\  dQdi 


j_:i§r(xi'+'Mxo)  x+  (yi+MV  [ dxdy 


(2  -23a) 
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(x^+Mx0) x  j  dx 
y*Py  (y)expj-j^-(yi+My0)yJdy  (2-23b) 

=  hx(xi;xO)hy (yi;y0}  (2-23c) 

Note  that  when  M=-l  we  obtain  Class  9  imaging. 

As  a  second  example,  consider  the  Fourier  transforming 
properties  of  lenses.  Figure  2-10  presents  the  optical 
system  for  this  illustration.  Here  the  output  image  (i.e., 
the  Fourier  transform  of  the  object)  is  given  by 

g(xf'yf)  =  3Xfy/f(xO'yo)expj"=>xT(xoxf+yoyf)|dxodyo 

(2-24a) 

Consequently 

h(xf,yf;x0,y0)  =  y^xp  J- j  ^xQxf  Jexp  j  -  j^yQyf  j  (2-24b) 

=  hx(x0'xf )hy  (y0'yf)  (2-24c) 

Case  8:  Separable  space  invariant  in  one  dimension 

For  this  case  we  return  to  our  SAR  stripping  mode 

radar  and  make  some  additional  simplifying  assumptions.  If 

the  difference  of  Y^max  and  y^m^n  on  ground  (Figure  2-3) 

is  fairly  small,  then  let  y.  be  the  nominal  value  of  y.  in 

In  1 

h{4  ^x2_xl//yl ^  °-  EcI-  (2-20a)  and  after  some  manipulations  we 
obtain  [2-9] 

h(x2,t;x1,y1)  =  hN(t;y1)hM(x2-x1)  (2-25a) 


X2d„d . 


/ 


px (x) exp 


jilL 
1  Ad . 
i 
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gUpYj)  image 


,y0)  object 

|* —  do 


Fig.  2-9.  Square  aperture  magnification 


Fig.  2-10.  Fourier  transforming  properties 


Case  9;  Separable  space  invariant  in  both  dimensions 

Pursuing  our  SAR  example,  a  further  simplification  of 

2  2  2 

the  range  offset  frequency  equal  to  zero  and  y^+z^  « 
results  in  the  following  PSF  [2-9] 


'  hn(t~4i) 


^*M  ^ X  2  -  x  P  "  Vt~t2,h>1tX2'Xl) 

(2-26) 


A  ^  1 

where  t0  =  -  . 

2  c 


A  second  example  of  a  separable  space  invariant 
imaging  system  is  a  television  raster  display.  The  horizon¬ 
tal  dimension  is  electronically  filtered  (due  to  the  frequ¬ 
ency  response  of  the  components)  and  the  PSF  of  the  electron 
gun.  The  vertical  dimension  is  blurred  or  "interpolated" 
by  the  PSF  of  the  electron  guns.  Both  dimensions  are 
space  invariant  in  a  well  tuned  monitor.  The  system  will  be 
separable  if  the  PSF  of  the  electron  gun  is  separable.  This 


is  usually  the  situation  for  a  Gaussian  CRT  spot. 

The  above  examples  are  tabulated  in  Table  2-4  as  a 
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summary.  It  is  clear  that  many  other  examples  exist,  but  in 
the  interest  of  brevity,  these  will  not  be  presented  here. 

2 . 5  Conclusions 

This  chapter  has  attempted  to  present  an  analysis  of 

the  numeric  computational  aspects  of  linear  imaging  systems 

from  the  viewpoint  of  the  computational  degrees  of  freedom 

of  such  imaging  systems.  A  matrix- lexicographic  structure 

was  developed  and  the  DOFc  of  the  PSF  matrix  was  then 

analyzed.  These  degrees  of  freedom  were  shown  to  range  from 

2  2 

the  most  complex  imaging  situations  (N  M  )  down  to  the 
simplest  imaging  systems  with  (2N-1) + (2M-1)  degrees  of 
freedom.  The  inherent  structure  of  each  PSF  matrix  was 
analyzed  for  9  cases  ranging  from  the  most  general  block 
structure  down  to  the  simplest  form  of  kronecker  products  of 
Toeplitz  matrices. 

Following  the  section  on  the  analysis  of  the  algebraic 
structure  of  each  PSF  matrix,  a  set  of  illustrations  from 
the  physical  world  of  imaging  systems  was  developed. 

Examples  from  x-ray,  radar,  coherent,  and  incoherent  optical 
systems  were  presented  as  illustrative  material  to  exemplify 
the  various  cases  developed  in  the  earlier  sections  of  the 
chapter.  It  is  hoped  that  this  chapter  has  aided  in  the 
ever  growing  marriage  of  computational  algorithms  and  their 
use  for  complex  image  reconstruction.  As  future  imaging 
systems  are  developed,  it  should  be  possible  to  utilize  the 
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Case  No. 


Reference 


Imaging  System 


Arbitrarily  Complex  Imaging 
System 


SAR  Stripping  Mode  with  Range 
Curvature 


[2-9] ,  [2-11] , 
[2-12] 


CAT  Scanner 


[2-4] 


Astigmatism  and  Curvature  of 
Field 


[2-14] 


OTF , MTF , CTF 


[2-5]  ,  [2-10] 


None  Presented 


Pushbroom  Sweep  Arrays 
SAR  no  Range  Curvature 

Astigmatic  Processor 
Rotational  Motion  Blur 

Lens  with  Square  Aperture 
Fourier  Properties  of  Lenses 

SAR  Stripping  Constant  Range 

SAR  Stripping  Low  Altitude 
TV  Displays 


[2-9] 

[2-10] 

[2-15] 

[2-10] 

[2-10] , [2-16] 
[2-9] 

[2-9] 


Table  2-4.  Physical  examples  of  imaging  systems 
for  the  9  cases 


framework  developed  here  to  readily  anticipate  the  comple 
xity  of  computational  reconstruction  of  imagery  once  the 
structure  of  the  PSF  matrix  is  discovered. 


i 

I 


Chapter  3 

PRINCIPLES  OF  RADAR  IMAGING 


3 . 1  Imaging  Systems 

The  general  two-dimensional  imaging  system,  as  depicted 
in  Fig.  3-1  is  a  transformation  of  the  target  or  object 
function  f  ( C  ,  n )  into  the  observation  function  g(x,y)  through 
a  point  spread  function  h(x,y;£,n),  which  could  be  linear 
or  nonlinear  and  its  variables  discrete  or  continuous  or  a 
mixture  of  both  [3-1] .  The  purpose  behind  the  imaging 
system  is  to  reconstruct  or  estimate  the  object  function  f 
as  close  as  possible  in  some  sense  from  the  observation  g. 
There  are  various  imaging  systems  in  the  real  world,  e.g. 
optical,  infrared,  x-ray,  and  radar  systems;  with  a  wide 
range  of  characteristics  and  limitations,  the  most  conven¬ 
tional  one  being  the  passive  optical  system  in  which  a 
camera  records  the  reflectivities  or  transmissivities  of  the 
object.  In  this  case  relatively  little  reconstruction 
effort  is  required  because  the  point  spread  function  is 
almost  space  invariant  and  in  most  situations  has  an 
impulse-like  shape  which  gives  a  g  closely  resembling  f  as 
one  usually  experiences  from  taking  photographic  pictures 
in  his  daily  life.  On  the  other  hand,  the  radar  imaging 
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Fig.  3-1.  Imaging  system  and  reconstruction 
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Fig.  3-2.  Antenna  illumination  geometry 
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system  uses  an  active  device  to  radiate  signals  whose 
returns  usually  have  to  be  processed  by  complicated 
procedures  to  reconstruct  the  object  function. 

3 . 2  Principles  of  Radar  Imaging  Systems 

There  are  three  major  physical  parts  in  a  radar 
system:  the  transmitter,  the  receiver  and  the  target.  The 

transmitter,  which  is  required  for  all  active  imaging 
systems,  is  used  to  illuminate  the  target.  The  receiver 
is  a  sensing  device  which  records  the  signals  reflected 
from  the  target  and  is  followed  or  accompanied  by  a  data 
processing  unit.  The  target  is  the  object  some  physical 
properties  of  which  are  to  be  calculated  from  the  received 
data.  In  fact  the  functions  of  the  radar  systems  can  be 
dichotomized  as  target  detection  and  parameter  estimation. 
Detection  of  a  target  is  the  determination  of  its  presence 
in  the  unavoidably  noisy  situation,  and  parameter  estimation 
is  the  measuring  of  the  characteristics  of  the  targets,  e.g. 
their  ranges,  velocities,  angular  positions,  sizes,  etc.  by 
the  extraction  of  available  information  from  the  received 
echoes  when  the  presence  of  the  targets  has  been  determined 
or  assumed  [3-2].  By  this  token,  the  radar  imaging  in  which 
shape  and  size  of  the  target  are  of  concern  belongs  to 
the  second  category. 

In  order  to  reconstruct  a  radar  image  of  some  target 
from  its  signal  returns,  two  prerequisites  have  to  be 
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satisfied  from  the  system  point  of  view.  First,  the 
returned  data  has  to  have  a  two-dimensional  format. 

Second,  radar  imaging  geometry  must  be  such  that  the  return 
from  each  pulse  or  signature  contains  different  information 
about  the  target.  In  fact,  degrees  of  freedom  (DOF) 
analysis  on  the  radar  returns  as  explored  in  some  detail 
in  this  dissertation  provides  an  attempt  at  evaluating  the 
information-carrying  capability  of  the  system  by  analyzing 
the  extent  to  which  the  above  conditions  are  satisfied. 

The  two  geometrical  coordinates  associated  with  the 
radar  or  the  targets  in  a  2-D  radar  imaging  system  are 
usually  called  range  and  azimuth  (or  cross-range) .  (In  3-D, 
there  is  one  more  called  elevation) .  Range  is  the  direction 
along  which  the  radar  signal  is  transmitted,  reflected  and 
received.  Range  information  is  provided  by  the  return  from 
a  single  pulse  resolved  by  timing  (when  a  short  pulse 
waveform  is  used)  or  range  compression  technique  (when  a 
linear-FM-like  waveform  is  used) .  Azimuth  is  the  direction 
orthogonal  to  the  range  direction  in  the  surface  of  interest. 
It  is  the  component  of  the  relative  motion  between  the  radar 
and  target  through  which  an  effective  target  rotation  is 
created. 

In  the  simplest  geometry  of  radar  imaging,  the 
processing  of  the  data  is  separable  along  the  azimuth  and 
the  range  directions.  Although  the  coupling  of  the  two 
dimensions  has  to  be  taken  into  account  in  a  complicated 
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radar  imaging  geometry,  this  chapter  will  discuss  processing 
principles  when  the  system  is  assumed  separable  and  hence 
decomposable  into  two  consecutive  1-D  processes  because  the 
understanding  of  processing  principles  along  individual 
dimensions  is  basic  to  comprehending  sophisticated  radar 
imaging  systems.  In  the  separable  radar  imaging  systems  the 
operations  in  the  two  dimensions  could  be  any  combination  of 
incoherent  and  coherent  processing  depending  on  the  wave¬ 
form  and  antenna  pattern  used.  Incoherent  processing 
assumes  incoherence  of  those  signal  returns  from  a  single 
target  point  and  makes  no  attempts  at  processing  the  phase 
part  of  the  returns,  whereby  coherent  processing  uses  the 
phase  relations  among  different  returns  to  "compress"  the 
pulse  to  achieve  high  resolutions,  as  will  be  explained. 

We  now  describe  the  principles  of  the  four  possible  opera¬ 
tions  along  the  range  and  the  azimuth  directions  indivi¬ 
dually. 

(A)  Incoherent  range  operation. 

The  returns  of  a  single  pulse  "rom  a  particular  target 
point  are  assumed  incoherent  in  phase.  Thus  no  compression 
technique  (as  this  would  require  the  knowledge  of  the  phase 
relation)  could  be  applied  and  the  only  range  resolving 
ability  comes  from  the  actual  timing  of  the  strong  magnitude 
pulse  return.  If  t  is  the  time  duration  of  the  pulse  then 
the  ideal  range  resolution  PT  will  be  half  the  time  on 
target  (TOT)  of  that  pulse  translated  into  the  spatial 
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domain,  i.e. 


1 

2 


CT 


(3-1) 


where  c  is  the  propagation  speed  of  the  EM  wave.  For 
example,  to  achieve  a  range  resolution  of  1  ft.  x  must  be 
less  than  2  nanoseconds.  Thus  for  high  range  resolving 
ability  very  narrow  pulses  have  to  be  used,  which  requires 
very  high  peak  powers  to  keep  a  sufficiently  large  signal 
noise  ratio.  Obviously  the  requirements  of  an  extremely 
short  pulse  and  high  peak  power  could  not  easily  be  met 
from  a  technical  viewpoint  as  the  resolution  pursued  becomes 
higher  and  higher.  The  main  difficulty  lies  in  the  fact 
that  the  duration  of  a  rectangular  pulse  is  inversely 
proportional  to  its  bandwidth  making  a  simultaneous  large 
bandwidth  and  a  high  SNR  very  difficult.  One  way  to  get 
around  with  this  conflict  is  to  use  phase  modulated  wave¬ 
forms  which  results  in  what  we  call  a  coherent  range  opera¬ 
tion  as  follows. 

(B)  Coherent  range  operation. 

In  fact,  it  is  the  bandwidth  B^,  and  not  the  time 
duration  of  a  pulse  waveform  that  determines  the  range 
resolution,  i.e. 


Equation  (3-2)  is  a  universal  relation  which  applies  to  the 
special  case  of  rectangular  waveforms  of  incoherent  range 
operation.  Combining  Eqs.  (3-1)  and  (3-2)  gives 
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1 


(3-3) 


for  a  rectangular  pulse  with  time  duration  x  and  bandwidth 
Br  and  without  any  phase  modulation.  The  quantity  Br  •  t  is 
called  the  time-bandwidth  product  (TBP)  which  is  an  indica¬ 
tion  of  the  compression  capacity  of  the  pulse  12-31. 

Coherent  range  processing  uses  phase  modulation  on  the 
transmitted  pulses  to  enlarge  its  TBP.  For  example  in  the 
case  of  a  linear-FM  modulation  it  can  be  shown  that  for 
large  bandwidth  the  signal  bandwidth  is  proportion  to  k  •  t 
where  k  is  the  linear  FM  rate  which  determines  how  fast  the 
frequency  is  changing  within  a  pulse  [3-3]  and  x  is  the 
pulse  duration.  Thus  in  contrast  to  Eq.  (3-3)  the  bandwidth 
and  the  pulsewidth  are  now  in  proportion  to  each  other  and 
the  S/N  will  not  suffer,  but  benefit  instead,  from  using  a 
larger  bandwidth  pulse  which  is  accompanied  by  a  longer 
pulse  duration.  Other  forms  of  phase  modulation  are  also 
possible  to  "stretch"  the  pulse  duration  while  increasing 
its  bandwidth,  but  usually  the  linear  FM  or  chiro  modulation 
is  preferred  because  of  its  high  efficiency  in  terms  of  TBP, 
easy  analysis  and  implementation.  Upon  receiving  the  echo 
of  a  modulated  pulse,  the  data  processor  "compresses"  or 
"decodes"  it  to  get  high  resolution  profiles.  The  resultant 
de-modulated  signal  is  in  the  decoded  time  domain  and 
consists  of  pulses  whose  width  or  resolution  is  determined 
by  the  bandwidth  of  the  transmitted  signal  and  whose 
locations  determined  by  relative  ranges  of  the  target 
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points.  The  compression  ratio,  or  the  ratio  between  the 
original  pulse  duration  and  the  duration  of  the  compressed 
pulse  from  a  point  target,  is  equal  to  the  TBP  in  the  ideal 
case  [3-3].  An  intuitive  explanation  of  the  high  resolution 
advantage  of  a  coherent  radar  processing  is  as  follows: 
Conceive  that  the  point  spread  function  in  the  range  direc¬ 
tion  of  e.  radar  system  is  essentially  a  scaled  form  of  the 
transmitter  signal  waveform.  Although  the  linear  FM  is  rela¬ 
tively  long  in  time  duration  and  hence  the  PSF  wide  in  the 
range  domain,  high  range  resolution  is  possible  because  of 
the  high  bandwidth  of  the  signal  due  to  the  phase  modulation. 

This  is  the  key  point  of  coherent  radar  signalling. 
Hence  the  simultaneous  achievements  of  a  high  resolution 
and  a  large  signal  energy  are  possible  by  a  phase  modulated 
long  duration  pulse.  It  is  noted  that,  however,  because  the 
transmitted  signal  is  in  a  coded  form  (e.g.,  linear  FM 
modulated),  a  decoding  scheme  (e.g.,  pulse  compression) 
which  is  a  coherent  processing,  has  to  be  applied  upon  the 
received  echoes,  in  contrast  to  noncoherent  processing 
(simple  range  grating)  for  straight  short  pulse  signalling. 

The  compression  or  decoding  of  the  return  of  a  linear 
FM  signal  from  the  target  can  be  decomposed  into  two  concep¬ 
tual  steps  [3-4]:  a  quadratic  phase  removal  followed  by  a 
Fourier  transforming.  The  former  is  to  compensate  for  the 
linear-FM  modulation  of  the  signal  and  the  latter  is  to 
localize  or  resolve  the  various  target  points. 
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Because  of  the  strong  mathematical  similarity  between 
a  linear  FM  signal  and  the  quadratic  phase  delaying  property 
of  the  lens  [3-8],  optical  or  hologram  processing  techniques 
can  be  used  to  compress  the  radar  returns  and  to  reconstruct 
the  target  images.  The  analogy  between  a  linear  FM  wave¬ 
form  and  the  phase  delay  of  a  lens  is  that  both  of  them  are 
quadratic  in  phase,  with  the  linear  FM  rate  analogous  to  the 
inverse  of  the  focal  length  of  the  lens.  Thus  a  lens  of 
suitable  focal  length  (or  a  conical  lens  if  the  focal 
lengths  change  in  a  linear  fashion,  as  described  in  a  later 
chapter)  can  be  used  to  remove  the  quadratic  phases  associa¬ 
ted  with  the  linear  FM,  followed  by  a  Fourier  transforming 
to  complete  compressing  the  signal  [3-3] . 

(C)  Incoherent  azimuth  operation. 

Before  the  development  of  the  synthetic  aperture  radar 
(SAR)  techniques,  incoherent  processing  was  the  only  way  to 
obtain  the  azimuthal  information  of  the  target  by  a  radar 
technique.  At  that  time,  the  traditional  imaging  radar 
achieved  azimuth  resolution  by  using  an  antenna  with  illu¬ 
mination  pattern  very  narrow  in  azimuth.  Physically  speak¬ 
ing,  the  narrower  the  beam  width,  the  "narrower"  the  point 
spread  function  and  thus  the  better  the  resolution  obtain¬ 
able.  There  is  no  coherent  processing  necessary  along  the 
azimuth  direction  because  the  illumination  patches  are 
independent  azimuthally  and  no  fixed  phase  relations  are 
maintained  among  different  pulses. 
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From  the  antenna  theory  it  is  well  known  that  the  half 
power  beam  width  8  in  radians  of  a  physical  antenna  of 
length  L  is  (See  Fig.  3-2). 

8  -  X/L  (3-4) 

Thus  fine  azimuth  resolution  (8  small)  demands  a  very  long 
physical  antenna  (L  large) ,  not  practically  available  in 
many  situations  as  higher  and  higher  azimuth  resolutions 
are  required  in  recent  years. 

(D)  Coherent  azimuth  operation. 

From  the  analysis  of  the  coherent  range  operation  in 
(B)  it  can  easily  be  conjectured  that  a  high  azimuth  resolu¬ 
tion  could  be  achieved  by  phase  coding,  e.g.,  a  linear  FM 
modulation  along  the  azimuth  direction  just  as  in  the  range 
direction.  However  since  the  object  of  interest  is  two 
dimensional  (range  and  azimuth)  while  the  signal  is  onl'  one 
dimensional  (a  function  of  time  only) ,  and  because  only  one 
dimensional  (i.e.,  the  range)  information  can  be  inferred 
from  a  pulse  return,  no  explicit  phase  modulation  upon  the 
transmitted  signals  along  the  azimuth  can  be  realized. 
Fortunately,  an  implicit  modulation  is  made  possible  by 
the  relative  motion  of  the  target  and  the  radar  which 
creates  quadratic  phase  history  to  each  target  point  in  the 
azimuth  direction.  In  fact,  the  range  history  of  any 
rotating  target  point  as  a  function  of  the  time  generates 
quadratic  or  linear  FM  phases  for  small  angle  changes  as  a 
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first  order  approximation  [3-5]  .  This  is  the  baseline  of 
the  principle  of  synthetic  aperture  radars  (SAR) . 

The  azimuth  resolving  ability  of  the  SAR  can  also  be 
explained  from  the  antenna  point  of  view.  The  radar 
radiates  pulses  at  different  target  aspect  angles  and 
receives  their  echoes  from  the  targets  shortly  after. 

Because  of  the  relatively  short  TOT  of  any  pulse  it  can  be 
assumed  that  during  any  single  TOT  the  motion  of  the  targets 
or  the  radar  is  negligible.  The  returns  of  different  pulses 
from  the  same  target  point  are  azimuthally  modulated  by  the 
relative  motion  between  the  radar  and  that  particular  point. 
The  induced  linear-FM  returns  can  be  compressed  to  yield 
high  azimuthal  resolution  as  in  the  range  case.  This  is 
analogous  to  the  case  of  a  physical  antenna  array  where  the 
received  signals  at  each  array  element  are  coherently 
processed  and  summed  to  yield  a  sharp  effective  antenna 
pattern  [3-6] .  The  analogy  between  a  physical  antenna  array 
and  a  synthetic  antenna  is  depicted  in  Fig.  3-3. 

From  Fig.  3-3 (b)  it  is  observed  that  the  whole  physical 
azimuthal  beam  width  at  range  R  gives  the  length  of  the 
synthetic  aperture  available  at  that  range.  If  L  and 
are  the  physical  and  synthetic  antenna  width  along  azimuth 
then 


L 


eff 


R*  B 


RA 

L 


(3-5) 


where  B  is  the  physical  antenna  width  in  radians.  Since  the 
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operation  of  SAR  utilizes  the  two-way  beam  pattern  in  the 
sense  that  the  phase  shift  is  introduced  on  both  the  paths 
to  and  from  the  targets,  the  round-trip  phase  shift  effec¬ 
tively  reduces  the  wavelength  by  a  factor  of  2.  Thus  the 
effective  antenna  pattern  width  is,  from  Eq.  (3-1)  , 


eff 


2L 


(3-6) 
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The  azimuth  resolution  p  is  the  effective  beam  width 

cl 


projected  on  the  target  at  range  R 


Pa  =  BeffR  = 


L 

2 


(3-7) 


which  is  proportional  to  the  azimuth  size  of  the  physical 
antenna,  and  is  independent  of  X  and  R.  Thus  in  the  SAR  to 
achieve  higher  azimuth  resolution  a  shorter  antenna  has  to 
be  used,  in  contrast  to  a  single  antenna  in  a  traditional 
incoherent  case  as  in  (C) .  This  is  because  a  smaller 
antenna  has  a  wider  illumination  pattern  which  generates 
a  wider  bandwidth  available,  a  characteristic  of  linear-FM 
signals  [  2-7 ] . 


3 . 3  Conclusions 

This  chapter  analyzed  one-dimensional  radar  signal 
processing  principles  assuming  system  separability.  In  a 
radar  system  coherent  in  both  the  range  and  azimuth  dimen¬ 
sions,  linear  FM  modulation  exists  in  both  dimensions.  In 
the  range  direction  it  is  created  by  actually  modulating 
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the  transmitted  waveform  whereby  in  the  azimuth  direction 
by  an  effective  rotation  of  the  target.  However,  the 
returns  are  processed  in  the  receiver  end  by  the  same 
compression  technique  to  resolve  the  target  points  in  both 
dimensions . 

As  the  analysis  goes  from  one  dimensional  processing  to 
a  two-dimensional  one  as  will  be  done  in  the  following 
chapters,  new  problems  could  generate  many  difficulties. 
Coupling  between  the  two  dimensions,  e.g.  range  curvature, 
range  alignment,  relative  scaling,  is  a  typical  example. 
Also,  motion  compensation  and  focusing  problems  adds  the 
complexity  to  the  system  analysis. 

System  classification  and  decomposition  principles 
proposed  in  Chapter  2  will  prove  very  useful  tools  in  the 
subsequent  chapters  in  obtaining  practical  and  efficient 
radar  image  reconstruction  procedures  under  the  various 
imaging  geometries. 


Chapter  4 


SAR  STRIPPING  MODI' 


4 . 1  Introduction 

Synthetic  aperture  radar  can  map  two-dimensional  areas 
with  high  azimuth  resolution  achieved  by  the  use  of  coherent 
signal  processing  upon  the  taraet  point  phase  histories 
induced  by  a  relative  motion  between  radar  and  taraet.  In 
imaging  a  large  sized  target  which  is  usually  stationary, 
e.g.,  the  ground  terrain,  the  radar  is  carried  on  a  moving 
vehicle,  e.g.,  an  airplane  or  spacecraft,  to  get  radar 
echoes  from  different  aspect  angles  of  the  target  points. 
Based  on  the  same  synthetic  aperture  princicle,  there  are 
basically  three  modes  of  SAR  ground  mapping  depending  on 
the  imaging  geometry:  spotlight  mapping,  strip  mapping,  and 
Doppler  beam-sharpening  [4-1,2],  This  chapter  concerns  the 
stripping  mode  of  SAR  terrain  imagine  which  is  the  most 
interesting  mode  having  obtained  great  attention  since  its 
birth  in  early  1950's. 

In  its  operation  an  airplane  or  spacecraft  flies 
straight  over  the  ground  of  interest,  radiates  pulses  at 
different  locations  and  records  the  returned  echoes.  As 
shown  in  Fig.  4-1,  let  (x,y,z)  be  a  rectangular  coordinate 
system  with  (x,y)  being  the  ground,  and  assume  unless  other- 
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wise  stated  that  the  flight  is  along  the  x  axis,  with  y=0 
and  at  a  constant  velocity. For  simplicity  a  side-looking 
radar  is  also  assumed,  although  it  could  as  well  be  squint 
[4-2].  In  the  geometry  depicted  in  Fig.  4-1,  we  use 
(Xi#yi,zi)  and  (x2,y2,z2)  to  denote  the  coordinates  of  an 
arbitrary  target  point  and  the  radar  receiver  on  the  air¬ 
craft,  respectively.  We  will  assume  that  the  effect  of  the 
height  of  mountains  and  structures  on  the  ground  are  negli¬ 
gible  and  that  the  curvature  of  the  earth  surface  can  be 
ignored  such  that  z^=0  for  all  the  ground  points  to  be 
mapped.  Furthermore  we  assume  that  the  time  origin  coin¬ 
cides  with  the  x  origin  of  the  aircraft.  Thus  x2=vt  where 
v  is  the  velocity  of  the  aircraft.  As  is  generally  the 
case,  the  antenna  is  assumed  to  be  shared  by  the  transmitter 
and  the  receiver.  This  necessitates  the  pulsed  nature  of 
the  signal  waveform  and  inevitably  creates  blind  ranges 
[4-3].  Echoes  from  targets  in  blind  ranges  reach  the  radar 
while  it  is  transmitting  and  not  receiving  and  thus  are 
lost.  Also,  the  pulse  repetition  frequency  (PRF)  sets  an 
upper  bound  to  the  maximum  range  without  range  ambiguities. 
For  simplicity  of  analysis  we  assume  that  the  signal  pulse 
train  consists  of  pulses  of  identical  waveform  at  a  constant 
PRF.  Thus  let  fs(t)  be  the  modulation  function  of  a  single 
pulse  centered  at  t=0  then  the  infinite  lengthed  pulse  train 
wave  function  will  be,  in  its  analytic  form, 
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00 


f(t)  =  ^2  fg  (t-nTs)exp  ( joJc  t )  (4-1) 

n=-» 

where  Tg  is  the  pulsing  period  and  <ii  the  angular  frequency 

of  the  carrier.  As  described  in  Fig.  4-2,  if  the  "effective" 

time  width  of  f  (t)  is  T  then  the  length  of  the  time 
s  p 

during  which  the  transmitter  is  not  in  use  between  consecu¬ 
tive  pulses  is  T  -T  which  decides  the  maximum  range 
s  p 

deviation  of  the  radar  returns  without  range  ambiguities. 

In  addition  to  above  timing  factors,  the  depression 
angle  ip,  or  the  angle  between  the  horizontal  plane  and  the 
radiated  beam,  and  the  antenna  pattern  are  major  parameters 
of  the  system  with  its  geometry  depicted  in  Figs.  4-3  and 
4-4.  The  point  spread  function  is  naturally  a  complicated 
expression  with  many  parameters  interacting  mutually  and 
hence  precise  image  reconstructions  demand  formidable 
efforts.  However,  depending  on  the  degree  of  accuracy  pur¬ 
sued,  if  we  model  the  PSF  of  the  system  in  some  desirable 
way  by  appropriate  geometrical  considerations  and  appro¬ 
ximations,  we  will  be  able  to  simplify  the  description  of 
the  system,  making  the  system  evaluation  work  the  purpose 
of  which  will  be  described  in  the  next  section  relatively 
easier  and  the  reconstruction  more  feasible.  Of  course, 
by  so  doing  we  inevitably  distort  the  system  by  using  an 
inexact  model,  and  an  incomplete  or  nonoptimal  (in  some 
sense)  reconstruction  of  the  target  image  is  to  be  expected. 
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Fig.  4-3.  Flight-path  geometry  on  vertical  plane 


Fig.  4-4. 


Flight-path  geometry  on  slant  range  plane 


It  is  obvious  that  the  more  approximations  one  makes,  the 
more  degradation  will  result  in  the  reconstructed  image. 

In  this  chapter  we  will  give  a  hierarchy  of  the  SAR  system 
models  in  progressive  inaccuracies.  We  will  show  that  these 
models  fit  the  system  classification  categories  in  Chapter  2 
nicely  so  that  the  reconstructions  can  be  done  readily  by 
using  the  appropriate  algorithms  summarized  in  Chapter  2. 

We  will  also  tabulate  the  approximations  made  and  their 
justifications  at  each  step.  It  will  be  shown  that  in  its 
simplest  form  the  system  is  separable  and  space-invariant. 

An  imaging  system  has  two  kinds  of  DOF's;  one  of  them 
is  computational  (DOFc)  as  studied  in  Chapter  2  and  the 
other  is  informative  (DOF) .  We  now  explain  the  latter  and 
differentiate  the  two. 

4 . 2  System  Evaluation 

In  general,  gathering  more  data  provides  more 
information  to  solve  for  the  unknowns  at  the  expense  of 
increased  requirements  in  storage  and  complexity  in  compu¬ 
tation.  On  the  other  hand,  intuition  suggests  that  after 
some  "threshold  amount"  of  data  is  obtained,  the  additional 
observations  do  not  always  provide  equal  amounts  of  new 
information.  This  is  due  to  the  inherent  "blurring"  of  the 
imaging  systems  and  observation  noise,  etc.  Thus  the 
concept  of  degrees  of  freedom  (informative  DOF)  has  arisen 
to  measure  the  number  of  truly  independent  samples  of 
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data  one  gathers  under  a  particular  imaging  system  [4-4,5]. 

In  the  system  evaluation  the  concept  of  eigenvalues  of 
a  correlation  matrix  or  the  Gramian  matrix  is  usually 
adopted  [4-6,4-8],  For  example,  in  the  continuous-discrete 
case,  we  equate  the  degrees  of  freedom  with  the  number  of 
eigenvalues  of  the  Gramian  matrix  larger  than  some  threshold 
determined  by  the  noise  level  of  the  system.  This  is  equi¬ 
valent  to  the  singular  value  analysis  of  the  system.  For 
the  purpose  which  will  be  clear  later  we  now  show  that 
orthogonal  transforming  the  input  and/or  output  data  of  a 
linear  system  will  not  change  its  eigenvalue  spectrum  and 
hence  its  system  performance,  other  than  a  possible  reorder¬ 
ing  of  its  singular  values.  Consider  the  discrete-discrete 
case  for  the  sake  of  ease  in  proof:  Let  [H]  be  the  matrix 
of  the  linear  system  and  [P] ,  [Q]  be  orthogonal  matrices 

multiplied  with  the  output  and  input  vectors,  respectively; 
PfcP  =  PPfc  =1  and  QtQ  =  QQt  =  I  where  the  superscript  t 
denotes  a  transposition.  The  set  of  eigenvalues  of 
PHQ (PHQ) =  PHQQtHtPt  =  PHHtPt  is  the  same  set  of  eigen¬ 
values  of  P^PHH*"  =  HHfc  except  for  additional  zeroes  due  to 
possible  size  difference  of  P  and  Q  [4-8].  Note  that  the 
sets  of  eigenvalues  of  PHQ  and  H  differ  in  general,  though. 

To  differentiate  the  DOF  and  DOF,  we  note  first  that 

c 

both  of  them  are  derived  from  the  system  PSF.  However,  the 
DOFc  is  usually  determined  before  DOF  and  it  need!  only 
the  structural  information  of  the  PSF  and  can  be  determined 
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by  identifying  the  PSF  with  the  simplest  ease  in  Chapter  2. 


The  DOFc  is  used  to  determine  the  amount  of  effort  required 
to  estimate  the  DOF  which  requires  the  actual  numeric  values 
of  the  PSF  in  the  matrix  diagonalization  manipulation. 


4.3  Derivation  of  Point  Spread  Function  for  the  SAR 
Imaging  System 

Referring  to  Fig.  4-1,  =  0,  X2  =  vt,  ^2  ~  0  and  z2 

is  the  flight  height.  Define  the  ground  range 

Rg  =  [(xrx2)2+(yry2,2]% 

=  [<xL-vt)2+y^]%  (4-2) 

and  the  slant  range 


=  [(xrx2,2+(yry2)2+(zrz2,2]J 

=  [(x1-vt)2+y2+z2]% 

- 


(4-3) 


The  propagation  delay  associated  with  a  point  target 

2  R 

at  (x.fy.fZ.)  with  range  R  is  —  where  the  factor  2  is 
lix  c 

because  of  the  round  trip  of  the  wave  propagation.  Let 
^Xl'yl^  reactivity  function  of  the  terrain  and 

A  (X1 ' yf» x2 ' z2^  k®  the  illuminating  intensity  of  the  antenna 
beam  on  the  terrain  point  (x^,y^)  when  the  receiver  is  at 
(X2*0,Z2).  If  the  antenna  pattern  remains  the  same  during 
the  flight,  it  is  easily  seen  that 
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A  (X1 '  Y}. '  x2 ' 2  2  ) 


A(x  -x2,y1,z2)  . 


The  received  echoes  are  the  product  of  the  illuminating 
pattern  A,  the  terrain  reflectivity  p  and  the  delayed 
signal  function  f,  summed  over  the  ground  coordinates 
(xl»yi),  where  the  propagation  decay  has  been  assumed 
uniform  in  the  (x^,y^)  plane  and  neglected 


z(t)  =  J  A(x^-x2,y^,z2)  p  (x^  ,y^)  f  ( t-^)  dx^dy^  (4-4) 

—  OO 

Substituting  Eq.  (4-1)  into  Eq.  (4-4), 

OO  OO 

z(t)  =  j  j  A(x1-x2,y1,z2)p(x1,y1)  ^  fs(t-^-nTg) 

—  OO  ^ 

exp|jwc(t-~)|dx1dy1  (4-5) 


If  we  interpret  Eq.  (4-5)  as  a  two  variable  to  one 
variable  system  with  p(x^,y^)  as  its  input  and  z(t)  as  its 
output,  it  is  obvious  that  the  system  is  linear  with 
point  spread  function 


OO 

h(t;x1,y1)  =  ^  A(x1-x2,y1,z2) fj t-nTg- 

n=-°°  ' 


2[(X1-X2)2+Yl  +  Z2] 


) 


exp(  jur 


[(x1-x2)2+y2+z^]2 


)l 


(4-6) 


OO 

=  exp{jw  t} 2  A(x  -nvTg,y  ,z  ) f (t- 
n=-°°  \ 


nTs“ 


2  (x1~nvTg) 2+y^+z2 


— ^exp  j-ji 


[(xrnvTs)2ty2+z2] 
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where  x_=nvT  is  a  discrete  radar  coordinate  and  it  has  been 
2  s 

assumed  that  during  the  transmission  and  receiving  of  a 
single  pulse  the  aircraft  is  approximately  stationary  so  that 
x2  is  substituted  by  nvTg  in  Sq.  (4-6).  This  is  valid  if 
(sufficient  conditions) 

a)  A(x1-nvTs,y1, z2 )«a (x1~nvTs+Ax2 ,y  , z^) 

b)  (|x1-nvTs|+Ax2)2-(x1-nvTs)2<<(y1+Ay1)2-y^  and 

c)  2  J^{  |  xx-nvTs  |  +Ax2)  2+y2  +  z2  j^-2  ^(x1-nvTg)  2+y2+z2J^<<Xc 

for  all  x-^  and  y^  under  illumination,  where  Ax2  is  the 
maximum  distance  the  aircraft  travelled  during  the  receiving 
of  a  single  pulse.  Ay1  is  the  range  resolution  desired. 

See  Fig.  4-5.  Note  that  Ax2<vTg  as  assumed  earlier.  The 
above  relations  assume  that  the  movement  of  the  radar  during 
the  receiving  of  a  single  pulse  is  so  small  that  a)  illu¬ 
mination  pattern  is  effectively  fixed,  b)  no  range  walking 
within  one  pulse  and  c)  induced  phase  error  is  much  less 
than  one  wavelength,  respectively. 

RX 

(a)  is  easily  satisfied  by  noting  that  Leff='L~  is  °f 
the  order  of  hundred  or  thousand  meters,  and  is 
therefore  greatly  larger  than  Ax2,  which  is  of 
the  same  order  of  the  azimuth  resolution  desired. 

(b)  (  |  x^nvTg  |  +Ax2)  2- (x1-nvTg)  2 

2 

=  2Ax2 I x^-nvTg | +Ax2 

£  2Ax2  (Lg ff+Ax2 )  «s  2Ax2Lg££ 

while 
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(yl+Ayl)2_yl  =  2Ayi‘y+Ayl  85  2Ayx*y 

Since  AX2  is  of  the  same  order  as  Ay^ (azimuth  resolution 
vs.  range  resolution) ,  (b)  will  be  valid  if 


Jef  f 


<<  Y, 


(4-7) 


(c)  2 


|  £(  | x1~nvTg |  +A^ )  2+y^+z^  -  [(Xj-nvT,,)  2+y2+z2] 
2,_2,%,1,(IVnyrsl+Ax2) 2  _ ,,22.%,  (xl'nVTs) 


2  j  (yi+z2)  +%- 


,  2  2,% 

y!+Z2 


(yi+z2}  -t—TTTY 

(yl+z2> 


2  2 
=  — a — T-rr  ( I  x.  -nvT  j  • Ax9+Ax~) 
(y2+z2)*  i  s  z  2 

2 1 Xj-nvTg |  •  AX2 

(yI+4^ 

2Leff*Ax2 

1  -Trr&t <<  Ax2 

yl+Z2  ^ 


where  validity  of  Eq.  (4-7)  has  been  assumed.  Thus  if  we  let 
Ax2  be  of  the  same  order  as  Xc,  (c)  will  be  satisfied.  In 
fact,  the  achievable  r  :uth  resolution  is  of  the  same  order 
as  Xc  [4-11] . 

Because  the  sinusoidal  phase  term  exp{jwct}  in  Eq.(4-6) 
does  not  carry  any  information  on  p(x^,y^),  a  preprocessing 
can  shift  it  to  any  lower  frequency  desired.  In  optical 

processing  upon  SAR  data,  the  "offset"  frequency  u>Q^0  is  to 
separate  reconstructed  twin  images  from  each  other  and  from 
other  useless  images  [4-9,10]. 

Thus , 
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h (t;x1 ,y  ) 


00 


exp{jw  t}  A(x,-nvT  (y,,z2) 
n=~°° 

expi-jv2[(^~n^’M+Z^| 


(4-8) 


l  C  1 

Although  the  return  of  the  pulse  train  from  the  two 
dimensional  target  field  is  one  dimensional  -  i.e.,  function 
of  t  only,  the  recording  of  data  is  usually  two  dimensional. 
For  example,  because  of  its  huge  capability  for  2-D  data 
storage,  film  has  been  widely  used  for  SAR  data  recording. 
This  is  conveniently  done  for  latter  processing  because  the 
signal  returns  from  different  pulses  do  not  overlap  in  time. 
The  returns  of  each  transmitted  pulse  are  arranged  side  by 
side  as  in  Fig.  4-6 (b)  in  which  the  data  lines  are  contained 
between  0  and  Tg  in  the  new  t  axis  and  extends  along  the 
flight  path  axis  x  which  is  perpendicular  to  the  t  axis. 
Equivalently,  the  original  one  dimensional  t  axis  denoting 
the  data  collecting  time  is  transformed  into  two  dimensions, 
discrete  X2  and  continuous  t.  Because  the  coordinate  trans¬ 
formation  above  is  a  reordering  of  the  data  which  is  an 
orthogonal  transformation,  nothing  has  been  changed  from  a 
DOF  point  of  view  from  the  analysis  in  section  4.2. 

Now  S(x2»t)  is  nonzero  only  for  0  <  t  <  Tg  and  X2  is  a 
discrete  variable  occuring  at  nvTg  only.  Note  that  X2  and 
t,  which  are  the  variables  in  the  data  domain,  have  dimen- 
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J — — I — vk<U — L 


•  • 


(a)  received  signal 


S(x2,t) 


(b)  reordering  of  (a) 


Fig.  4-6.  Reordering  of  received  signal 


sions  of  length  and  time,  respectively. 


The  PSF  expressed  in  (x2,t)  variables  now  becomes 
h(x2,t;x1,y1)  =  exp  ( t+nTg )  J  A  (x^x.,  ,  y1 ,  z  ,  ) 

f  /t  2[(xrx2)2+yi+z2]%\ 

'  t-  ' 


exp  <- jcj 


[(x1-x2)2+y2+z2]i 


(4-9) 


where  0  £  t  £  Tg,  x2  =  nvTs‘ 

Eq.  (4-9)  can  easily  be  identified  with  the  kernel  of 

a  two-dimensional  transformation,^,  v,  Tg,  z2,  c,  A 

and  f  being  parameters.  Sometimes  it  is  desirable  that  the 

data  obtained  from  the  system  kernel  Eq .  (4-9)  be  multiplied 

x2 

with  exp  (- ja)QnTg)  =  exp  (- j  wQqj-)  ,  or  exp(-jw0t)  yielding 
range  offset  case  and  azimuth  offset  case,  respectively 
[4-12].  Although  the  multiplication  is  not  necessarily  an 
orthogonal  transformation  and  thus  would  affect  the  true 
DOF  of  the  imaging  system,  the  exact  effect  on  the  DOF  will 
not  be  investigated  here.  As  will  be  explained  in  Chapter  7 
this  is  one  of  the  preprocessors  designed  to  make  subsequent 
processings  or  presentation  easier.  The  offset  data  has 
equivalent  kernels 

h(x2,t;x1,y1)  =  exp  (jojQt)A(x^-x2,y^,z2) 


LI 

2  2  2 
(x1-x2)  +y^+z2 

IS  ( 

i  2 

2  2  2 
(X1~X2)  +yl+Z2 

V 

c 

r p! 

)  3a)c 

c 

i 

range  offset  case  (4-10) 
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/  wo  \ 

h(x2,t;x1,y1)  -  exp  (  j— x2  1 A  (x^-x2  ,  ,  z2  ) 

/  2[<xrx2)2+yi+zIP\  j  .  2 [ow 

£s^t - jeXpj-:»c - j- 


azimuth  offset  case 


(4-11) 


In  the  following  analysis  we  shall  assume  that  the  range 
offset  case  is  used. 

4 . 4  Simplification  of  PSF 

In  this  section  we  try  to  simplify  the  PSF  (Eq.  (4-10)) 
of  the  SAR  imaging  system  from  several  physical  considera¬ 
tions.  Along  with  the  simplifications,  a  hierarchy  of 
models  of  PSF's  with  decreasing  complexity  will  be  derived 
along  with  their  associated  assumptions  and  approximations. 

We  start  by  noting  the  strong  relations  between 
variables  x1  and  x2 ,  y1  and  t,  respectively:  in  Eq .  (4-10), 

the  argument  of  fg,  t- j^2  |\x^-x2)  2+y2+z2J^J^/ c,  were  it  not 

2 

for  the  factor  (x^-x2)  ,  would  yield  a  simple  one  dimen¬ 

sional  relation  which  connects  y^  with  t  to  provide  the 
range  information  independent  of  azimuth  modulation.  This 
is  the  only  way  which  makes  h  (x2 , t ; x^ , y^ )  nonseparate  in 
azimuth  and  range  in  Eq .  (4-10).  If  this  fact  can  be 

ignored,  e.g.,  if  the  propagation  delay  induced  by  the 
variation  in  (x^-x2)  is  much  smaller  than  the  range  resolu¬ 
tion  interested,  then  the  PSF  can  be  considered  separate  in 
azimuth  (from  x^  to  x2)  and  range  (from  y^  to  t) : 
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Physically  this  means  that  the  range  resolution  cells 
under  antenna  illumination  do  not  move  to  overlap  each 
other  as  the  flight  continues.  The  situation  is  depicted 
in  Fig.  4-7,  where  the  range  y  of  a  target  point  is  plotted 
as  a  function  of  the  azimuth  distance  xy-x2  of  that  point. 
Various  ways  have  been  proposed  to  alleviate  the  problem 
of  range-azimuth  coupling  [4-13,14,15]  where  the  range 
curvature  is  not  to  be  ignored.  However,  if  pr  is  the 
range  resolution  pursued  and  8  the  effective  beam  width 
then  this  range  curvature  will  be  negligible  if 


|  ^(x1-x2)  2+v‘i+z^\}%-  (yi+z2)^l  1  Pr 


for  all  x2.~x2  an<^  ^1  by  the  radar.  Assuming 

2  2  > 

^X1~X2^  <<  ^l+z2^  2  as  ;’'n  Eci'  anc1  8  small,  we  have 

[(x1-x2)2+yJ+Z2^-(y2+z2)%  (y2+z2)-V%_|_|_i-  i?\+z\)h 


(xrx2} 


(y  l+zp 


.  2 ,  2»: 
(yl+Z2 


*  f  Leff 
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So  we  require  that 


pr  >  8  Leff 


(4-13) 


for  negligible  range  curvatures.  I.e.,  PSF  can  be 
approximated  by  Eq.  (4-12) ,  which  is  separate  in  azimuth 
and  range,  if  the  inequality  (4-13)  holds. 

To  see  that  the  system  with  kernel  Eq.  (4-12)  is 
separate  in  azimuth  and  range,  we  rewrite  Eq.  (4-12)  as 


2  . 


h(x2,t;x1,y1)  =  exp ( jwQt) f g 


A(x1~x2  ,y1  ,z2)  exp<-ja) 


2  [(xrx2)2+yi+z2]% 


! 


=  ht(t;y1)hz(x2;x1,y1) 


(4-14) 


(t  _  ,  2  2 .  i: . 

2 (y1+z2 ) -\ 

t - - -  I  is  the  response 

of  a  unit  point  source  (impulse  response)  at  y^,  which  is 
independent  of  azimuth  dimension;  and 


hz(x2;xl'yl)  =  A(xi~x2'yl'z2)exp)  ^ 


i'3' 


_(Xl-x2)2+yJ+^] 


‘I 


is  the  azimuth  response  of  a  unit  point  source  at  (x^,y^). 
Note  that  hz (x2;x1,y1)  =  hz<xi"x2;yi)  is  sPace  invariant  in 
and  x2,  but  varies  its  form  as  y^  changes.  If  we  use 
Eq.  (4-14)  as  the  kernel  of  the  SAR,  then  the  input-output 
relation  will  be 
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z(x2,t)  =||  ht (t;y1)hz (x1-x2;y1)p (x1,y1)dx1dy1 

—  oo 

OO  00 

=  hfc(t;y1)  h2  (x1-x2;y1)p  (x1,y1)dx1"ldy1 

—  OO  —OO  -* 

(4-15) 

Thus  equivalently  the  imaging  system  of  SAR  is  a  transforma¬ 
tion  of  p(x^,y^)  into  z(x2,t)  in  a  sequential  order: 
azimuth  transformation  followed  by  range  transformation. 
Because  of  the  dependence  of  hz  upon  y^ ,  h  and  hz  are  not 
separable  and  thus  their  order  cannot  be  interchanged  in 
modelling  the  system.  Accordingly,  the  reconstruction  of 
the  ground  reflectivity  function  pfx^.y^)  from  its  image 
z(x2,t)  has  to  follow  the  reversed  order. 

To  put  Eq.  (4-14)  in  a  more  practical  form,  we  note 

2  2 

from  inequality  (4-7)  that  |x^-x2|  <<  y1+z2  and  hence 


exp  <-]  oj 


[(xr*2)2+yi+z2? 


exp<-30) 


2  2  ^  <xi-x2) 

2(yl+z2>  +7T7T3 

(y1+z2) 


=  exp 


-Dm 


7  2  3- 

2(y^+z2) 2 


exp 


-D<*> 


(x1~x2) 

c(y\+z\)h 


(4-16) 


and  thus  hz (x2; ,y^)  can  be  approximated  by 

,  -  ,  2  2 .  h 

2 (yi+z2 )  # 

hz(x2;xl'yl)  155  exP]"3a)c - |A  (x^-x2  ,  Y  x  ,  z  2 ) 


(x1-x2) 

eXPl'j^  ^2-0 


(yl+z2) 


(4-17) 
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Because  the  phase  term  exp<-ju) 


~  ,  2  2.  % 
2 (y1+z2)^ 


is  independent  of 


c  c 
v  / 

azimuth  its  effect  can  be  taken  out  and  absorbed  in  the 
reflectivity  function  p(x,y).  In  terms  of  the  system 
block  diagrams,  we  now  have  Fig.  4-8  where 


ht  (t;Yl)  =  exP  < j“0t)  f< 


hz(xl_X2'yl)  =  A(x1“x2'yl'Z2: 


(4-18) 


(  .wc  (x 
expj-3-  _ 


rx2>  ) 


Eq.  (4-18)  is  the  form  assumed  for  most  SAR  processing. 

The  equation  clearly  expressed  the  separability  of  the  PSF 

in  which  the  range  information  is  provided  by  the  timing 

of  the  returns  whereby  the  azimuth  information  is  provided 

by  the  phase  history  induced  by  the  motion  x^-x2. 

All  the  assumptions  required  so  far  can  usually  be 

justified  in  practical  SAR  systems  of  stripping  mode.  We 

proceed  to  approximate  Eq .  (4-18)  by  a  PSF  separable  in 

azimuth  and  range,  i.e.,  such  that  h  is  independent  of  y. 

z  X 

in  Eq.  (4-18) .  This  will  be  true  if 
(A)  A(x1-x2,y1,z2)  =  A(x1~x2,z2) 


(B) 


( x^ — x  2 ) 


(yT+Zo)  2 


v 


lxrx2)2  ) 

(Y'l+z22)h  ) 


where  y^  and  yj^  are  azimuth  coordinates  for  any  two  target 
points  under  the  antenna  illumination.  (A)  can  be  made 
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approximately  true  by  an  appropriate  antenna  pattern  tail¬ 
ored  for  the  imaging  geometry. 


.  V  +2 

(B)  will  be  true  if  3 lmax  2  «  1  where  y.  and  y.  are 

- «-  3 1mm  3  lmax 

yf  .  +zZ 
3 lmin  2 

the  y^  coordinates  of  the  target  points  at  maximum  and 
minimum  ranges  covered  by  the  antenna  beam,  respectively. 

In  that  case,  the  azimuth  modulation  would  be  same  linear 
FM  at  all  range  bins.  We  shall  not  elaborate  on  the  exact 
requirements  for  the  inequality. 

It  is  pointed  out  that  in  general  (B)  cannot  hold. 
However,  if  it  could,  then  the  system  of  Eq.  (4-18)  would 
be  separable: 

h  (x2» t;x^,y^)  =  hfc  (t.-y^)  hz  (x2~x1)  (4-19) 


Further  theoretical  reduction  of  the  PSF  is  still 

possible:  if  we  make  the  offset  frequency  0  and  assume 

2  2  ^ 

that  y.  >>  z_  for  all  valid  y.  such  that  z(y1+z2)2/c  can 

2y  1  A  0 

be  approximated  by  — — ,  then  by  changing  variable  ty  = 
we  will  have 

h 


A  2yl 


or 


so 


ht(t;t2)  =  fg(t-t2) 


h (x0, t ; # y ^ )  =  h1  (x^/trx,  ,t1  ) 


1  '2  '  '  i f  1 


=  hfc  (t-t2)  hz  (Xj-Xj^)  (4-20) 


89 


which  is  a  separable  space  invariant  PSF  (SSIPSF) . 

A  summary  of  the  properties  of  the  PSF's  under  various 
assumptions  is  listed  in  Table  4-1. 

4 . 5  Conclusions 

This  chapter  presents  the  analysis  of  the  stripping 
mode  of  the  SAR  from  a  PSF  point  of  view.  A  series  of  PSF's 
of  the  imaging  system  is  derived  based  on  various  geometri¬ 
cal  and  mathematical  assumptions.  Depending  on  the  degree 
of  precision  one  pursues  the  PSF's  vary  from  the  most 
complicated  case  of  a  non-separate  one-dimensional  invariant 
processing  to  a  separable  two-dimensional  invariant  proce¬ 
ssing.  The  requirements  for  most  of  the  approximations  are 
derived.  Many  parameters  are  present  in  the  radar  equations 
which  look  formidably  complicated  in  some  of  its  most 
precise  forms.  By  classifying  the  equations  according  to 
the  method  in  Chapter  2,  parameters  can  be  clearly  separated 
from  the  variables  associated  with  the  imaging  equations  and 
great  insight  into  the  image  formation  methods  is  readily 
provided . 

The  concept  of  preprocessing  is  briefly  mentioned  in 
this  chapter.  Details  of  processing  techniques  for  the 
reconstruction  are  not  presented  until  later  chapters. 


91 


Chapter  5 

TURNTABLE  RADAR  IMAGING 


5 .  I  Introduction 

Synthetic  aperture  radars  achieve  high  azimuth  resolu¬ 
tion  by  coherently  processing  the  phase  histories  of  the 
target  points.  Based  on  this  principle  there  are  three 
modes  of  SAR:  the  stripping  mode,  the  doppler  beam 
sharpening  mode  and  the  spotlight  mode.  In  the  previous 
chapter  the  stripping  mode  has  been  analyzed  and  many  radar 
parameters  introduced.  In  this  chpater  we  study  a  radar 
imaging  geometry  closely  resembling  the  spotlight  mode  and 
in  a  well  controlled  environment  where  the  relative  motion 
between  the  radar  and  target  is  a  strict  circle. 

We  will  show  that  the  underlying  radar  imaging  system 
is  very  similar  to  a  computer  aided  tomographic  (CAT)  system 
from  the  PSF  point  of  view.  Unfortunately,  physical 
limitations,  e.g.  the  aspect-angle-dependence  of  the  target 
reflectivities  and  the  shadowing  effect  from  3-D  obscura¬ 
tion,  discourage  one  from  applying  a  tomography-like  recon¬ 
struction  to  the  reflected  signals.  Hence,  instead  the  SAR 
principles  will  be  applied  directly  to  small  angle  looks 
and  several  looks  will  then  be  registered  and  incoherently 
summed  to  give  the  full  reconstruction  of  the  object 


reflectivity  function.  A  DOF  as  well  as  Nyquist  rate 
analysis  in  the  frequency  domain  will  be  derived  to  show 
the  minimum  number  of  data  points  required  for  image  recon¬ 
struction  under  specified  physical  constraints  and  require¬ 
ments.  Basic  relations  between  bandwidth  and  resolution 
will  also  be  discussed. 

5 . 2  The  RAT  SCAT  Facility  and  Data  Acquisition 

The  original  design  purpose  of  the  RAT  SCAT  (standing 
for  RAdar  Target  SCATter  Site)  facility  is  to  measure 
radar  cross  section  (RCS)  of  various  targets  at  different 
distances,  angles,  elevations  and  frequencies  [5-1].  In 
operation,  the  target  (say  a  model  airplane)  is  placed  on  a 
rotator  at  a  distance  rQ  from  the  radar  to  its  rotation 
center  as  shown  in  Fig.  5-1.  A  reference  sphere  S  is 
sitting  at  distances  r^  from  the  radar  R  and  ^  from  the 
rotation  center  C.  The  angle  between  line  RS  and  the  target 
line  of  sight  RC  is  a.  Let  (£,ri)  ,  (x,y)  be  two  rectangular 
coordinate  systems  with  origins  at  C.  Let  (£,n)  be  fixed 
on  the  target  and  (x,y)  on  the  ground  at  an  angle  0  from  the 
former  coordinates,  as  depicted  in  Fig.  5-2.  At  discrete 
angle  9^  the  radar  radiates  continuous  waves  (CW's)  single 
frequencies  f^  for  various  k.  The  radar  receiver  at  the 
reference  sphere  S  takes  the  signal  directly  from  R  to  S  as 
a  reference  and  beats  it  with  the  signal  reflected  from  the 
target.  The  resultant  in-phase  and  quadrature  phase 


components  become  the  data.  This  process  continues  for 
different  k  and  i  to  form  a  2-D  data  array.  For  simplicity 
we  shall  assume  that  at  each  aspect  angle  the  radar 
radiates  the  same  set  of  M  frequency  waves  at  a  fixed 
frequency  step  Af.  We  shall  assume  that  the  angle  step  A0 
is  also  constant  as  one  advances  the  aspect  orientation. 

5 . 3  Hypothetical  Target  Reflectivity  Function 

Referring  to  Fig.  5-2,  let  f(£,n)  be  the  reflectivity 
function  of  the  target,  where  by  reflectivity  function  f(£,n) 
we  mean  the  ratio  of  the  received  signal  due  to  a  point 
target  at  (£,n)  with  the  radiating  signal.  At  wavelengths 
X  small  compared  with  the  curvature  of  the  target  body,  as 
used  in  our  experiments,  the  target  looks  specular  to  the 
radar  [5-2]  so  that  only  those  surfaces  at  appropriate 
orientations  to  the  radiation  path  reflect  strong  energy 
back  to  the  radar  receiver.  In  addition,  wherever  a  point 
(£,n)  is  blocked  by  some  other  points  or  surfaces  in  the 
line  of  sight  (LOS)  to  the  radar,  shadowing  occurs.  In 
other  words,  the  shadowing  effect  occurs  because  of  the 
non-convexity  of  the  surface  of  the  target.  Thus  f(£,n)  is 
actually  a  function  of  aspect  angle  9.  Nevertheless,  for 
ease  of  analysis  we  assume  that  f(£,n)  is  independent  of  0 
and  we  shall  see  a  close  resemblance  of  this  imaging  system 
to  that  of  a  tomography.  A  great  deal  of  insight  can  be 
obtained  by  this  theoretical  assumption.  Even  if  we 


release  this  assumption,  as  we  shall  do  later,  the  DOF 
analysis  based  on  fixed  f(£,n)  is  still  valid  in  the  real 
situation . 


5 . 4  PSF  of  the  Imaging  System 

In  Fig.  5-1  let's  assume  that  the  angles  a  and  3  are  so 
small  (in  our  experiments  they  are  «  1°)  that  the  signal 
reflected  to  S  would  essentially  be  the  same  as  if  S  were  in 
the  line  RC  at  the  same  distance  from  C.  The  distance 
between  the  object  point  at  (£,n)  and  the  radar  is  a 
function  of  rQ ,  0,  £  and  n . 


r  (rn , 0 , £  #  n )  =  [ (rn-x) 2+y2]^ 


-  (rQ-x), 


(5-1) 


where  we  have  assumed  that  rQ-x  >>  y  for  all  target  point 
coordinates  (x,y)  or  (£,n)  and 


x  =  £cos0+nsin9 


£  =  xcos9-ysin9 


y  =  ncos9-£sin0 


n  =  ycosO+xsinB  (5-2) 


where  0  is  the  angle  from  £  axis  (on  the  tarqet)  to  the 
x  axis  (on  the  imaging  device) .  Thus 

3r(rg,0,£,n)  9  ( (r^-CcosS-nsinB ) 2+ (ncos9-£sin9 ) 2 ] ^ 


2  (rQ-x)  (£sin0-ncos9)+2y (-nsin9-£cos9) 


r  (r q  ,  0  ,  £ ,  n ) 


=  £sin9-ncos6  =  -y 


(5-3) 


r>j 


where  we  have  assumed  Tq-x  «  r (rQ , 8 , £ , n  )  . 

Equations  (5-1)  and  (5-3)  reveal  that  the  lines  of 
constant  range  and  lines  of  constant  doppler  are  parallel 
to  the  y,x  axis,  respectively.  Lines  of  constant  range  can 
be  separated  or  resolved  by  range  or  timing  processing  while 
lines  of  constant  dopplers  can  be  separated  by  azimuth  or 
doppler  processing,  roughly  speaking. 

We  proceed  to  find  the  PSF's  of  the  imaging  system. 

Let  a  be  the  maximal  radial  extent  of  the  target. 

Define  line  projection  at  angle  0  and  range  x  by 


g( 


C 

x,0)  =  | 


f (C,n)dy 


(5-4) 


-a 


which  is  the  integration  of  the  reflectivity  function  f(£,n) 
over  all  target  points  at  distance  rQ-x  to  the  radar  when 
the  x  axis  is  at  an  angle  0  from  £  axis.  The  radar 
transmits  signals  of  the  forms 


f^(t)  =  Acos  (2Ttf^t+$)  k  =  1 , 2  ,  .  .  .  ,  M  (5-5) 


at  discrete  0^,  i  =  1,2,...,N,  where  A  is  an  amplitude 
factor  and  <I>  a  phase  term.  The  signal  reflected  to  the 
reference  sphere  from  the  "line  mass"  at  range  r^-x  is 


z(x,i,k,t)  =  Bg(x 


,e.)COs[2TTfk(t--^-! - y*] 


(5-6) 


where  the  signal  power,  propagation  decay  and  reflectivity 
phase  have  been  absorbed  into  the  complex  constant  B.  The 
return  from  the  whole  object  is  a  summation  of  the  contri- 
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butions  from  different  ranges: 


z(i,k,t)  =  I  z(x,i,k,t)dx 
-a 


=  B  |  g  (x 

-a 


,6i)cos  27ifk^t 


(5-7) 


ax 


which  is  again  a  CW. 

The  signal  received  by  S  along  path  RS  assumes  the 

form 


y(k,t)  =  C  cos 


•5 


(5-8) 


When  signals  in  Eqs .  (5-7)  and  (5-8)  are  mixed  or  multiplied 

at  S ,  they  become 


z  (i,k,t)y(k,t)  =  BC 


BC 


a  T  (  r0+r2-2x\  1 

g  (x ,  9  ^  )  cos  |^2 tt f k  ft - - - J  + 

cos  |^2tt f k ^t 

a  \  (  r0+rl+r2'2x\ 

-a  1  L  V  k  '  ) 

•Vcos[2"£k(^T^+¥)]dx 


dx 

+2$  dx+ 


BC  fa  , 

T  g(x 

1  -a 

(5-9) 

The  first  term  in  Eq.  (5-9)  is  a  high  temporal  frequency 
term  which  can  be  filtered  out  to  leave  only 


I  (i  ,  k)  =  D 


J  _ 


g  (x 


,  9^  )  cos 


,  c  2x,.,  -  rl  1 0  ~ "  2 
2u  f,  — +2tt  t,  — — — 
k  c  k  c 


^  dx 


(5-10) 


Similarly  by  beating  z(i,k,t)  in  Eq.  (5-7)  with  90  shifted 
version  of  y(k,t),  one  obtains 

r. -r  -r . 


Q  (i , k)  =  -D 


fa  ,  „  ,  .  c  2x,n  £  rl'r0  2\ 

g(x,Gi)sin^27Tfk-E-+27rfk - - - J 


dx  (5-11) 
98 


Define 


?k(g(x,ei)>  =  g{x 


.  (  27rfkx\ 

, 0 i ) exp  I - 3  -  j 


and  compare  Eqs.  (5-10)  and  (5-11)  with  (5-12).  2tt  f , 


(5-12] 


rl'r0_r2 


in  (5-10)  and  (5-11)  is  a  linear  phase  term  with  resoect 

to  f  which  if  compensated  for  would  make  I(.i,k)  and  Q(i,k) 

the  Fourier  component  of  the  shadow  gram  g(x,0.)  at  angle 

2  f  1 

0^  and  frequency  — .  According  to  the  circular  projection 

theorem  [5-5],  the  data  [I (i , k) , Q ( i , k) ]  takes  on  the  value 

in  the  ring  in  Fourier  frequency  domain  as  shown  in  Fig.  5-3 

with  the  extra  linear  phases  mentioned. 

Mathematically  the  structure  of  the  data  obtained  has 
a  considerable  similarity  to  that  of  the  Fourier  transformed 
data  of  an  infinitely  narrow  parallel-beam  tomographic 
imaging  system  [5-4,5]: 

In  the  computer  aided  tomographic  system,  let  p(y,6) 
denote  the  integral  of  the  object  transmission  (in  the  radar 
system  it  would  be  the  reflectivity)  function  f(£,g)  or 
f(x,y)  along  the  line  L(r,0),  as  shown  in  Fig.  5-4.  As 
before,  coordinates  (£,n)  are  assumed  fixed  on  the  object 
or  target  and  (x,y)  fixed  on  the  data  gathering  device  (in 
the  radar  system  it  would  be  the  transmitter  and  receiver) . 
The  two  coordinate  systems  are  related  by 


x  =  £cos0+risin9 
y  =  ncos0-£sin8 


£  =  xcos9-ysin6 
n  =  ycos0+xsin8 


(5-13) 


N  angle  steps 


as  in  Eq.  (5-2).  The  equation  for  the  line  L(r,9)  is 
x-r=0  or  £cos0+nsin0-r=O .  Thus 


P(r,9)  =  J  f(C,n)ds 
L(r,0) 


(5-14) 


where  s  is  along  line  L(r,0) 


The  set  S={p  (r,  0  )  | -°°<r<°° ,  O<0<ir)  is  called  the  shadow- 
gram  which  is  a  transformation  of  the  object  function  f(£,n) 
or  f(x,y),  and  is  the  data  one  obtains  from  a  parallel 
beam  tomographic  system.  A  famous  way  to  reconstruct  f(£,n) 
from  p(r,0)  is  to  use  the  back  projection  convolutional 
method  [5-6].  Other  reconstruction  algorithms  can  be  found 
in  [5-7] .  However,  if  we  take  the  Fourier  transform  of 
p(r,0)  with  respect  to  r, 

oo  oo 

P(u>,0)  =  |  p  (r ,  0  )  exp  (- jwr)  dr  =  j  j  f  ( £  ,  n )  exp  (- jwr )  dsdr , 


-°°  L  ( r ,  9  ) 


and  compare  with 


F  (oj ,  0 )  =  f  (?  ,  n)  exp  [- jw  (^cos0+nsin0  )  ]  d£dn  (5-15) 


where  F(w,0)  is  the  two-dimensional  Fourier  transform  of  f 
in  polar  coordinates,  it  is  obvious  that 


P  (w ,  9 )  =  F(u>,6) 


Hence  f(^,n)  can  be  reconstructed  by  taking  a  2-D  inverse 
FFT  upon  F(oj,9)  or  P(w,9): 


;,n)  =  -iy 

4  7T 


d0  P  (oj,  0 )  exp  [  jo)  (^cos9+nsin0  )  udu]  (5-16) 


0  - 00 


Thus  it  seems  the  RAT-SCAT  image  could  be  obtained  by 
a  2-dimensional  inverse  Fourier  transformation  upon  the 
data,  just  as  in  the  transformed  tomographic  system.  There 
are,  however,  several  differences  which  should  not  be  over¬ 
looked  : 

(1)  The  validity  of  the  hypotheses  made  in  Section 

5.3. 

(2)  A  linear  phase  factor  (if  not  compensated  for) 
in  RAT  SCAT  as  in  Eq.  (5-11). 

(3)  Bandpass  (instead  of  low  pass)  property  of  data 
gathering  ability  of  RAT  SCAT  as  in  Fig.  5-2. 

Because  of  (1)  and  (2)  one  has  to  seek  algorithms 
other  than  those  from  tomography  to  reconstruct  our  RAT  SCAT 
images.  Becuase  of  (3)  the  radar  images  have  no  "DC"  term 
and  magnitudes  or  intensity  of  complex  reflectivities  have 
to  be  taken  in  order  to  display  the  final  stages.  We'll 
next  derive  the  minimal  number  of  data  required  to  recon¬ 
struct  a  target  of  given  site  and  required  resolution,  show 
the  relation  between  resolution  and  bandwidth  and  then 
propose  a  reconstruction  algorithm  which  turns  out  basic  to 
all  modes  of  SAR's. 

5.5  DOF  of  the  Data 

The  DOF  problem  of  an  imaging  system  arises  because  the 
data  available  for  image  reconstruction  are  not  independent. 
In  other  words,  there  are  overlappings  among  different  PSF's 
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of  the  data  acquisition  system.  A  common  method  to  tackle 
the  DOF  problem  is  to  form  the  correlation  matrix  [T]  of  the 
imaging  system  and  try  to  diagonalize  it.  The  eigenvalue 
spectrum  and  the  rank  of  [T]  give  a  quantitative  measure 
of  the  DOF. 

The  degrees  of  freedom  problem  of  a  tomographic  projec¬ 
tion  system  has  been  analyzed  in  detail  by  McCaughey  and 
Andrews  [5-5] .  Because  of  its  mathematical  similarity,  we 
will  apply  some  of  their  results  to  the  RAT  SCAT  data. 

First,  for  the  purpose  of  clearity  and  convenience, 
let's  normalize  the  spatial  extent  of  the  target  as  in  [5-5]. 
We  assume  that  the  object  has  a  maximal  radial  extent  of  a 
meters  and  that  the  minimal  and  maximum  radar  frequencies 
applied  at  each  aspect  angle  are  f  .  and  f  ,  respectively, 
which  give  respective  2  way  wavelengths 
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(5-17) 


where  c  is  the  speed  of  EM  waves.  The  normalized  frequen¬ 
cies  have  a  spectrum  from 


fl  X 

max 

Refer  to  Fig.  5-5. 
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For  the  circular  imaging  system  on  the  unit-radius 
object  it  has  been  shown  in  [5-5]  that  the  (i,m,k,l)th 
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(a)  spatial  domain 


(b)  frequency  domain 


Fig.  5-5.  Spatial  and  frequency  domains  of  a 
spatially  normalized  target 


entry  r*'™  of  [ r ]  is 


i,m 
r  = 


exp[-j2mj^(£)cx>s0^+nsin9^)  ]exp[+j2Tru^  (£cos8m+nsin9m)  ]dE)dn 


J-j_  [2tt  (u^— 2u^u^cos  (0i"em)+ui5  3^ 
tu^-2uku1c°s(0i-0m)+u^“ 


(5-19) 


where  (i,k)  and  (m,l)  are  indexes  of  any  two,  possibly 
same,  data  points;  i  and  m  are  azimuth  angle  indices  and  k 
and  1  are  frequency  indices  along  radial  direction.  J^t*) 
is  the  first  order  Bessel  function  of  the  first  kind,  R  is 
the  unit  circle  and  O<0  .  ,  0  <2r,  f  <u,  ,u.<f0.  Ea.  (5-19) 
can  be  rewritten  as 


•  J  [ 2 it p  (k,  1 ,  i  ,m)  ] 

ri,m  _  1 _ 

k ,  1  p  (k ,  1 ,  i  ,m) 


(5-20) 


2  2  J- 

where  p(k,l,i,m)  =  [Uj^^u^u-^cos  (0^-9^  h-u^]  2  ,  which  is 

nothing  but  the  distance  between  two  points  (i,k)  and  (m,l) 

in  the  frequency  plane.  In  other  words,  r*'™  is  only  a 

function  of  the  distance  between  the  two  data  points. 

Therefore,  one  can  think  of  the  data  as  samples  of  a  two- 

dimensional  stationary  process  with  correlation  function 
J,  (2rp ) 

— i -  which  is  the  Fourier  transform  of  the  unit  circle 

P 

which  defines  the  object  spatial  extent.  The  correlation 
pattern  is  the  same  no  matter  where  the  points  are  consi¬ 
dered  and  manipulation  methods  of  stationary  processes  such 
as  the  Nyquist  sampling  theorem  can  be  applied  to  easily 
solve  the  DOF  problem. 
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‘"1  ■  1  ■  • ' 

A  stationary  process  with  correlation  function  — - - 

P 

has  a  power  spectrum  (Fourier  transform  of  the  correlation 
function)  circ(r)  and  is  bandlimited  with  radial  cutoff 
frequency  1,  as  shown  in  Fig.  5-6.  According  to  the 
Nyquist  sampling  theorem,  the  sampling  interval  in  the 
frequency  domain  must  be  <  1/2  in  order  to  avoid  any  alia¬ 
sing  in  the  spatial  domain.  This  coincides  with  the  more 
intuitive  fact  that  to  adequately  represent  an  object  of 
maximum  spatial  extent  S,  the  sampling  interval  in  its 
frequency  domain  should  not  be  larger  than  1/S  (note  that 
S=2  in  our  case  of  a  normalized  target) .  This  information 
can  be  used  to  determine  the  maximum  allowable  frequency 
step  in  the  radial  direction  and  maximum  ancle  step  in 
azimuth  without  ambiguities.  As  an  alternative  way  of 
analysis  tne  "effective  correlation  region"  in  the  data 
domain  can  be  determined  from  Eq.  (5-20)  and  approximate  DOF 
of  the  data  band  in  Fig.  5-5b  can  be  determined  by  dividing 
the  ring  area  by  the  effective  correlation  region. 

Because  of  the  stationarity  of  the  random  process  in 
the  data  space,  an  interesting  phenomenon  results  when  one, 
say,  doubles  the  initial  frequency  in  Fig.  5-5b  while 
keeping  other  parameters  unchanged.  The  area  in  the  ring 
will  be  approximately  doubled  and  thus  the  same  to  the  DOF. 
Although  the  range  (along  the  radial  direction)  resolution 
doesn't  change  by  such  a  new  arrangement,  the  azimuth  or 
along  angle  direction  does  double  its  DOF  by  the  doubling 
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Fig.  5-6. 


Correlation  of  the  data  and 
its  Fourier  transform 


uy, 

target  point 


Ar  -x2-x,  =  y,<£  for  small  <f> 


radar 


H 

radar 


Separability  of  azimuthal  and  range  processing 
at  small  angle  variation  4>.  (Target 
coordinates  fixed  on  paper) 


of  its  circumference  and  hence  number  of  independent 
samplings.  This  means  twice  (better)  the  azimuth  resolu¬ 
tion  although  the  radial  resolution  remains  the  same. 

This  phenomena  does  not  exist  in  one-dimensional  radar 
pulsing  where  the  range  data  is  the  only  information  sought. 
In  that  case  the  only  purpose  of  the  carrier  frequency  is 
to  carry  a  signal  of  certain  bandwidth  and  it  is  the  band¬ 
width  and  not  the  carrier  frequency  which  determines  the 
range  resolution.  However,  in  our  2-D  radar  system  the 
carrier  frequency  comes  into  the  radar  imaging  equation 
directly  and  manifests  itself  as  a  major  azimuthal  para¬ 
meter.  In  terms  of  radar  terminology,  a  raise  in  carrier 
frequency  increases  the  doppler  rate  and  thus  will  widen 
the  potential  azimuth  bandwidth  of  the  frequency  spectrum 
generated  by  a  single  target  point.  However  in  some 
situations,  e.g.  the  stripping  mode  of  SAR  in  which  the 
antenna  pattern  width  is  inversely  proportional  to  the 
wavelength  the  two  factors  cancel,  leaving  the  doppler 
spectral  width,  and  hence  the  system  performance  independent 
of  the  carrier  frequency.  [5-2] 

5 . 6  Actual  Reconstruction  Method 

Physically  the  radar  imaging  system  has  lots  of  differ¬ 
ences  from  the  tomography  projection  system  because  of  the 
many  widely  different  imaging  characteristics.  Some  of  the 
differences  have  been  described  in  the  end  of  Section  5.4. 
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In  the  radar  imaging  system  range  resolution  is 
obtained  by  the  ''timing"  of  the  signal  return  from  the 
target  point.  In  RAT  SCAT  multiple  frequency  signature  (MFS) 
is  used  and  the  timing  information  is  resolved  by  taking 
Fourier  transformation  along  the  stepped  frequencies. 
Ideally,  the  relative  motion  between  the  target  points  and 
the  radar  should  be  zero  to  obtain  range  resolution  with 
any  high  degree  of  precision  desired.  On  the  other  hand, 
the  azimuthal  resolution  is  introduced  by  the  created 
different  doppler  histories  at  different  azimuthal  points  by 
way  of  a  relative  motion  between  the  target  points  and  the 
radar.  The  seemingly  conflicting  requirements  are  resolved 
in  the  RAT  SCAT  system  in  which  different  frequency  compo¬ 
nents  at  an  aspect  angle  6  are  obtained  during  which  there 
is  no  relative  target  motion.  Azimuthal  information  is  then 
provided  by  the  phase  differences  of  the  same  frequency 
components  at  different  9's  due  to  the  range  change  of 
target  points  induced  by  the  target  rotation. 

In  Fig.  5-7  let  ( x1 , y 1 )  and  (x2,y2)  be  the  coordinates 
of  the  same  point  with  the  aspect  angle  difference  being  <J>. 
Similar  to  Eq.  (5-13)  we  have 

x2  =  x1cos4>+y1sin<I’  x-^  =  x2cos<t>-y2sin<i> 

or  (5-21) 

y2  =  y^osS-x^inG  y^^  =  y2coslJ)+x2sin4> 
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Thus  the  range  of  the  point  undergoes  a  change  of 


Ar  (x1  fY-^2  'y2  )  =  x2-xl 

=  x^cos<j>+y^  sin<j>-x^ 

=  x^ (cos$-l) +y^sin$  (5-22) 

which  creates  different  phse  histories  to  each  target 
point  in  going  from  one  radar  coordinate  (x^,y^)  to  the 
other  (J<2»y2^-  In  EcJ-  (5-22)  if  i>  is  small  enough,  then 

cos$  ~  1  and  Ar (x^ , y^ , x2 , y2 )  ~  y^sin$  x  y^$  (5-23) 

which  is  independent  of  x^  or  the  range.  The  phase 
variation  as  a  function  of  $  induced  by  r($>)  creates  a  dop- 
pler  frequency  ~  =  y.  Because  different  y  are  associated 
with  different  doppler  shifts,  the  azimuth  targets  can  be 
resolved  by  taking  Fourier  transforms  on  the  data  along 
small  azimuthal  angles,  as  long  as  the  data  are  taken  from 
appropriate  range  bins.  Usually  separability  of  azimuthal 
and  range  processings  are  desired  because  of  its  high 
efficiency  in  computation  [5-13],  in  which  the  two  process¬ 
ings  can  even  be  interchanged  if  so  desired.  To  insure 
the  separability  range  curvature*  has  to  be  avoided  by 
limiting  Ar  in  [5-22]  to  be  less  than  one  range  bin  width 
or  max  Ar<l  range  bin  width.  This  amount  to  a  coherence 
aspect  angle  txone  range  bin  width  in  the  target  with 
normalized  size.  Since  both  the  range  and  azimuth  require 

*The  term  "range  curvature"  in  the  subsequent  chapters  bear 
the  same  meanina  as  "motion  through  resolution  cells”  in 

[4-13]. 


1 10 


1 


Fourier  transformation  as  explained,  separability  implies 
using  2-D  Fast  Fourier  Transformations  (FFTs ' )  which  trans¬ 
forms  the  data  in  frequency  domain  to  the  reflectivity 
function  in  spatial  domain. 

Fig.  5-8  shows  a  section  of  data  for  such  a  small  fan 

angle  <t> .  B  is  the  range  bandwidth  with  range  frequency 

step  Af  and  B  and  B  are  the  minimum  and  maximum  azimuth 
Z1  z2 

bandwidths  at  f^  and  f^  respectively,  with  azimuth  frequen¬ 
cies  step  Af  and  Af  .  Note  that 
Z1  z2 


B  =  N  x  Af  =N*  f.A0 
z^  a  a  1 


B  =  N  x  Af  =  N  x  f  A0 
z2  3  Z2  3  Z 


(5-24) 


where  A9  is  the  azimuth  step  angle  in  radian  and  N  is  the 

a 

number  of  azimuth  points  coherently  processed. 

Discrete  Fourier  transforms  on  M  number  of  data  along 
range  frequency  with  frequency  step  Afr  and  bandwidth  B^ 
give  the  range  profile  with  unambiguous  spatial  dimension 
R  and  range  bin  width  Ar  where 


and 


Eqs. 

theorem. 


(5-24)  and 
Similarly 


r 

(5-25)  are  dual  relations 
in  the  azimuth  direction 


of 


(5-25) 


(5-26) 

sampling 
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where  Z^,Z2  are  the  unambiguous  spatial  dimensions  of  the 
Fourier  transformed  data  along  azimuth  direction  at  the 
lowest  and  highest  range  frequency  (i,e.,  f^  and  f2)  , 
respectively  and  Az^  and  Az2  are  the  azimuth  bin  widths  at 
frequencies  f^,f2  respectively.  Note  that  Az^  ss  Az2  or 
Z  «  Z2,  if  the  2-D  FFT  which  is  separable  is  valid  and 
creates  little  image  distortion. 

In  summary,  the  reasons  for  taking  the  2-D  FFT  on  small 
angle  data  are: 

1.  The  reflectivity  function  is  a  function  of  aspect 
angle. 

2.  Satisfactory  phase  compensation  for  the  propagation 
between  the  radar  and  the  target  center  is  extremely  diffi¬ 
cult,  if  not  impossible. 

3.  However,  the  reflectivity  function  can  be  assumed 
constant  over  small  aspect  angle  during  which  the  azimuth 
and  range  processing  can  be  separated  and  FFT  techniques 
can  be  employed. 

4.  The  shadowing  effect  can  be  reduced  to  a  minimum  by 
adopting  the  technique  of  coherent  processing  over  small 
angles  and  then  incoherent  summing  over  large  angles,  as 
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described  in  the  next  section. 


5 . 7  Incoherent  Summation  of  Different  Looks 

Coherent  processing  should  be  applied  to  only  those 
radar  returns  from  target  points  of  fixed  reflectivity 
during  the  imaging.  As  the  azimuth  aspect  changes  by  more 
than  a  certain  amount,  speckle  and  blocking  effect  change 
the  target  reflectivity  by  so  much  that  the  data  can  no 
longer  be  deemed  as  coherent.  Coherent  processing  on  the 
data  of  a  small  aspect  angle  span  gives  the  "look"  of  the 
target  from  that  particular  angle  and  shows  only  those 
target  surfaces  approximately  normal  to  that  particular 
LOS.  Those  looks  could  be  integrated  incoherently  to 
enhance  the  final  target  image.  The  procedure  will  be 
further  explained  in  the  following  chapter. 

5 . 8  Conclusions 

This  chapter  presents  the  theory  of  a  high  resolution 
radar  imaging  from  both  a  radar  systems  viewpoint  and  a 
degrees  of  freedom  or  numerical  analysis  viewpoint. 
Similarity  with  the  computer  aided  tomographic  scanner 
imaging  technology  is  pointed  out.  However  the  differences 
between  the  two  systems  are  emphasized  and  a  radar  unique 
reconstruction  algorithm  is  developed  for  combined  coherent 
and  noncoherent  imaging.  The  actual  reconstruction  method 
is  presented  which  will  be  used  in  Chapter  6  on  the  data 
obtained  from  the  RAT-SCAT  facility.  The  pictorial  images 
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resulting  from  the  computational  procedures  as  shown  in 
Chapter  6  are  surprisingly  recognizable  and  suggest  that 
these  techniques  may  have  some  practical  application  in  the 
future - 
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Chapter  6 

EXPERIMENTAL  RESULTS  OF  RADAR  TURNTABLE  IMAGING 

6 . 1  Introduction 

This  chapter  presents  the  experimental  results  of  the 
radar  turntable  imaging  in  Chapter  5.  A  complete  sequence 
of  processing  techniques  is  introduced  starting  with  the 
photographic  presentation  of  the  raw  data.  Explanations  as 
well  as  conjectures  of  the  intermediate  results  are  given 
along  the  way  towards  the  final  images.  Many  processing 
parameters  are  tested  against  the  qualities  of  images 
resulted . 

Three  aircraft  targets  are  individually  imaged  using 
actual  data  which  are  the  complex  returns  (in-phases  and 
quadrature-phases)  of  the  radar  signals  from  the  targets. 
Two  of  them  are  model  airplanes  (F-102,  F-5E)  of  reduced 
sizes  and  the  other  is  a  full-sized  real  F-102  plane. 
Experiments  are  done  for  two  setups  of  the  targets,  L.e., 
horizontal  and  vertical  positions , resulting  in  azimuth  and 
vertical  images,  respectively.  All  notations  bear  the  same 
meanings  as  in  Chapter  5. 

6.2  Imaging  a  Model  F-102  Airplane 

6.2.1.  Physical  properties  of  the  imaging  geometry 
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A  real  F-102  has  an  overall  length  and  a  wing-to-wing 
width  approximately  68  ft.  and  38  ft.  respectively.  The 
model  F-102  is  0.29  times  the  real  size.  Other  numerical 
values  used  in  the  experiments  are  listed  in  Table  6-1.  The 
imaging  geometry  is  depicted  in  Fig.  5-1.  The  bandwidth  in 
the  range  direction  is  approximately  870  MHz  which  has  an 
ideal  range  resolution  Pr 


c  ~  9. 84xl0^ft/sec  ~  ~  c_, 

P  r  =  T7r-  =  - jrL -  =  °*57 

r  2x870x10  /sec 

The  frequency  step  Af  =  3.4  MHz  corresponds 
guous  range  interval  lr 


ft.  (6-1) 

to  an  unambi- 


*r  =  =  132  ft-  (6'2) 

which  is  approximately  6.7  times  the  largest  dimension  of 
the  target.  This  means  that  the  data  has  been  oversampled 
by  a  factor  of  6.7  in  the  MFS  frequency  domain.  The  same 
conclusion  could  also  be  reached  from  the  DOF  analysis  in 
the  previous  chapter  by  applying  Eq.  (5-12).  The  step 
frequencies  in  the  azimuth  direction  are  between 


Af  ~ 

f  ■  xA0  =  31.87 

MHz 

(6-3) 

Z1 

mm 

Af  ~ 

Zn 

f  xA0  =  36.91 

max 

MHz 

(6-4) 

with  corresponding  unambiguous  intervals  12.9  ft.  and 
14.1  ft.  Refer  to  Fig.  5-8.  The  unambiguous  azimuth  inter¬ 
vals  are  longer  than  the  target  dimension,  which  means  that 
A9  =  0.2°  is  adequate  to  sample  the  data  in  the  azimuth 


117 


Table  6- 


f  . 
min 


max 

Af 

r 

M 


A  9 


1652  ft. 

1112  ft. 

540  ft. 

9.13  GHz 
10  GHz 
3 . 4  MH  z 
256 

0.2°  ss  0.0035  rad. 


Model  F  102  imaging  parameters 
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direction.  However,  the  sampling  is  barely  dense  enough 
because  the  unambiguous  interval  is  only  slightly  larger 
than  the  target  size. 

6.2.2  Source  data  from  model  F-102  airplane 

The  source  data,  which  are  the  returns  of  each  single 
frequency  waveform  at  each  particular  azimuth  angle  of  the 
target,  are  complex  numbers  consisting  of  both  the  in-phase 
and  quadrature-phase  components,  which  can  be  transformed 
into  polar  coordinates  to  become  the  magnitude  and  phase 
components.  Fig.  6-1  shows  the  phase  and  square  root  of  the 
magnitude  of  the  source  data  from  a  horizontally  mounted 
model  F-102  target.  They  are  256  data  points  along  lines  of 
constant  azimuth  angle  and  approximately  900  points  along 
the  azimuth  direction  from  0°  to  180°  of  the  target  orien¬ 
tation.  Note  that  only  180°  of  data  is  sufficient  for  the 
whole  target  reconstruction  since  the  target  was  presumed 
symmetric  horizontally.  The  square  rooting  is  necessary  to 
reduce  the  gigantic  dynamic  range  of  the  radar  cross  section 
( RCS )  of  an  airplane  as  a  function  of  the  aspect  angle  or 
radar  frequency  [6-1],  at  wavelengths  (ssl.5  cm  in  our  case) 
much  smaller  than  the  average  curvature  of  the  airplane. 

Note  that  at  approximately  60°  of  the  data  there  were  strong 
returns  probably  due  to  the  strong  reflectivities  of  the 
delta  wing  structure.  Similar  phenomenon  happened  close  to 
90°  or  broadside  position  where  the  fuselage  served  as 
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a)  Square  root  of  maqnitudc  b)  Phase  (0°-90°) 

(0-90°) 


c)  Square  root  of  maanitude  d)  r’hase  (90n-130  ’ 

(90  -1 8  0  ) 


Fiq.  6-1.  Source  data  of  the  horizontal  model  F-l  >2; 

horizontal  direction:  ranqe  step  frequency 
vertical  direction:  azimuth  anqle 
200  x  469  points  each  imauo 


120 


strong  reflectors.  Starting  from  approximately  120 u  the 
gain  of  the  radar  system  seemed  raised  resulting  in  higher 
brightness  of  the  magnitude  part  of  the  photographs.  The 
phase  part  of  the  source  data  are  linearly  coded  from  -n  to 
it  for  the  photograph.  The  data  is  in  the  frequency  domain 
in  both  the  range  and  the  azimuth  dimensions.  Because  of 
the  linear  phase  due  to  y,-Yq-Y0  as  explained  in  Chapter  5 
the  oversampling  phenomenon  along  the  range  is  not  very 
obvious.  It  turns  out,  as  shown  in  Fig.  6-2,  that  this 
linear  phase  contributes  a  fixed  amount  of  shift  in  the 
spatial  domain.  The  highly  structured  phases  are  probably 
a  result  of  a  large  signal  noise  ratio  of  the  radar  system. 

6.2.3  Fourier  transforming  in  the  range  direction 

Because  the  source  data  was  in  two-dimensional  frequ¬ 
ency  domain,  some  form  of  two-dimensional  transformation 
was  necessary  to  bring  the  data  into  the  spatial  domain.  As 
explained  in  Chapter  5,  we  segment  the  data  according  to 
their  azimuth  angle  into  sectors  and  then  apply  2-D  Fourier 
trans formations  on  each  sector  which  is  approximately 
rectangular  or  separable.  Fourier  transforming  the  data 
along  a  fixed  aspect  angle  gives  the  "projections"  of  the 
target  onto  that  line  of  sight  (LOS)  at  that  angle.  If  the 
processing  is  assumed  separable  between  the  range  and  the 
azimuth  then  the  magnitude  parts  of  the  projections  provide 
the  range  information  (range  profiles)  while  the  phase 


parts,  when  processed  with  the  phases  at  neighboring  aspect 
angles,  will  resolve  the  azimuth  information.  Toward  this 
assumption  our  source  data  was  Fourier  transformed  along 
the  range  direction  for  each  azimuth  angle  and  the  phase 
and  the  magnitude  parts  of  the  resultant  data  are  then  shown 
in  Fig.  6-2.  Note  that  now  the  data  is  in  spatial  domain 
along  the  range  and  still  in  frequency  domain  along  the 
azimuth.  Again  for  a  better  presentation  square  rooting  is 
applied  to  the  magnitude  parts,  which  assume  a  dynamic  range 
of  more  than  30  dB.  The  magnitude  parts  clearly  show  the 
location  of  the  target  which  is  off  the  center  of  the  pic¬ 
ture  because  of  the  linear  phase  and  the  target  occupies 
giy  of  the  whole  range  dimension  as  a  result  of  the  over- 
sampling  in  the  range  frequency  domain  of  the  source  data. 
Note  the  strong  returns  curving  in  and  out  in  the  figure 
along  the  azimuth  angle.  The  bright  lines  in  the  center  of 
the  magnitude  pictures  account  for  the  DC  bias  of  the 
radar  receiver  and  processing  unit.  The  faint  mirror 
images  of  the  target  with  respect  to  the  central  DC  line  are 
due  to  the  imbalance  of  the  in-phase  and  quadratic-phase 
circuitry  of  the  radar  system.  To  show  this,  let  D^  and  D^ 
be  the  gains  of  in-phase  and  quadrature  phase  circuitries 
of  the  mixer  of  the  radar  receiver.  From  Eq.  (5-10),  the 
actual  complex  data  becomes 
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I (i,k)+j8 (i,k) 


exp^-j2iTf 

where 

A  f3 

G(k,i)  =  g  (x,  i)  exp  ( j  2irkx)  dx  (6-6) 

-a 

The  imbalance  of  the  second  term  in  Eq.  (6-5)  causes  the 
mirror  images  in  Fig.  6-2.  It  is  reasoned  that  in 
Fig.  6-2  (a)  the  strong  returns  closest  to  the  radar  (which 
is  to  the  right  of  the  photographs)  are  due  to  the  nose  and 
those  in  Fig.  6-2  (c)  are  due  to  the  tail  of  the  target.  In 


Fig.  6-2  (a)  observe  a  bright  tiny  spot  very  close  to  the 
nose  line  at  approximately  60°  azimuth.  Again,  it  is  due  to 


the  returns  from  the  delta  wing. 

The  phase  figures  6-2 (b)  and  (d)  are  correlated  with 
their  counterparts  (a)  and  (c)  in  that  prominant  structures 
exist  only  in  those  regions  where  there  are  strong  magni¬ 
tudes  due  to  target  points  of  high  reflectivity.  In  other 
regions  the  phases  look  random  and  noise  like. 
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There  are  two  ways  to  remove  the  redundancy  of  the 
data  in  the  range  dimension.  One  is  to  choose  a  larger 
frequency  step  ,\f  from  the  DOF  analysis  in  Chanter  5  to  keep 
only  minimum  frequency  components  required  ir  the  source 
data.  The  other  is  to  retain  in  the  spatial  domain  the 
portion  containing  the  target.  rn,e  latter  option  has  the 
advantage  of  offering  a  higher  S/N  in  the  reduced  data.  In 
our  experiments  we  always  keep  64  range  bins  in  the  spatial 
domain  which  gives  a  data  reduction  of  4:1.  The  windowing 
of  the  data  is  visually  determined  from  Fig.  6-2 ia)  and  (c). 
Fine  tuning  of  the  target  center  is  necessary  because  of  the 
unavailability  or  uncertainty  of  and  y^.  This  will  be 

discussed  in  subsection  6.2.5. 

6.2.4  Azimuth  processing 

If  the  reconstruction,  which  is  basically  a  two- 
dimensional  transformation  as  explained  before,  can  be 
assumed  separable,  then  since  the  range  transformation 
has  been  accomplished  in  the  previous  subsection,  only 
azimuth  transformation  remains  to  be  done.  In  our  case  of 
RAT  SCAT  imaging,  separability  means  no  range  curvature  or 
range  bin  moving  problems  induced  by  the  rotation  of  the 
target  necessary  for  azimuth  resolution.  Since  the  maximal 
radial  extent  of  the  target  is  assumed  to  be  no  more  than 
10  feet,  the  range  change  will  be  less  than  10*<J)feet  for  all 
the  target  points  when  the  target  undergoes  a  rotation  of 
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<p  radians  assumed  so  small  that  t  anj  «  p .  If  the  range 
walking  is  to  be  prohibited  at  all,  then 

10x0  <  pr  =  0.52  ft.  (6-7) 

or 

4>  £  3°  (6-8) 

However,  it  is  sometimes  desirable  to  process  coherent 
intervals  larger  than  the  maximum  interval  without  range 
walking  even  if  separable  processing  is  assumed.  This  has 
the  advantage  of  a  better  azimuth  resolution,  although  it 
will  be  subjected  to  some  blurring  due  to  the  range  walking 
introduced.  The  quantitative  analysis  of  the  interaction 
between  the  range  walking  and  azimuth  bandwidth  will  not  be 
pursued  here.  It  is  pointed  out,  however,  that  if  one  were 
willing  to  use  a  non-separate  model  to  process  the  data, 
then  as  long  as  the  target  reflectivity  function  could  be 
reasonably  assumed  to  be  constant  within  the  coherent  azi¬ 
muth  interval,  the  interval  processed  could  be  as  large  as 
possible  without  actually  suffering  from  range  walking 
problem.  This  approach  would  yield  optimum  DOF  of  the  final 
images  and  hence  better  resolutions.  It  would,  of  course, 
demand  much  more  computational  requirement  as  described  in 
Chapter  2  and  was  not  adopted  in  our  experiments. 

In  our  experiments  the  data  in  the  ring  area  in 
Fig.  5-3  was  divided  into  possibly  overlapped  azimuthal 
sectors.  Each  sector  is  individually  Fourier  transformed 
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along  azimuth  (recall  we  have  performed  range  Fourier 
transformation  in  6.2.3)  and  only  the  magnitude  part  which 
is  of  our  interest  is  kept  to  give  target  "looks"  at  the 
respective  aspect  angles.  Figure  6-3  (a)  shows  the  image 
looks  centered  at  0°,  30°,  60°,  120°  and  150°  with  the 
azimuthal  coherent  interval  being  3.2°  or  16  points. 

Figure  6-3 (b)  and  6-3  (c)  are  looks  at  the  same  angle  with 
6.4°  and  12.8°  coherent  intervals,  respectively.  Note  the 
different  relative  scales  between  azimuth  and  range  in 
Fig.  6-3  (a),  (b)  and  (c)  in  which  figure  6-3  (a)  and  6-3 (b) 

have  been  two-d imensionally  linearly  interpolated  at  a 
factor  of  4  and  2  respectively  to  fully  utilize  the  dimen¬ 
sion  of  the  displaying  monitor  from  which  all  of  our  2-D 
images  were  shot.  In  all  the  image  looks  of  Fig.  6-3  the 
radar  is  to  the  right.  Salient  features  of  the  airplane, 
e.g.,  the  nose,  the  tail,  the  delta  wing  and  the  fuselage 
can  be  readily  identified  at  appropriate  looks.  For 
example,  at  0'J  the  target  is  sitting  with  nose  on  the  radar. 
The  images,  which  show  reflectivities  of  the  target,  have  a 
very  large  dynamic  range  of  brightness  and  the  sidelobe 
effects  are  very  serious  which  manifest  themselves  as 
cross-like  structures  running  horizontally. and  vertically, 
corresponding  to  the  two  dimensions  of  processing.  Theore¬ 
tically,  images  without  sidelobes  could  be  obtained  only 
from  all  of  the  data  points  in  the  whole  frequency  domain 
which  extends  from  -«°  to  °°  in  both  dimensions.  In  reality 
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a)  Azimuth  coherence 
intervals  =  3.2 
(16  points) 


b)  Azimuth  cohereng 
intervals  -  6.4 
(32  points) 


c)  Azimuth  coherencg 
intervals  =  12. R 
(64  points) 


Model  F-102  azimuth  images  at  0°,30°, 
90  , 120°  and  160°  asoect  angles 


F  i  q .  6  -  3  . 


our  data  for  each  look  comes  only  from  an  approximately 
rectangular  area  of  the  frequency  space.  This  means, 
equivalently,  a  rectangular  function  has  been  multiplied 
with  the  data  on  the  otherwise  infinitely-extended  data 
space.  In  the  spatial  domain  the  effect  is  as  if  the  ideal 
images  had  been  convolved  with  the  Fourier  transformation 
of  a  2-D  rectangular  function,  which  is  a  sine  function. 

The  sine  function  has  a  first  sidelobe  of  13.4  dB  below 
the  main  lobe  level.  Compared  to  the  30  dB  dynamic  range 
of  the  strength  of  the  target  reflectivities,  the  -13.4  dB 
sidelobe  could  completely  mask  the  weak  signals  close  to  the 
strong  ones.  To  alleviate  the  sidelobe  problem,  a  well¬ 
shaped  function  should  be  multiplied  to  the  data  before 
Fourier  transformation  is  taken  [6-2] .  Figure  6-4  shows  the 
image  looks  when  a  two-dimensional  Hanning  window,  which  has 
a  first  sidelobe  less  than  -30  dB,  is  incorporated  into  the 
data.  Comparing  Fig.  6-3 (b)  and  Fig.  6-4  it  is  seen  that 
the  sidelobes  have  been  reduced  and  the  main  lobe  width 
slightly  increased  as  predicted  in  [6-2]. 

To  verify  that  the  Hanning  windowed  Fourier  transforma¬ 
tion  is  a  valid  image  reconstruction  technique,  Fig.  6-6, 
a  sequence  of  28  abutting  looks  with  each  look  occupying 
6.4°  of  data  azimuthally,  is  presented.  The  rotation  of  the 
fuselage  is  very  clear  as  it  goes  from  0°  to  180  '  in  the 
consecutive  looks.  In  all  the  following  experiments  the 
Hanning  window  is  alv/ays  incorporated  in  the  FFT's. 
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6.2.5  Incoherent  summation  of  the  looks 

Because  of  the  specularity  of  the  target,  each  look  in 
Fig.  6-5  shows  only  those  target  points  which  contribute 
coherently  to  the  radar  returns  at  the  particular  aspect 
angle.  It  is  reasoned  that  a  suitable  combination  of 
different  looks  should  give  a  much  enhanced  overall  target 
image.  Indeed,  we  used  the  following  procedure  to  accom¬ 
plish  this: 

(1)  determine  the  center  of  rotation  for  each  look 
(and  hence  all  of  the  looks  since  the  distance  between  the 
radar  and  the  rotation  center  is  fixed  upon  the  data 
collecting) . 

(2)  rotate  each  look  according  to  their  respective 
aspect  angle. 

(3)  incoherently  sum  the  rotated  and  registered  looks 
i.e.  only  the  magnitude  part  of  each  look  is  used  in  the 
summation. 

An  additional  advantage  of  the  incoherent  summation  is  the 
reduction  of  the  "speckles"  due  to  the  phase  errors  in  cohe¬ 
rent  systems. 

In  the  case  of  "azimuth  imaging",  (i.e.  horizontal 
setup  of  target) ,  the  images  obtained  from  180°  orientation 
are  always  flipped  for  the  other  180°  and  added  for 
symmetry  after  the  incoherent  summation. 

There  are  three  parameters  independently  "tuned"  to 


reach  final  summed  images: 


(1)  the  coherence  interval  in  the  azimuth  direction. 

(2)  the  amount  of  sector  overlapping  for  adjacent 
looks . 

(3)  the  center  of  rotation  assumed  in  all  looks. 

Each  of  the  above  parameters  is  tested  in  our  experiments  to 
show  how  they  affect  the  image  quality. 

Figure  6-6  through  6-8  test  the  parameters  (1)  and 
(2).  Overlapping  factors  of  abutting,  50 %  overlap  and  90% 
overlap  and  coherence  intervals  of  3.2°,  6.4°  and  12.8°  are 
adopted  in  these  figures.  Figure  6-9  test  different  rotation 
centers  (different  range  bins)  when  the  other  parameters 
are  fixed.  The  image  quality  seems  to  be  relatively  insen¬ 
sitive  to  the  rotation  center  of  each  look  compared  to  the 
case  of  F-5E  imaging  to  be  discussed  later.  Generally 
speaking,  the  experimental  results  show  considerable 
improvement  over  the  single-look  images  and  clearly  show  the 
outline  of  the  characteristic  delta  wing  of  the  F-102 
aircraft. 

6.2.6  Imaging  the  vertical  model  F-102 

The  imaging  principle  of  the  RAT  SCAT  is  that  the 
information  of  the  target  reflectivity  on  the  plane  of  the 
rotation  can  be  resolved.  In  other  words  the  images  should 
outline  the  vertical  projections  of  the  target  on  the 
ground.  Thus  by  using  the  data  gathered  from  the  target 
airplanes  which  sat  vertically  (i.e.,  roll  angle  =  90°)  one 
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a) 


Abutt inq 


b)  50*  overlapping 


c)  90  overlapping 


Fig.  6-6.  Incoherently  inteqrated  azimuth  images  of 
model  F- 102  with  different  overlapping 
factors 

aspect  coherence  intervals  =  3.2 

centers  of  rotation  =  range  bin  number  183 


1  35 


a)  Abutting 


b)  50%  overlapping 


Fig.  6- 7  . 


Counterparts  of  Fig.  6-6  with 
f>.4°  aspect  coherence  intervals, 
(all  other  parameters  same  as  Fig 


a)  Abutting 


b)  50"  overlapping 


c)  90  overlapping 


Fig.  6-8.  Counterparts  of  Fiqs.  6-6  and  6-7 
12.8°  aspect  coherence  intervals, 
other  parameters  same  as  Fiqs.  6- 

6-7) 


with 
(all 
6  and 


a)  Centered  at  range  bin 
187 


b)  Centered  at  range  bin 
188 


Fig.  6-9.  Incoherently  intenrated  azimuth  images  of  model 
F-102  with  different  centers  of  rotation 
aspect  coherence  intervals  6.4 
overlapping  factor:  abutting 


1 18 


fc-  • 


should  come  up  with  the  "side  view"  or  elevation  image  of 
the  airplane.  The  processing  principles  are  much  the  same 
as  described  before.  However  in  this  setup  the  radar  images 
are  not  symmetric  anymore  so  the  data  from  the  full  360° 
rotation  is  necessary  for  our  reconstruction  algorithm  of 
incoherently  summing  different  looks. 

Figure  6-10  shows  the  result  of  the  model  F-102  from 
the  vertical  setup.  The  azimuth  coherence  intervals  were 
6.4°  ,  the  looks  were  abutting  and  the  center  of  rotation 
was  at  range  bin  number  178. 

6 . 3  Imaging  the  Model  F-5E 

The  same  processing  techniques  were  applied  to  an  F-5E 
model  airplane  which  has  a  different  shape  and  size  from 
those  of  the  F-102.  The  size  of  the  model  target  is  0.30 
times  that  of  an  actual  plane  which  is  approximately  48  feet 
in  length  and  26  feet  from  wing  tip  to  wing  tip.  and  Y2 

were  unknown  constants  while  other  parameters  were  the  same 
as  those  in  the  case  of  the  model  F-102. 

In  the  followinq  experiments  the  target  was  assumed  in 
a  horizontal  position  (0°  roll)  unless  otherwise  stated. 

Figure  6-11  shows  a  resultant  sequence  of  14  looks  with 
abutting  6.4°  coherence  intervals,  covering  0°  to  90° 
aspect  angles.  Note  the  rotation  of  the  fuselage  along 
the  consecutive  looks.  Aqain  the  radar  was  to  the  right  for 
each  look  so  the  right  wingtip  can  be  identified  in  the 
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a)  Abutting 


b)  50*  overlapping 


Elevation 
di f ferent 
coherence 
centers  of 


images  of  model  F-102 
overlapping  factors; 
intervals:  6.4 

rotation:  range  bin 


w  i  t  h 
aspect 

numbo 


last  few  looks  in  which  the  target  was  approximately 
broadsided  to  the  radar.  The  other  14  looks  covering  90° 
through  180°  aspect  angles,  which  are  not  shown  here,  have 
the  corresponding  orientation  of  the  previous  looks.  These 
28  looks  were  incoherently  combined  and  flipped  for  symmetry 
with  the  rotation  center  being  at  range  bin  number  177.25 
to  become  Fig.  6-12.  Note  the  protrudent  wings  of  this 
aircraft  as  opposed  to  the  delta  wings  of  an  F-102.  As 
mentioned  before,  the  quality  of  the  image  for  the  model 
F-5E  seems  more  sensitive  to  the  rotation  centered  assumed; 
at  approximately  4  bins  (center  at  #181)  off  the  true 
center  the  target  image  looks  like  a  big  blub  as  in 
Fig.  6-13.  Figure  6-14  is  a  fine  tuning  of  the  rotation 
center  at  0.25  range  bin  resolution. 

Figure  6-15  shows  the  effect  of  the  overlapping  factor 
between  adjacent  looks. 

Figure  6-16  are  the  elevation  images  of  the  target 
with  a  vertical  position  of  the  target.  The  images  are 
centered  at  range  bin  number  178  with  different  amounts  of 
overlapping . 

6.4  Comparison  Between  the  Model  F-102  and  F-5E  Images 

Figure  6-17  presents  a  summary  of  photograDhs  for  the 
model  F-102  and  model  F-5E  airframes  for  both  azimuth  and 
elevation  plots.  The  photographs  show  the  distinct  shapes 
of  the  wings  of  the  two  aircraft.  Except  for  a  slight 
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Fig.  6-12.  Model  F-5E  azimuth  image  aspect 
coherence  intervals  3.2° 
overlapping  factor:  abutting 
center  of  rotation:  range  bin  177.25 


Centered  at  ranqe  bin 
177 


b)  Centered  at  range  bin 
177.25 


c)  Centered  at  range  bin  177.5 


Fig.  6-14.  Fine  tuninq  of  the  rotation 

center  for  model  F-5E  azimuth 
images 


a)  50%  overlapping 


b)  Abutting 


Fig.  6-15.  Effect  of  overlanpina  factors  tor  model 
F-5E  azimuth  images 
aspect  coherence  intervals: 6 . 4 
centers  of  rotation:  range  bin  178 
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a)  50%  overlapping 


b)  Abutting 

Model  F-5E  elevation  imaqes  with 
different  overlapping  factors, 
aspect  coherence  intervals  6.4 
centers  of  rotation:  range  bin  178 


a)  F-102  azimuth  image  b)  F-102  elevation  image 

(28  looks)  (56  looks)  (1/2  scale) 


c)  F-5E  azimuth  image  d)  F-5E  elevation  image 

(28  looks)  (56  looks)  (1/2  scale) 


Fig.  6-17.  Model  F-102  and  modeloF-5F,  azimuth  and 
elevation  images  (6.4(  coherence) 
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difference  of  the  physical  scales  of  these  models  (0.29  vs 
0.30)  all  other  parameters  are  fixed  for  these  images. 
Consequently,  it  is  clear  that  the  F-5E  is  a  smaller  air¬ 
craft  and  naturally  has  a  different  azimuth  and  elevation 
projection  than  does  the  F-102. 


6 . 5  Imaging  a  Real  F-102  Airplane 

The  third  target  imaged  was  a  full-sized  real  F-102 
airplane . 

5,552  ft.,  2,744  MHz  and  1,000  kHz,  respectively  to  account 


In  this  case,  Yq,  an^  Afr  were  changed  to 


for  the  large  size  of  the  real  plane.  and  y£  were  un¬ 

known  constants  and  other  parameters  remained  the  same  as  in 
the  previous  cases  From  the  parameters  it  is  found  that 
the  unambiguous  range  interval  Ir  is 


Ir  =  2x A  f 


=  449  ft. 


and  the  range  bandwidth  is 


B  =  256  x  1,000  kHz  =  256  MHz 
r 


which  ha..  an  ideal  range  resolution  of 


P>-  = 


r  2B 


1.76  ft. 


The  data  ring  in  Fig.  5-3  has  a  smaller  radius  and 


width,  hence  smaller  area  than  in  the  cases  of  model  air¬ 
craft.  However,  because  of  the  larger  size  of  the  target, 
the  correlation  field  in  the  data  space  has  a  smaller 
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"effective"  autocorrelation  region  and  thus  the  DOF  in  the 
small  data  ring  is  approximately  the  same  as  in  the  model 
airplane  cases  because  of  the  cancellation  of  the  scaling 
factors.  Thus  it  is  expected  that  the  quality  of  the  target 
images  should  be  about  the  same.  Roughly  speaking,  a  target 
of  larger  spatial  extent  needs  a  narrower  sampling  interval 
in  the  frequency  data  domain  to  avoid  any  aliasing  effect. 
Consequently  the  area  in  the  sampled  data  domain  would  be 
smaller  if  the  same  number  (DOF)  of  independent  data  points 
which  determines  the  quality  of  the  final  images  is  to  be 
present . 

Figure  6-18  presents  some  photographs  of  the  real 
F-102  radar  images  at  different  rotation  centers.  As 
predicted,  the  images  qualities  are  about  the  same  as  those 
in  Fig.  6-9.  Again,  here  the  delta  wing  structure  is 
clearly  visible. 

6 . 6  Summary  and  Conclusions 

Reconstruction  of  the  reflectivity  function  of  three 
airplanes  is  attempted  experimentally  in  this  chapter. 

Three  major  parameters  of  the  reconstruction  algorithms 
proposed  in  Chapter  5  are  tested.  They  are:  aspect 
coherence  interval  (the  optimal  value  of  which  is  deter¬ 
mined  by  the  interaction  between  the  azimuth  resolution  and 
the  range  walking  when  separable  processing  is  assumed) , 
factor  of  overlapping  between  adjacent  looks,  and  location 
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of  the  rotation  center  which  has  a  "focusing"  effect.  As  it 
turns  out,  some  targets  are  sensitive  to  some  of  the  para¬ 
meters  while  insensitive  to  others.  It  is  very  difficult 
to  predict  what  parameters  should  be  used  in  a  given 
imaging  condition.  For  our  purpose  of  obtaining  the  "best" 
images  judged  by  human  perception,  trial  and  error  seems 
to  be  the  only  feasible  way  at  the  present  time. 

In  the  final  stage  of  our  experiments  to  take  pictures 
off  the  TV  monitor  for  display,  most  factors  of  the  shooting 
conditions  were  well  controlled  and  calibrated.  In  fact, 
most  of  our  images  were  taken  from  the  same  camera  at 
fixed  focal  length,  shutter  speed  and  aperture.  The 
brightness  and  the  contrast  on  the  monitor  were  fixed  within 
the  accuracy  of  the  monitor  operation.  The  background 
intensity  has  always  been  kept  as  low  as  possible.  The 
major  nonuniformity  of  the  pictures  so  presented,  if  any, 
is  conjectured  due  mainly  to  the  films  because  of  the 
film  grain  noise  and  possible  difference  in  the  sensitivity 
to  the  exposure. 
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Chapter  7 

RADAR  IMAGING  WITH  TARGET  MOTION 

7 . 1  Introduction 

Imaging  from  ground-based  (stationary)  radars  of 
moving  targets  is  often  possible  by  utilizing  a  "synthetic 
aperture"  developed  from  the  target  motion  itself,  or 
sometimes  called  "inverse  SAR"  for  the  obvious  reason. 

This  chapter  addresses  the  theory  about  a  ground-based 
radar  imaging  system  in  which  a  target  aircraft  is  imaged 
by  its  own  motion  induced  doppler.  The  aircraft  is  imaged 
from  both  a  straight  flight  and  a  turn  with  recognizable 
results.  Analysis  shows  that  two  phase  components  exist 
in  the  radar  return,  one  being  gross  velocity  induced,  the 
other  being  interscatterer  interference  within  the  target 
itself.  The  former  phase  must  be  removed  prior  to  imaging 
and  techniques  are  developed  for  this  task.  Preprocessing, 
range  curvature,  range  alignment,  motion  compensation,  and 
presuming  are  all  addressed  in  this  chapter. 

7.2  Radar  Imaging  Equations 

Figure  7-1  shows  the  flight  path  of  a  target  aircraft 
which  has  an  overall  length  of  approximately  80  feet  and 
wing  span  of  about  70  feet.  Two  portions  of  the  flight  path 
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along  which  the  data  were  obtained  for  imaging  will  be 
called  interval  1  and  interval  2,  as  shown  in  Fig.  7-1. 

The  first  interval  is  when  the  airplane  was  flying  straight, 
at  angles  approximately  30°  to  15°  off  broadwide,  whereby 
the  second  interval  occurs  when  the  airplane  was  making  a 
standard  left  turn. 

7.2.1  First  interval 

Figure  7-2  is  a  reproduction  o'  the  first  interval  of 
Fig.  7-1.  The  aspect  angle  of  the  target  center  viewed 
from  the  radar  undergoes  a  change  ,  which  in  this  case  is 
the  same  as  the  change  of  the  aspect  angle  of  the  target 
body  with  respect  to  the  radar  line  of  sight  (LOS)  because 
of  the  straight  flight.  In  fact,  it  can  be  shown  that  it 
is  the  latter,  and  not  the  former  angle  change,  which 
provides  the  azimuthal  information.  For  this  reason,  we 
redraw  Fig.  7-2  using  the  target  center  as  the  origin  of 
the  coordinate  system.  This  becomes  Fig.  7-3.  Observe  a 
close  resemblance  to  the  rotational  geometry  of  a  turntable 
imaging  system  as  in  chapters  5  and  6  (7-1]. 

In  Fig.  7-3  let  coordinates  (£,n)  be  fixed  on  the 
target  and  (x,y)  be  rotating  with  the  radar  and  with  the 
same  origin  0  as  (£,n).  Let  the  angle  from  axes  (  to  x  be 
0 .  Then  we  have 

x  =  £cos0+nsin8  £  =  xcos0-ysin0 

or  (7-1) 

y  =  -£sin0+ncos8  n  =  xsin0+ycos0 
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Fig.  7-3.  Radar  trajectory  relative  to  target 
center  in  first  interval 
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Let  (xr,yr)  or  (£r*nr)  be  the  coordinates  of  the  radar  and 
(x  .y  )  or  (E  *n  )  be  those  of  any  target  point.  Note  that 

3  3  3  3 

Yr  =  0  identically.  Let  £r  be  the  £  coordinate  of  the 

o 

radar  had  the  airplane  flown  to  the  position  broadwide  to 
the  radar  and  define  R  =  | £  | .  Let  s (x  ,y  ,x  ,y  )  be  the 

O  IT  q  3  3  r  XT 

distance  between  (x  ,y  )  and  (x  ,y  )  .  Then 

33  it  r 


(xa  ,ya,xr 'yr^  =  ^xa“xr)  2+ *ya"yr' 2] 


=  f  x^ -  2 x  x  + 
\  r  r  a 

-  (“r^Va) 

X  X 

~  r  a 

"  xr  x 

r 


2  2' 
x  +x  +y 
a  a  a  j 


(7-2) 


for  target  distance  large  enough  such  that  x  +y  <<  x  . 

33  r 

Note  that  x„  and  x=  are  both  functions  of  time  t. 
r  3 

As  the  flight  goes  by,  the  change  in  s (x  ,y  ,x  ,y  ) 

4  a  3  r  r 

induces  a  Doppler  frequency  of  the  target  point  at  (x  ,y  ) 

3  3 


fd(xa'ya: 


2  ds  2  (  i 
X  dt  X  V  dt 


dx  dx  \ 

- i 

dt  dt 


(7-3) 


where  X  is  the  wavelength  of  the  transmitted  signals. 

Equation  (7-3)  has  a  very  interesting  interpretation. 


The  Doppler  frequency  of  the  target  point  (x  ,y  )  is 

3  3 

comprised  of  two  part: 


.  _  dx  / 
4  2  _£  (, 
\1  X  dt  y 


when  -g—  =  Ol,  due  to  the  trajectory  of 
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the  target  center,  and 


.  t  dx 
A  2  a 

X  ~3F 


(when  =  Oj , 


due  to  the  relative  motion 


of  the  target  inter-scatterer  points  with  respect  to  the 

target  center.  Recall  that  in  the  last  two  chapters  the 
dx 

term  -g—  was  zero  because  the  target  was  rotating  on  a 

dx 

table  and  hence  x  was  a  constant.  Also  in  RAT  SCAT  -rr 

r  dt 

is  a  constant  because  the  rotation  rate,  or  the  step 
increment  in  aspect  angle  of  the  target  is  a  constant.  If 
we  are  to  apply  the  techniques  developed  in  turntable  ima¬ 
ging  to  the  underlying  system,  these  two  differences  have 
to  be  compensated  first. 


Because 


2  4 

2  2U  _  1  V  _  1  ]r 

°  7  0  2  Ro  8  R3 


n  dn  ,  n"5  dn 

r  r  1  r  r 

R~  ~dF  2  13  “dF 

0  Ro 


(7-4) 


(7-5) 


Assuming  a  constant  flight  speed  v,  from  Fig.  7-3  we 

dnr 

have  j-.  =  -v  and  hence 

at 

dxr  r|rv  .  1  v  3  „  ,, 

aF  T-  *  I  j  \  (7‘6> 

o  R 

o 

Thus  the  Doppler  frequency  induced  by  xr  is  a  linear 
function  of  nr  plus  higher  order  terms  and  the  phases  of 
the  return  due  to  xr  are  quadratic  plus  other  higher  order 
terms,  the  first  two  being  of  opposite  signs.  It  is  clear 


from  Eq.  (7-6)  that  as  the  aircraft  flies  closer  to  being 


broadside  or  normal  to  the  radar  LOS  corresponding  to  a 
decreasing  nr»  the  linear  relation  between  and  is 

more  precise.  In  that  situation  the  "rate"  in  Doppler 
frequency 


,2  .  _  -  2  2  . 

d  *r  .  2  /  v2  3nrv  \ 

at2  X  \  Ro  2R3  / 

o' 


(7-7) 


becomes  larger  and  larger  and  converges  to  a  constant  pro- 
2 

portional  to  ,  which  in  terms  of  signal  processing  is  a 
Ro 

linear  FM  rate. 

Rough' y  speaking,  the  history  of  xr  does  not  provide 

useful  imaging  information  but  only  serves  as  a  carrier  to 
dx_ 

carry 


dt 


which  contributes  to  coherent  imaging  ability  by 


providing  a  wide  spectral  bandwidth  inherent  in  a  linear  FM 

signal.  Upon  reception,  xr  phases  should  be  removed  first 

(sometimes  called  de-chirping  or  motion  compensation)  to 

ease  the  processing  afterwards  to  extract  the  relevant 

dx„ 

Schemes  to  remove  these 


azimuthal  phase  information  in 


dt  ' 


phases  will  be  described  in  Section  7-5.  We  now  analyze  the 

dx 

higher  order  effects  of  which  does  not  exist  in  RAT  SCAT. 


Since  x^  =  £  cos0+n  sin8 , 

a  a  a 


Also 


t  dx 
2  _ a 

X  dt 


2 .  r  •  q .  Q , d9 

I(-£asine+nacos0)3r 


2y 


a  d_9 
dt 


tanQ 


or 


(7-8) 
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so 


ae 

dt 


1  + 


v 

R 


( — v ) 


1  + 


(7-9) 


d0 


Again,  ^  approaches  maximum  as  the  airplane  flies 

2y. 


normal  (nr=0)  to  the  radar 


_  .  d0  v  , 

In  that  case  -3-r-  =  ana 
at  R 

o 


a  v 
R 


:7-io) 


2  o 

This  means  that  at  positions  close  to  broadside  the  doppler 
frequency  of  a  single  target  point  at  (x.y  )  due  to  change 
in  x_  is  proportional  to  y_  and  is  independent  of  x  .  In 

a  a  a. 

this  situation  there  is  a  one-dimensional  (along  y  )  Fourier 

a 

transform  relation  between  the  returns  and  the  target 

reflectivities.  Over  a  wide  range  of  the  flight  path,  the 

factor  — i-j  could  change  so  much  as  to  affect  the  perfor- 
nr 

o 

mance  of  the  Fourier  transform.  However,  over  a  small 
interval  where  — is  approximately  constant  and  can  thus 


1+ 


3 


be  replaced  by  a  nominal  value,  Fourier  transform  techniques 
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can  still  be  applied.  Thus  if  one  does  not  c urtuensnte  for 
the  above  abnormality  factor,  there  is  a  tradeoff  between 
the  resolution  (determined  by  the  aspect  anqle  span  cohe¬ 
rently  processed)  and  the  fidelity  of  a  Fourier  transform. 
One  way  to  compensate  for  this  is  to  interpolate  the  azi¬ 
muthal  samples  to  obtain  the  data  uniformly  soaced  in  the 

azimuthal  frequency  domain.  Another  possibility  is  to  use 

dO 

a  variable  PRF  which  follows  or  the  rate  of  rotation  of 

dt 

the  target  as  seen  from  the  radar.  Note  from  Eq.  (7-10) 
that  an  increase  in  v,  or  a  decrease  in  .1  or  Rq  'will 
increase  the  obtainable  resolution  accompanied  with  a 
higher  minimum  PRF  required.  The  A -dependence  will  be 
further  explained  later.  The  "speckles”  of  the  final  images 
from  the  coherent  processing  could  be  removed  somewhat  by 
incoherently  averaging  different  coherent  looks,  as  in 
Chapters  5  and  6. 

7.2.2  Second  interval 

The  flight  path  of  the  second  irferv-'.  .s  depicted  in 

Fig.  7-4  in  which  (C,',n')  are  coordinates  centered  and  fixed 

on  the  target  and  (x,y)  are  fixed  on  the  radar  with  the 

origin  o  being  the  center  of  the  turn.  Let  ( '  ,  •  i )  be  the 

coordinates  parallel  to  (t',n'),  with  origin  at  o.  Again 

let  the  coordinates  of  the  radar  and  the  target  points  have 

subscripts  r  and  a,  respectively.  Because 

■r  =  xcosP-ysinfi  (7-11) 

n  =  xsinv+ycos) 

It  1 


j5=5'+r°. 

{ n  =  n ' 


and  <  ,  where  r  is  the  radius  of  the  turn,  we  have 

I - -  i  o 


5 ' +rQ  =  xcos9-ysin9 
n  1  =  xsin9-t-ycos9 

(7-12) 

£ '  =  xcos9-ysin9-rQ 
n'  =  xsin9+ycos9 


which  is  the  coordinate  transformation  relation  between 
U  '  ,n  '  )  and  (x,y) . 

The  distance  between  a  target  point  at  (£  ,n  )  and  the 

cl  5L 

radar  at  (Sr,n  )  is 

s(£  ,n  ,£;  ,n_)  =  [  (£-£_)  2+  (n,-nr)  2]^ 

d  d  L  L  d  L  dl. 

=  U2+n2-2£  E  -2n  n  +C?+n?l% 

L  L  di  didd 
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(7-13) 


A  2  2  7' 

where  R  =  (x  +y  )  2  is  the  distance  between  the  center  of 
o  r  r 

the  turn  and  the  radar.  In  Eq.  (7-13)  we  have  assumed  that 


the  radius  of  the  turn  is  small  compared  to  Rq.  The  Doppler 
frequency  of  the  point  target  at  (£  , n  )  is  then 

d  cL 
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fd^a' 


na)  =  X 


ds  (£. 


,r|a'^r  ,nr^ 


dt 


2  ds  d9 
X  d0  dt 


(7-14) 
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ds 
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,  d[L (x  cos9-yrsin6)+n  (x  sin0+y  cosQ)] 

—  di  i  a.  IT  IT 


de 
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(7-15) 


Let 


(£-,,0,,)  =  r  (cos0  ,sin0  ) 

d  d  d  d  d 


and 


(^r'nr)  =  rr  (cos0r,sin0r)  =  Rq  (cos0r , sin9r) 

where  r  and  0  are  the  polar  coordinates  with  respect  to  the 
coordinates  (£,n)  with  origin  at  0.  Then 


3?  =  IT  ra‘Rosin(9a-0r)  =  ra’ sin (0a"9r)  (7“16) 

o 

Equation  (7-15) ,  which  assumes  that  r  <<  R  as  in 

o  o 

—  d  s 

Eq.  (7-13),  has  the  following  interpretation:  g-g-  is  propor¬ 
tional  to  r  ,  the  radial  distance  of  the  target  point  from 

d 

origin  0  and  is  independnent  of  0,  modulated  by  a  sinusoidal 

function  of  0-0  .  While  0  is  independent  of  0 ,  0  is  a 
a  r  a  r 

linear  function  of  9.  Hence  each  target  point  has  a  sinu¬ 
soidal  range  history.  In  the  extreme  case  when  rQ  =  0,  this 
would  be  the  same  as  a  turntable  geometry.  The  effect  of 
0^  is  to  make  the  family  of  constant  doppler  lines  rotated 


164 


by  an  angle  0  .  Note  that  0=0  for  the  target  center  and 

L  cl 

1  =  -rosinSr  <7'17> 

is  the  range  velocity  of  the  target  center.  Obviously 
II  =  0  at  0r  =  0°  and  180°  corresponding  to  broadside 
positions  and  is  maximum  in  absolute  value  corresponding  to 
head-on  or  tail-on  positions. 

In  the  case  rQ  is  not  small  compared  to  Rq  one  can 
still,  in  principle,  decompose  s (£  ,n  ,  £  ,n  )  into  two  parts, 

cl  3  IT  r 

one  consisting  of  the  pure  target  center  motion  and  the 
other  a  rotating  of  the  target  as  seen  from  the  radar,  to 
separate  the  relevant  imaging  information  from  the  irrele¬ 
vant  trajectory  information  using  the  precise  coordinate 
relation  in  Eq.  (7-12).  However,  derivation  of  the  exact 
expression  is  very  complicated  and  tedious  so  we  decided  not 
to  elaborate.  The  comment  that  is  necessary  is  that  the 
effective  target  rotation  rate  which  determines  the  azimuth 
resolution  is  the  sum  of  the  aspect  rate  due  to  the  turn 
itself  (thus  is  uniform)  and  the  aspect  angle  rate  of  the 
aircraft  center  as  seen  from  the  radar.  If  the  maximum  rate 
of  target  center  aspect  change  which  occurs  at  broadside  is 
much  smaller  than  the  uniform  turn  rate  ,  the  effective 
aspect  angle  rate  can  be  approximated  by  a  constant,  assum¬ 
ing  constant  radius  and  velocity  of  the  turn. 

As  in  Interval  1,  before  applying  a  Fourier  transform 
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to  the  data  returns  with  small  aspect  angle  change*,  the 
motion  of  the  center  of  target,  which  is  much  more  compli¬ 
cated  in  interval  2,  has  to  be  compensated.  This  will  be 
discussed  in  subsequent  sections. 

7 . 3  Preprocessing 

For  most  practical  purposes,  radar  imaging  systems 

which  determine  the  relation  between  the  data  returns  and 
the  reflectivities  of  the  target  can  be  considered  linear 
[7-2]  and  the  system  classification  method  developed  in 
Chapter  2  [7-4]  can  be  used  to  decide  ways  to  reconstruct 
the  reflectivities  directly  from  the  raw  data.  This  situa¬ 
tion  is  depicted  in  Fig.  7-5.  The  data  return  g(x,y)  is  a 
linear  transformation  of  target  reflectivity  function  f  { £ , n ) 
through  the  radar  signal  radiation  and  the  echo  reception. 
For  ease  of  presentation  we  will  assume  that  both  g  and  f 
in  Fig.  7-5  are  discrete  so  that  the  system  can  be  repre¬ 
sented  by  a  matrix  [H]  and  g  and  f  by  vectors  as  in 
Chapter  2.  Depending  on  the  waveforms  of  transmitted 
signals,  (e.g. ,  short  pulse,  linear  FM  pulse,  or  step- 
frequency  waveforms)  and  the  imaging  geometries  (e.g., 
shape  and  size  of  target,  direction  of  relative  motion, 
resolution  required,  etc.),  the  radar  imaging  systems 
represent  a  wide  spectrum  of  the  cases  in  [7-4].  Once  the 
relation  [H]  between  the  reflectivity  and  data  is 

Assuming  no  range  curvature  problem,  as  discussed  later. 
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(precisely)  decided  by  the  flight  or  radar  data,  a  straight¬ 
forward  reconstruction  of  f  and  g  can  be  achieved  by 
applying  the  pseudoinverse  of  [H]  to  g  yielding  a  minimum 
square  error  reconstruction.  Methods  of  taking  [H]  are 
readily  provided  by  each  of  the  nine  cases  in  Chapter  2. 

The  above  reconstruction  scheme,  although  straight¬ 
forward  in  theory,  usually  involves  a  great  deal  of  compu¬ 
tation  because  of  the  complexity  of  [H] .  In  the  worst 
case,  one  would  expect  to  resort  to  a  full  singular  value 
decomposition  (SVD)  to  find  IH]  1.  Certainly  a  decompo¬ 
sition  of  [H]  such  that  the  structure  of  the  imaging  geome¬ 
try  can  be  better  utilized  would  warrant  the  efforts  in 
many  cases. 

A  perceivable  way  to  accomplish  this  is  to  do  some 
preprocessing  upon  the  raw  data  such  that  the  resultant 
data  has  a  much  simplified  relation  to  the  reflectivity 
than  the  raw  data  itself.  Diagrammatically ,  [H]  can  be 

replaced  by  a  cascaded  system  of  [H^]  and  [H21  as  in 
Fig.  7-6  or  Fig.  2-2  (d)  and  f  can  be  estimated  by  multiply¬ 
ing  [H2]  followed  by  [ H ^ ^ ,  to  g  with  the  hope  that 
[ H ^ ]  would  be  so  simplified  in  structure  or  so  small  in 
size  compared  to  [H]  that  the  extra  effort  on  obtaining 
[H2 ]  1  would  be  warranted.  For  this  purpose  [H2]  *  is 
called  preprocessing.  Examples  of  [H2]  *  are:  range 
alignment,  presumming,  de-chirping,  and  general  motion 
compensation.  Some  of  them  will  be  discussed  in  the 
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target  reflectivity 
f  (  5,  n) 


linear 
radar 
system  [H] 


radar  return 

g(x.y) 


g  =  [h]  f 

[H] :  Point  spread  function  matrix  (PSF) 

[H][H]t:  correlation  matrix 

-  determines  systems  degrees  of 
freedom  (DOF) . 


Fig.  7-5.  Linear  radar  imaging  system 


target 


[^3  Preprocessing 

-  Range  alignment 

-  Range  curvature  elimination 

-  Motion  Compensation 

-  Presuming 

[H^J  ^ :  Ideally  a  Fourier  transform 


Fig.  7-6.  Decomposition  of  [H]  in  Fig.  7-5 


following  sections. 


7 . 4  Range  Curvature  and  Range  Bin  Alignment 

In  general,  the  radar  return  of  the  signal  pulse  from 
the  target  provides  the  range  information  while  the  history 
of  the  returns  along  some  range  bin  provide  azimuthal 
information.  These  two  sources  of  information  could  have 
been  coupled  such  that  a  separable  or  even  separate  pro¬ 
cessing  [7-4]  would  not  be  adequate  to  recover  the  infor¬ 
mation  to  the  extent  of  accuracy  one  pursues.  There  are  two 
major  sources  of  non-separability  in  the  radar  imaging 
system:  range  curvature  and  range  misalignment.  We  will 
describe  the  phenomena  and  propose  methods  to  avoid  or 
correct  them. 

7.4.1  Range  curvature 

A  single  radar  pulse  return  contains  the  information 
about  the  surfaces  or  lines  whose  points  are  equi-distant 
from  the  radar  transmitter  or  receiver.  These  surfaces 
or  lines  can  be  resolved  by  the  timing  (for  short  pulse)  or 
range  compression  (for  long  duration  linear  FM-like  pulse) 
or  even  synthetic  range  profiles  (for  MFS  pulses)  tech¬ 
niques.  Since  the  range  direction  has  been  compressed  and 
resolved  in  our  source  data,  the  simplest  way  to  resolve 
the  azimuth  would  be  to  do  one-dimensional  processing  along 
cross  range  direction.  This  requires  that  each  particular 
point  has  contribution  to  only  those  range  bins  which  are 


aligned  for  azimuthal  processing.  Such  is  the  case  for  low 
or  medium  resolution  SAR  imaging  with  aligned  returns.  As 
the  resolution  requirement  becomes  greater  and  greater  one 
is  usually  forced  to  reduce  the  range  bin  width  and/or  to 
increase  the  azimuthal  interval  over  which  the  data  are  to 
be  processed  coherently.  Both  of  these  would  eventually 
create  range  curvature  problem  since  the  surfaces  of 
constant  range  as  mapped  on  the  target  move  further  away 
as  the  relative  motion  between  the  radar  and  the  target 
continues. 

There  are  two  effects  as  one  increases  the  azimuthal 
coherence  interval:  increased  bandwidth  and  range  curvature. 
In  addition,  the  amount  of  computation  increases  with  the 
signal  bandwidth.  Assuming  separable  processing  there  is  a 
tradeoff  as  to  how  far  one  should  increase  the  coherence 
interval  without  too  serious  a  range  curvature  problem. 

There  have  been  few  quantitative  criteria  proposed  for 
optimizing  this  tradeoff,  especially  when  the  computational 
factor  is  also  considered. 

It  is,  however,  much  easier  to  find  the  critical 
condition  under  which  the  range  curvature  starts  to  occur. 
This  is  whenever  any  target  point  of  interest  has  a  range 
variation  of  more  than,  say,  1  range  bin  width,  so  that 
the  point  has  contribution  to  more  than  1  range  bin.  This 
condition  can  be  derived  from  the  geometry  of  imaging,  as  in 
Chapters  5  and  6. 
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Since  range  curvature  problems  are  introduced  by 
crossing  the  range  bins,  ideally  one  can  avoid  it  by 
increasing  the  range  bin  width.  This  means  sacrificing 
range  resolution  for  the  azimuthal  resolution  in  the  case  of 
separable  processing.  It  is  not  true,  though,  that  the 
range  curvature  limits  the  width  of  coherent  azimuth 
processing  available.  In  fact,  in  the  range  curvature 
situation  one  can  do  some  limited  compensation  by  the 
techniques  described  in  [7-5, 7-6]  or  even  full  compensation 
by  resorting  to  a  non-separable  model  for  the  imaging 
system  [7-4]  and  relying  on  singular  value  decomposition 
(SVD)  techniques.  However,  all  of  our  experiments  in  the 
following  chapter  assume  separable  processing  for  ease  of 
computation  and  implementation. 

7.4.2  Range  alignment 

In  addition  to  the  range  curvature,  there  is  another 
problem  which  hinders  the  separability  of  the  processing: 
range  misalignment.  As  described  before,  azimuthal 
processing  must  operate  upon  the  returns  from  target 
points  at  equal  range.  Thus  precise  timing  or  other 
schemes  on  returns  of  individual  pulses  to  insure  correct 
range  bin  alignment  is  of  ultimate  importance  to  a  single 
azimuth  processing. 

In  obtaining  the  data  from  our  radar  system,  range 
tracking  is  provided  by  a  Poly/Kalman  estimator  which 
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tries  to  lock  the  first  strong  peak  of  each  pulse  return 
onto  a  specific  range  bin.  For  example,  if  the  target 
point  closest  to  the  radar  is  the  wing  tip,  then  the  wing 
tip  returns  of  different  pulses  hopefully  will  be  locked  in 
the  same  range  bins.  Because  of  scintillation  of  the 
reflectivities,  this  range  locking  method  is  not  always 
reliable  and  misalignment  occurs  from  time  to  time.  Two 
automatic  algorithms  are  proposed  to  realign  the  data, 
one  being  in  the  frequency  domain  and  the  other  in  the 
spatial  domain. 


7. 4. 2.1  spatial  domain  re-alignment 

Let  f  (r)  and  f  (r)  be  the  high  ranqe  reso 
Z1  t2 


lution 


( HRR)  returns  (or  our  source  data)  from  adjacent  pulses 
where  ~  the  Pu^-Se  repetition  interval  (PRI) 

and  r  is  the  recorded  range.  Because  of  the  tiny  aspect 
angle  change  in  one  PRI,  if  we  consider  only  the  magnitude 
of  the  returns,  then 


m  (r+Ar)  «  m  (r) ,  where  m  (r)  =  |f  ( r ) |  (7-18) 

Z1  z2  l  l 

for  some  Ar,  the  amount  of  misalignment  which  we  would 

like  to  estimate.  Define  a  correlation  function  between 

the  two  waveforms  m  (r)  and  m.  (r) : 

Z1  z2 


R(s)  = 


m  (r)m  (r-s)dr 
;  j  -co  z2 _ 

[|  m^(r)dr|  m^(r)dr] 


(7-19) 


Then  because  m.  (r+Ar)  =  m  (r) ,  from  the  Schwartz  inequality 
Z1  Z2 

we  have  that  R(s)  will  be  maximal  at  s  =  Ar  and  the  amount 
of  misalignment  can  be  determined  by  a  subsequent  peak 
detection  on  R(s).  It  is  observed  that  the  denominator 
of  Eq.  (7-19)  is  independent  of  R(s)  and  can  be  dropped 
without  affecting  the  peak  location.  Thus  we  could  as  well 
use 

r°o 

R' (s)  =  m  (r)m  (r-s)dr  (7-20) 

J  t2 

—  OO 

which  is  a  straight  convolution  relation. 

7. 4. 2. 2  frequency  domain  re-alignment 

There  are  three  differences  between  f  (r)  and  f  (r) , 

C1  c2 

the  returns  from  adjacent  pulses.  One  is  that  the  range 
bins  have  been  shifted  relatively.  The  other  two  are  an 
overall  constant  phase  change  due  to  the  motion  of  the 
target  center  and  a  phase  variation  due  to  the  effective 
rotation  of  the  target.  The  last  phase  variation  can  be 
negligible  if  the  sampling  rate  in  azimuth  direction  (i.e., 
the  PRF )  is  large  enough.  In  that  case 

f  (r)  =  <j>  (t.  ,  t9)  f .  (r-Ar)  (7-21) 

tl  1  ^  z2 

where  ^>  ( t ^ ,  t^ )  is  a  phase  factor  independent  of  r  and  again 
Ar  is  the  quantity  to  be  estimated. 

If  one  takes  Fourier  transforms  of  both  sides  of 
Eq.  (7-21)  with  respect  to  r,  then 
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r¥  (f)  =  «D(t, 

zi  1 


Let  <t>  (f)  and  4>  (f) 

C1  Z2 

?  (f) ,  respectively,  then 

Z2 


t9)J  (f )  exp  (- j2rrf  •  Ar)  (7-22) 
r2 

be  the  phase  part  of  !F  (f)  and 


A4>  (f)  =  <P  (f)-4>  (f) 

Z1'Z2  Z1  Z2 


$  (t^  ,  t2  )  exp  (- j27TfAr ) 

(7-23) 


and 

&  *{A<b  (f)  }  =  S  1  { 4>  (t  ,  t0)  exp  (- j2TTfAr)  } 

1 1  /  t 1  z 


=  0»(t1,t2)6(r-Ar)  (7-24) 

hence 

I^Uo,  „  (f)  }  I  =  6  (r-Ar)  (7-25) 

tl't2 


The  above  analysis  suggests  the  following  way  to 
estimate  Ar:  Take  the  Fourier  transform  of  each  pulse 
return  and  keep  only  the  phase  part.  Find  the  phase  differ¬ 
ence  of  adjacent  signatures.  Finally,  take  the  inverse 
Fourier  transform  of  that  phase  d':  “ference  and  keep  only  the 
magnitude  part.  The  location  of  the  peak  along  the  final 
range  direction  gives  the  estimate  of  range  bin  shift 
required  to  realign  the  range  bins. 

Clearly,  the  above  two  re-alignment  schemes  are  equiva¬ 
lent  in  the  noise-free  situation.  Under  noisy  condition, 
they  would  have  different  performances.  However,  a  detailed 
analysis  is  not  attempted  in  this  work. 


7 . 5  Motion  Compensation 
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As  described  in  Section  7.2,  there  are  two  components 
of  the  target  motion:  the  motion  of  the  target  center 
relative  to  the  radar  and  that  of  different  target  points 
relative  to  the  target  center  as  viewed  from  the  radar. 

Only  the  latter  contributes  to  the  imaging  ability  of  the 
radar,  and  the  relation  between  the  latter  phase  variation 
and  the  target  reflectivity  is  a  simple  Fourier  transforma¬ 
tion  in  the  azimuth  direction.  Thus,  a  motion  compensation 
[ H2  ]  which  removes  the  effect  of  the  motion  of  the  target 

center  is  highly  desirable  as  a  preprocessing.  Equivalently 
this  is  to  find  the  trajectory  of  the  target  center  and 
remove  its  effect  from  the  raw  data. 

Two  schemes  of  such  a  motion  compensation  are  proposed 
for  our  imaging. 

First  scheme:  The  flight  path  of  the  target  center 
can  be  inferred  from  the  timing  of  the  pulse  returns.  For 
example,  in  the  first  interval  after  the  flight  path  has 
been  decided  to  be  a  straight  line  and  the  azimuth  angle 
determined,  the  result  of  the  analysis  in  Section  7.2  can  be 
applied  to  determine  the  coefficients  of  the  quadratic  and 
other  higher-order  phases  and  remove  the  xr  effects  to 
leave  only  the  phases  associated  with  x  which  is  re] 2vant 
to  imaging. 

Second  scheme:  Since  the  trajectory  of  a  single  target 
point  is  very  close  to  that  of  the  target  center,  the 
returns  from  that  point,  if  available,  can  as  well  be  used 


as  a  reference  to  compensate  for  the  targer  center  motion. 

In  fact,  this  is  equivalent  to  considering  this  target 
point  as  the  rotation  center  of  the  target.  The  phases  of 
this  reference  point  as  a  function  of  azimuthal  signatures 
can  then  be  subtracted  from  those  of  all  range  bins  at  the 
corresponding  signatures.  Care  should  be  exercised  to 
assure  two  things  in  this  process:  first,  the  size  of  the 
reference  point  must  be  small  enough.  This  is  because  the 
size  of  the  reference  point  decides  the  best  possible 
azimuthal  resolution.  Second,  for  each  signature,  the 
reference  range  bin  must  correspond  to  the  reference  point. 
This  requires  range  alignment  as  described  before.  In 
addition,  reference  point  must  be  rigid  on  the  target. 

There  are  situations  where  other  forms  of  motion 
compensation  should  or  could  be  applied:  For  example,  in 
many  cases  mathematical  simplicity  of  the  imaging  equation 
can  be  guaranteed  only  if  the  relative  motion  between  the 
radar  and  the  target  is  very  simple.  In  the  case  where  the 
flight  speed,  for  example,  is  varying  due  to  the  weather  or 
maneuvering  an  interpolation  of  the  azimuthal  data  can  be 
applied  to  obtain  equally  spaced  samples.  Alternatively,  a 
varying  PRF  which  follows  the  speed  of  the  flight  can  be  an 
effective  tool  in  assuring  a  simple  imaging  equation.  On 
the  other  hand,  some  times  even  if  the  flight  speed  is  a 
constant,  the  interpolation  or  varying  PRF  schemes  may  still 
be  needed.  For  example,  in  the  first  interval  portion  in 
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d  0 

Section  7 . 2 , —  in  Eq .  (7-9)  is  not  a  constant  due  to  the 

ratio  of  nr/R Q,  and  thus  the  Fourier  transform  relation 
would  not  be  exact  if  there  were  no  compensation. 

7 . 6  Oversampling  and  Presummina 

Usually  the  radar  imaging  system  is  oversampled  in  the 
azimuth  direction  because  of  a  too  high  PRF .  The  purpose 
of  presumming  is  to  remove  the  oversampling  redundancy. 

In  the  case  of  terrain  imaging,  over-sampling  could  be  a 
result  of  not  processing  the  whole  antenna  illumination 
pattern  along  the  azimuth  direction.  In  that  case,  the 
azimuth  pattern  width  utilized  or  coherently  processed 
determines  the  resolution  of  the  image  which  in  turn  deter¬ 
mines  the  amount  of  data  redundancy.  In  our  aircraft 
imaging,  the  situation  is  subtly  different.  Here  the 
azimuthal  width  of  the  aircraft  is  so  small  that  we  always 
try  to  make  full  use  of  the  maximum  width  of  the  effective 
radar  illumination  pattern,  which  is  the  azimuthal  length 
of  the  aircraft  itself.  Under  this  condition  the  PRF 
required  is  decided  by  the  azimuth  dimension  on  the  air¬ 
craft  and  the  azimuth  resolution  again  by  the  interval 
coherently  processed.  Thus,  assuming  other  parameters 
fixed,  a  larger  aircraft  would  require  a  higher  minimum 
PRF  to  insure  that  no  aliasing  would  occur  in  the  final 
images.  Also,  since  the  effective  antenna  illumination 
(i.e.,  the  overall  aircraft  azimuthal  length)  is  independent 
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of  the  wavelength,  X,  the  minimum  PRF  or  the  resolution  in 
the  aircraft-imaging  case  would  be  function  of  X.  This  is 
in  contrast  to  the  ground  terrain  imaging  cases  where  the 
full  antenna  illumination  pattern  width,  which  is  propor¬ 
tional  to  X ,  is  to  be  used  so  that  the  resultant  resolution 
is  independent  of  X  because  of  a  cancelling  effect  [7-1,2] . 

Let  f  be  the  carrier  frequency  and  L  the  length  of  the 
aircraft  along  the  direction  normal  to  the  LOS  and  on  the 
imaging  plane,  as  shown  in  Fig.  7-7.  Let  A0  be  the  orien¬ 
tation  change  of  the  target  between  two  adjacent  pulses 
as  observed  from  the  radar,  then  the  azimuthal  frequency 
change  will  be 

Af  =  f  A0  (7-26) 

Z  U 

This  means  a  sampling  interval  of  Afz  in  the  azimuthal 

frequency  domain,  which  implies  a  non-ambiguous  azimuthal 

1  c 

time  interval  of  j-j—  ,  or  spatial  interval  of  2Kt~>  ;''n 

z  z 

accordance  with  the  sampling  theorem. 

The  oversampling  factor  o  can  be  determined  by 

o  «  j±-/L  (7-27) 

z 

To  remove  the  oversampling  and  leave  minimum  useful 
data,  a  coherent  low  pass  filter  followed  by  sampling  at  a 
correspondingly  low  rate  should  be  applied  along  the  azimuth 
following  motion  compensation.  The  effect  of  low-pass 
filtering  is  to  remove  the  high  frequency  noise  which  other- 
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Fig.  7-7.  Determination  of  azimuth 
oversampling  factor 
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wise  would  appear  in  the  resultant  image. 

7 . 7  Conclusions 

In  this  chapter  a  radar  imaging  system  more  practical 
than  the  previous  RAT-SCAT  is  researched.  Although  the 
baseline  coherent  processing  principle,  i.e.,  Fourier 
transforms  along  azimuthal  direction ,  is  the  same  as  that 
in  RAT  SCAT,  the  data  has  to  undergo  much  more  processing 
before  a  final  image  can  be  realized.  Motion  compensation, 
range  realignment,  and  presuming  are  additional  requirements 
compared  to  the  RAT  SCAT  case.  The  study  of  RAT  SCAT, 
which  is  in  a  well-control led  environment  is  of  much  help 
to  the  analysis  and  processing  in  this  and  the  following 
chapters. 
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Chapter  8 

EXPERIMENTAL  RESULTS  OF  RADAR  IMAGING 


WITH  TARGET  MOTION 

8. 1  Introduction 

Unlike  the  RAT- SCAT  case,  the  in-flight  data  represents 
complicated  flying  geometries  and  requires  much  more 
analysis  and  processing  to  reach  final  radar  images.  A 
series  of  processing  techniques  developed  from  the  previous 
chapter  will  be  presented  along  with  their  experimental 
results . 

The  mode  in  which  the  radar  operated  and  our  source 
data  was  acquired  was  a  wide  band  high  resolution  mode. 

The  transmitted  pulses  were  linear  FM  signals  whose  returns 
were  then  compressed  using  matched  filtering  techniques  in 
the  radar  receiver  and  aligned  side  by  side  by  using  a 
Poly/Kalman  centroid  estimator. 

8. 2  Experimental  Results  -  First  Interval 

A  condensed  overall  view  of  magnitude  oart  of  first 

interval  data  is  shown  in  Fig.  8-1  in  which  each  row 

corresponds  to  the  log  magnitude  of  the  compressed  return  of 

a  single  pulse.  Only  every  16th  signature  is  shown  in  this 

figure.  Recalling  that  this  interval  represents  the  radar 

returns  when  the  target  aircraft  was  flying  toward  a  broad- 
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side  position  (Fig.  7-1) ,  we  presume  that  the  first  high- 
intensity  range  bin  corresponds  to  the  left  wing  tip  and 
the  next  distinct  strong  returns  are  from  the  fuselage  and 
nose.  Note  that  the  radar  is  to  the  left. 

Then  it  can  be  perceived  from  Fig.  8-1  that  the  fuse¬ 
lage  is  at  a  greater  and  greater  range  distance  away  from 
the  wing  tip  along  the  signatures  as  a  result  of  closing- 
to-broadside  during  flight.  Also  observe  that  while  most 
portions  of  Fig.  8-1  seem  pretty  well  range-aligned,  other 
portions  do  need  re-alignment  before  a  separable  processino 
can  be  implemented. 

To  present  the  data  in  detail  all  of  the  first  512 
signatures  are  displayed  in  Fig.  8-2.  The  phase  image 
(Fig.  8-2 (b) )  indicates  clearly  that  the  target  points 
probably  lie  in  range  bin  number  50  to  200,  where  a  strong 
structure  of  phase  relationships  appears  as  a  result  of  the 
coherent  radar  pulsing.  This  is  also  shown  in  the  log 
magnitude  picture  Fig.8-2a,  although  with  less  clarity. 
There  is  a  transient  region  where  the  strength  of  the 
returns  decreases  gradually  with  the  range  or  time.  This 
is  conjectured  to  be  a  result  of  multiple  reflections  on 
the  target  which  take  more  than  before  re-radiating  to  the 
radar  receiver. 

To  investigate  further  the  behavior  of  the  returns,  we 
kept  only  the  regions  of  stronq  signal  returns  and  a 
sequence  of  80  second  data  or  8192  signatures  is  shown  in 
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Fig.  8-3  with  both  log  magnitude  and  corresponding  phase. 
Observe  the  quadratic-like  phases  along  the  flight  direction 
due  to  the  flight  geometry,  as  analyzed  in  Section  7.2. 

Since  the  radar  receiver  has  range  compressed  the 
signal  returns,  ideally  we  will  only  need  to  perform  some 
azimuthal  processing.  For  convenience  we  transpose  the 
data  so  that  the  horizontal  direction  now  denotes  the 
signature  or  along-flight  direction.  Figure  8-4  shows  the 
log  magnitude  and  phase  of  typical  signatures  (signature 
8001  to  8512) .  To  remove  the  quadratic  phases  from 
Fig.  8-4 (b) ,  three  options  exist,  the  first  two  being 
similar.  They  are  described  in  subsections  8. 2. 1-3: 

8.2.1  Linear  fitting  the  phase  differences  and  subse¬ 
quent  integration 

Since  the  differentiation  of  quadratic  phases  is  linear 
phases,  a  linear  fit  to  the  phase  differences  can  be  applied 
to  determine  the  quadratic  phase  curvature.  Figure  8-5 (a) 
shows  the  azimuthal  phase  difference  of  Fig.  8-4  (b) .  Note 
that  except  for  the  phase  wrap-around  in  the  right  half 
portion  and  the  small  variation  due  to  oversampling  and 
noise,  the  intensity  used  to  linearly  encode  the  phase 
between  -tt  to  tt  looks  quite  linear.  However,  before  a 
successful  linear  fit  can  be  obtained,  the  phase-wrap 
problem  has  to  be  solved  and  this  is  usually  not  a  very 
easy  task.  In  fact,  it  is  because  of  the  raoid  phase 
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modulus  phenomena  in  Fig.  8-4  (b)  that  causes  direct  unwrap¬ 
ping  Fig.  8-4  (b)  extremely  difficult.  We  used  a  simple- 
minded  scheme  to  unwrap  the  phase  differences  of  Fig.  8-5  (a) 
to  get  Fig.  8-5 (b) ,  from  which  the  linear  portion  of  phase 
variation  was  estimated  and  removed  to  leave  Fig.  8-5  (c). 
Since  Fig.  8-5  (c)  is  still  in  the  differentiation  domain, 
an  integration  brings  it  back  to  the  azimuthal  phase 
domain,  as  depicted  in  Fig.  8-5  (d). 

8.2.2  Linear  fitting  the  phase  difference  and  quadra¬ 
tic  subtraction 

An  alternative  to  applying  the  estimated  linear-phase- 
difference  is  to  subtract  the  estimated  quadratic  phase 
(from  integration  of  estimated  linear  phase  difference) 
from  Fig.  8-4 (b)  directly.  The  result  is  shown  in 
Fig.  8-5  (e) . 

Magnitudes  of  azimuthal  Fourier  transforms  of 
Fig.  8-5(d)  and  8-5(e)  are  shown  in  Figs.  8-6(a)  and  8-6(b), 
which  are  very  similar  visually. 

8.2.3  Target  point  referencing 

The  above  two  schemes  of  removing  phase  variation  due  to 
target  center  motion  are  based  on  an  assumption  that  the 
flight  path  is  relatively  straight  during  the  coherence 
time.  In  other  cases  where  the  range  trajectory  is  more 
complicated  than  a  low-order  polynomial  curve,  the  above 
schemes  are  expected  to  be  difficult.  Another  motion  corn- 
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Fig.  8-4.  Signature  number  8001-8512 

(range  direction  is  vertical) 


a)  Phase  difference  of 
Fig.  4b) 


b)  Unwraeped  version 
of  a ) 


Fig.  8-5. 


Motion  compensation  on  Fig.  8-4 (b)  usine 
linear  curve  fitting 


c)  Linear  phase  removed 
from  b) 


d)  Phase  integration 
of  c) 


e)  Quadratic  phase  removed 
directly  from  4b) 


Fig.  8-5.  (Continued) 
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pensation  scheire  somewhat  independent  of  the  flight  geometry 
and  simple  in  implementation  is  to  use  the  signal  returns 
from  a  reference  point  to  estimate  the  history  of  the  flight 
range  trajectory.  This  single  point  can  be  thought  of  as 
the  center  of  rotation  of  the  target  and  its  phases  subtrac¬ 
ted  from  all  other  range  bins  to  leave  only  the  phase 
histories  of  all  target  points  relative  to  this  reference 
point.  This  was,  in  fact,  the  techniques  used  in  subsequent 
imaging. 

Since  the  azimuthal  oversampling  factor  has  been 
determined  to  be  greater  than  50,  the  data  after  quadratic 
phase  compensation  can  be  reduced  by  a  factor  of  32  before 
Fourier  transformations  are  applied.  The  result  is  shown  in 
Fig.  8-6  (c).  A  comparison  of  Figs.  8-6  (a)  and  8-6  (c) 
confirms  the  validity  of  the  coherent  presumming.  Note  that 
in  Figs.  8-6  (a),  (b)  and  (c)  the  azimuthal  bin  width  is 

much  wider  than  the  range  bin  width  and  a  subsequent  inter¬ 
polation  has  to  be  done  to  properly  scale  the  images. 

Figure  8-7  is  a  series  of  processed  aircraft  images 
using  the  reference-point  scheme.  Consecutive  pictures 
represent  abutting  2048  signatures  or  20-second  flight 
time.  The  images  have  been  linearly  interpolated  in  azimu¬ 
thal  to  give  the  same  range  and  azimuthal  bin  width  such 
that  the  images  are  properly  scaled.  Visually  Fig.  8-7 (d) 
is  the  best  probably  due  to  the  best  range  alignment  of  the 
data  in  that  time  interval. 
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a)  Magnitude  of  Fourier 
transform  of  5d) 


b)  Magnitude  of  Fourier 
transform  of  5e) 


c)  Fourier  transform  of 

collapsed  version  of  5e) 


Fig.  8-6.  Cross-range  processing  Fig.  8-5 


a)  1st  20  seconds  or  2048  b)  2nd  20  seconds 

signatures  (~2 . 5  aspect 
change) 


c)  3rd  20  seconds 


d)  4th  20  seconds 


1  o  2 


i 


Fig.  8-7.  Aircraft  radar  images  with  abutting 
20  second  coherence  times 


Ideally,  an  increase  in  coherence  time  should  be 
accompanied  with  an  equally  increased  amount  of  resolution. 
This  is  not  the  case  in  Fig.  8-8,  where  coherence  times  of 
40  and  80  seconds  are  processed.  The  conjecture  is  that 
the  range  curvature  and  range  misalignment  which  tend  to 
blur  the  images  outplay  the  coherence  time  increase.  As 
described  in  the  previous  chapter,  one  way  to  alleviate 
the  range  curvature  problem  is  to  use  larger  range  bin 
widths.  To  test  this,  we  used  the  same  parameters  as  in 
Fig.  8-8 (b)  except  the  data  in  range  dimension  were  reduced 
by  a  factor  of  two  by  a  coherent  collapsing.  The  result 
shown  in  Fig.  8-9  is  to  be  compared  with  Fig.  8-8 (b) . 

The  phase  variations  of  target  points  induced  by  the 
target  motion  is  the  key  to  the  coherent  radar  imaging. 

As  one  can  see  from  Fig.  8-2,  the  magnitude  of  the  radar 
returns  which  provide  only  range  information  are  very 
similar  from  pulse  to  pulse  and  represent  a  great  deal  of 
data  redundancy.  From  the  DOF  point  of  view  one  would  like 
to  have  approximately  equal  amount  of  input  and  output 
data.  Since  only  the  magnitude  parts  of  the  processed 
images  are  needed  for  displays,  it  is  conjectured  that 
the  phase  portion  of  data  alone  is  sufficient  to  give  an 
image  of  a  comparable  quality.  This  would  achieve  a  factor 
2:1  in  data  reduction.  Experimental  result  shown  in 
Fig.  8-10  seems  to  support  this  conjecture.  Intuitively 
speaking,  the  range  bins  corresponding  to  no  strongly 
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a)  1st  40  seconds  1.0  .'did  4  > 

( «  5  asoect  chance) 


Fig.  8-9.  2:1  rancie  col]  aped;  coherenc 

time  ~  40  seconds 


Fig.  8-10.  Same  portion  as  Fig.  R-7d) 
phase  information  only 


reflective  target  points  have  a  random-like  phase  and  are 
likely  to  spread  their  energy  over  the  spectrum  after  the 
Fourier  transform  in  the  azimuthal  direction  is  taken.  On 
the  other  hand,  target  points  of  strong  reflectivities  give 
highly  correlated  azimuthal  radar  returns,  resulting  in 
clusterings  of  energy  corresponding  to  different  azimuthal 
target  points  in  the  processed  Fourier  transform  domain. 

In  this  way,  the  magnitude  of  the  returns  do  not  play  an 
important  role  in  determining  in  which  range  bins  lie  the 
strong  target  points. 

8 . 3  Experimental  Results  -  Second  Interval 

The  first  8000  signatures  of  the  second  interval 
source  data  taken  when  the  airplane  was  making  a  standard 
left  turn  are  shown  in  Fig.  8-11  and  Fig.  8-12.  Unlike  the 
straight  flight  of  interval  1,  the  phase  plot  here  has  a 
changing  azimuthal  structure  due  to  the  turning  motion  of 
the  target,  which  creates  complicated  range  and  Doppler 
histories.  In  addition,  there  are  several  occasions  when 
the  range  bins  are  seriously  out  of  alignment.  The 
overall  view  of  Fig.  8-11  shows  the  changes  of  relative 
positions  of  the  nose,  fuselage  and  wing  tip  due  to  the 
turn.  A  portion  of  data  was  taken  when  the  airplane  was 
approximately  nosed  into  the  radar  and  a  series  of 
resultant  images  are  shown  in  Fig.  8-13  using  the  reference- 
point  technique  as  a  phase  compensator.  In  this  case  the 
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a)  1st  2.5  seconds  og  256  b)  2nd  2.5  seconds 

signatures  { %  4 . 5°  asnect 
change ) 


c)  3rd  2.5  seconds  d)  4th  2.5  seconds 


Fig.  8-13.  Aircraft  radar  images  with  abutting 
2.5  second  coherence  times  in  second 
interval 
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nose  tip  serves  as  a  very  good  reference  point  as  shown  by 
the  degree  of  sharpness  of  the  nose  in  these  images. 

Figure  8-14  shows  images  of  different  coherence  times.  In 
Fig.  8-14  (b)  the  shape  of  the  fuselage  has  been  clearly 
imaged.  A  coherence  interval  of  18°rotation  of  the  target 
seems  too  much  to  give  a  satisfactory  image  as  a  result  of 
overwhelming  range  curvature. 

The  spread  patterns  close  to  the  nose  are  due  to  the 
aircraft  radar  which  was  constantly  scanning  during  the 
flight,  presenting  an  object  of  changing  reflectivity  and 
violating  the  assumption  that  the  target  was  a  rigid  body 
in  the  imaging  system. 

8 . 4  Range  Re-alignment  Results 

As  is  evident  from  Figs.  8-11  and  8-12  the  radar 
breaks  range  lock  quite  often  during  the  turn  of  the  tarqet 
aircraft.  This  is  to  be  expected  as  different  scatterers 
from  the  aircraft  dominate  the  leading  return  of  the  radar 
reflection.  Naturally  when  the  radar  broke  lock  seriously, 
one  would  not  expect  to  be  able  to  image  without  re¬ 
alignment  processing.  Section  7.4  presented  a  theoretical 
discussion  on  such  re-alignment  procedures  and  this  section 
will  present  some  experimental  results. 

The  realignment  procedures  presented  in  the  last 
chapter  are  based  on  a  sequential  comparison  between 
adjacent  signatures.  Because  of  the  discrete  nature  of  the 
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a)  1st  5  seconds 

(as  90°  aspect  change) 


h)  2nd  5  seconds 


c)  1st  10  seconds 

(as  18°  aspect  change) 


Fig.  8-14.  Aircraft  images  with 

different  coherence  times 
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data,  each  comparison  could  produce  a  misalignment  error  by 
up  to  one  pixel  width.  Hence,  although  the  scheme  is 
successful  on  local  realignment,  the  accumulated  errors 
from  this  primitive  recursive  procedure  make  a  global 
alignment  very  difficult  to  produce  satisfactorily  corrected 
images . 

A  remedy  to  this  is  to  compare  each  signature  with  a 
reference  which  has  been  aligned  globally.  However,  this 
reference  has  to  bear  considerable  similarity  to  the 
signature  being  compared  to  avoid  large  quantity  error. 

From  these  requirements  and  considering  the  slowly  varying 
nature  of  the  successive  signatures,  we  propose  to  use  for 
reference  an  exponentially  weighted  sum  of  the  magnitudes 
of  previous  realigned  signatures.  The  algorithm  is  easily 
implementable  by  assigning  a  weight  to  the  newly  aligned 
signature  and  a  complement  weight  to  previous  reference  to 
form  the  current  reference.  The  weight  can  be  chosen 
according  to  how  slowly  the  range  profiles  change  along 
the  signature  assuming  the  profiles  resemble  a  "Markov 
process."  A  small  weight  corresponds  to  high  correlation 
of  the  process.  On  the  other  hand,  a  weight  of  1  corres¬ 
ponds  to  very  low  correlation  and  coincides  with  the 
theoretical  realignment  scheme  described  in  the  last 
chapter. 

We  applied  the  modified  realignment  algorithm  using 
exponential  weighting  for  the  signature  magnitude  correla- 
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tion  and  obtained  a  series  of  realigned  data  as  in 
Fig.  8-15.  Sequential  pictures  contains  512  signatures 
or  a  flight  time  of  5  seconds  each.  Note  that  the  first  16 
pictures  correspond  to  the  original  data  in  Fig.  8-12. 

A  weight  of  0.05  was  used  in  all  21  data  segments  which 

comprise  the  whole  second  interval  of  the  flight.  Except 
for  very  few  signatures  the  realignment  algorithm  seems  to 
work  very  well. 

The  realigned  data  was  processed  during  the  same 
reconstruction  procedure  before  and  a  series  of  21  images 
produced  as  in  Fig.  8-16  in  which  the  effective  rotation 
rate  of  the  aircraft  was  assumed  constant.  Most  of  the 
aircraft  images  are  recognizable  and  clearly  show  the 
turning.  Fig.  8-16  (a)  and  (b)  correspond  to  the  motion 
before  completing  the  starting  of  the  turn  and  are  thus  not 
square  pictures.  The  nose  and  left  wing  tip  crossed  in  and 
out  respectively  at  the  later  portion  of  the  data  segment 
(e)  as  shown  by  Fig.  8-15 (e)  and  by  the  change  of 
Fig.  8-16  (e)  from  Fig.  8-16  (f).  In  Fig.  8-16  (r)  the  image 
was  a  superposition  of  two  images  with  the  same  orientation 
and  with  the  wingtip  of  one  image  being  on  top  of  the  nose 
of  the  other.  This  is  due  to  the  equi-distance  of  those 
two  target  points  from  the  radar  at  that  particular  target 
orientation.  In  Figs.  8-16  (t)  and  8-16  (u)  the  wingtip  was 
chosen  as  the  reference  range  bins  for  motion  compensation 
whereby  the  range  realignment  procedure  aligned  the  data 
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Fig.  8-16.  (Continued) 
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with  respect  to  the  fuselage.  Since  the  fuselage  trajectory 
is  not  parallel  to  that  of  the  winqtip,  the  motion  compen¬ 
sation  did  not  work  and  the  image  is  blurred. 

8 . 5  Conclusions 

This  chapter  has  presented  the  experimental  results 
associated  with  imaging  a  moving  aircraft  target  from  a 
stationary  ground-based  radar.  Preprocessing  techniques 
for  range  curvature,  range  alignment,  motion  compensation, 
and  coherent  presumming  are  all  presented.  Gross  phase  and 
fine  tuning  phase  are  separated  for  imaging  the  target. 
Coherence  processing  intervals,  range  collapsing  and  range 
re-alignment  are  all  examined.  Two  data  sequences,  interval 
1  and  interval  2,  are  taken  for  real  radar  data.  Many 
parameters  are  tested  and  their  results  are  shown.  Each 
resultant  image  has  been  tested  against  the  physical  air¬ 
craft  flight  information  before  any  interpretations  are 
attempted.  Images  of  encouraging  quality  confirm  the 
feasibility  of  the  coherent  radar  imaging  of  the  flying 
aircraft. 
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Chapter  9 


SUMMARY  AND  EXTENT IONS 


9 . 1  Summary 

This  dissertation  has  presented  analysis  and  experi¬ 
mental  results  of  several  radar  imaging  systems.  As  a 
preliminary  research  of  this  work,  a  classi f ication  method 
on  two  dimensional  systems  based  on  the  structure  of  their 
point  spread  function  matrices  has  been  proposed  to  aid  in 
understandina  imaging  systems.  The  classif ication 
principle  is  universal  and  cases  apply  to  the  specialized 
field  of  radar  imaging.  Many  physical  examples  are  shown 
to  fit  in  most  of  the  cases  nicely  and  it  is  hoped  that  more 
systems  could  be  identified  with  some  of  the  cases  and 
possibly  even  new  imaging  systems  developed  based  on  the 
classified  system  structure  as  the  image  reconstruction 
techniques  advance.  We  have  demonstrated  that,  depending 
on  the  accuracy  one  pursues,  the  stripping  mode  of  SAR 
imaging  can  range  from  the  non-separate  1-D  invariant  case 
to  the  separable  2-D  invariant  case.  While  going  to  a 
simpler  system  model  saves  some  computation  time  and  storage 
requirement,  part  of  the  information  becomes  irretrievable 
as  a  result  of  approximation.  Analysis  on  the  tradeoff 


between  the  DOF  loss  and  the  simplification  of  the  system  PSF 
is  a  possible  extention  of  the  system  classification  work. 

The  radar  system  analysis  and  the  processing  techniques 
in  this  dissertation  assume  discrete  linear  imaging  systems 
with  the  target  being  a  rigid  body  of  constant  reflectivi¬ 
ties.  Physical  phenomena,  e.g.,  target  point  blocking, 
reflectivity  change  as  a  function  of  aspect  angle,  multiple 
reflections,  parts  moving  on  the  target,  etc.  if  completely 
incorporated  in  the  PSF's  would  make  the  actual  imaging 
system  fairly  nonlinear.  Another  nonlinearity  comes  from 
our  data  processing  in  which  both  the  radar  returns  and  the 
resultant  image  are  transformed  on  a  point  by  point  basis, 
e.g.,  taking  the  log  of  the  magnitude  parts,  to  reduce  the 
dynamic  range.  The  effects  of  those  nonlinearities  remain 
to  be  analyzed  in  the  extention  work. 

In  the  RAT-SCAT  system,  image  reconstruction  was 
assumed  separable  on  the  individual  data  section  following 
considerations  of  physical  limitations  and  fast-algorithm 
desirability.  Larger  sectors  (e.g.,  25.6°)  could  have  been 
processed  if  we  had  chosen  to  nonseparably  process  the  data 
in  which  a  constant  target  reflectivity  function  could  be 
assumed.  The  quality  of  the  final  RAT-SCAT  images  are 
limited  by  target  size,  wavelength,  bandwidth  of  the 
transmitted  signal  set,  fidelity  of  the  radar  receiver 
processing  unit,  background  reflectivity,  sidelobe  effect, 
and  ability  to  find  the  true  rotation  center  of  the  compo- 
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nent  images. 


While  the  RAT-SCAT  imaging  is  performed  under  a  well- 
controlled  environment,  the  imaging  geometry  in  Chapters  7 
and  8  is  more  realistic  in  that  a  full-sized  flying  aircraft 
target  is  imaged.  However,  because  of  the  increased  comple¬ 
xity  of  the  PSF  ,  a  method  of  system  decomposition  is  pro¬ 
posed  in  which  the  preprocessing  simplifies  the  relation¬ 
ship  between  the  transformed  data  and  the  object  function. 
Range  realignment,  range  curvature  compensation,  motion 
compensation,  and  presuming  are  examples  of  preprocessing 
performed.  After  the  preprocessing,  the  remaining  system  is 
(hopefully)  a  separable  one  and  fast  processing  techniques 
like  the  FFT  can  be  applied.  Side  lobe  effects  in  the  final 
images  can  be  greatly  reduced  by  incorporating  a  weighting 
window  upon  the  data  in  the  presumming  and  FFT  processings. 

In  both  the  RAT-SCAT  and  the  motion-induced-aperture 
imaging,  the  range  bin  width  and  azimuth  bin  width  are 
equalized  by  interpolation  techniques  to  make  square  images 
suitable  for  human  interpretation.  This  involves  calculating 
the  relative  magnitudes  of  the  range  bandwidth  which  is 
determined  by  the  bandwidth  of  the  transmitted  radar  signals 
and  the  azimuth  bandwidth  which  is  proportional  to  the 
amount  of  rotation  of  the  target. 

While  the  range  resolution  is  independent  of  the 
carrier  frequency  fg  of  the  radar  signals,  the  azimuth  reso¬ 
lution  is  a  function  of  fg.  In  fact,  it  is  the  large  fg 
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which  makes  imaging  possible  along  the  azimuth  direction 
with  only  a  few  degrees  of  the  target  rotation.  Recall  that 
although  the  PSF's  of  the  tomographic  system  have  a  very 
similar  mathematical  form  as  those  of  the  RAT-SCAT  system, 
it  needs  a  full  180°  projections  to  reconstruct  the  image 
because  of  a  zero  f^. 

The  final  radar  images  are  enhanced  to  reveal  more 
details.  Most  of  the  postprocessing  enhancement  techniques 
used  are  nonlinear  and  belong  to  the  immeasurable  art  of 
image  enhancement  and  image  analysis  [9-1].  Nevertheless 
the  understanding  of  the  characteristics  of  the  radar  imaging 
systems,  in  particular  the  cascaded  decomposition  of  the 
processing  procedure,  could  aid  in  the  advancement  of 
techniques  for  image  storage,  encoding,  enhancement, 

quantization,  and  restoration.  For  example,  knowing  the 
dynamic  range  limitation  of  the  film  and  the  Fourier  trans¬ 
formation  property  of  the  lens,  we  could  store  the  radar 
images  on  the  films  in  the  frequency  domain  to  spread  the 
energy  over  the  film  [9-2].  In  the  image  encoding  domain, 
Chen  [9-3]  has  developed  a  radar  image  coding  method  based 
on  dual  mode  interpolation  for  images  with  strong  targets 
and  hence  high  dynamic  range. 

9.2  Sonar  Imaging 

Finally,  a  potential  application  of  the  coherent 
radar  imaging  principle  to  other  systems  will  be 


discussed . 
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As  described  earlier,  the  two-dimensional  imaging 

systems  must  collect  a  set  of  one-dimensional  data  about  the 

target  reflectivity  or  transmission  from  different  aspect 

angle  of  the  object  for  image  reconstruction.  In  the  MFS 

system,  the  gathered  data  is  in  the  frequency  domain  and 

represents  points  at  particular  target  aspect  angles  and 

with  radii  being  the  radar  frequency  f  translated  into  the 

2f  C 

Q 

spatial  frequency  domain  or  — where  c  is  the  speed  of 

light  in  the  case  of  radar  imaging.  One  advantage  of  the 

large  fQ  lies  in  the  fact  that  it  places  the  data  points 

away  from  the  origin  of  the  frequency  coordinate  system 

making  the  data  points  relatively  distant  from  each  other 

and  less  redundant,  as  analyzed  in  Chapter  4.  The  range 
2Br 

bandwidth  is  — —  in  the  radar  imaging,  where  Br  is  the 
effective  bandwidth  of  the  transmitted  signals  in  Hz. 

Because  of  the  exceedingly  large  value  of  c,  a  very  large 
bandwidth  Br  of  the  signal  is  required  to  achieve  a  fine 
range  resolution  comparable  to  the  target  size.  For  example, 
bandwidths  larger  than  1024  MHz  are  required  to  achieve  a 
range  resolution  of  0.5  foot.  This  large  bandwidth  require¬ 
ment  has  set  many  restrictions  on  the  system  design  and 
limitations  on  its  overall  performance. 

In  the  sonar  system,  the  situation  is  quite  different. 
Here  the  propagation  speed  v  of  the  sound  in,  say  the 
vacuum,  is  only  1086  ft/sec  approximately.  The  disparity 
between  the  light  speed  and  sound  speed  is  in  rne  6th  order 
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of  10' s  power.  To  obtain  the  same  data  points  in  the 
frequency  domain  of  a  sonar  system  as  in  the  radar  system 
requires 


c  v 

where  fg  and  Bs  are  the  carrier  frequency  and  the  bandwidth 

of  the  sonar  wave  in  Hz,  respectively.  Thus  f  =  6.6  kHz 

and  Bs  =  565  Hz  in  the  sonar  system  would  be  compatible 

with  fA  =  6.6  kH  and  B  =  512  MHz  in  the  radar  system, 

0  z  r 

_  ^ 

a  required  bandwidth  only  10  of  that  of  the  radar  system. 
An  interpretation  of  the  above  analysis  is  as  follows:  the 
amount  of  information  about  the  target  inferrable  by  a 
single  pulse  is  given  by  the  time-bandwidth-product  of  the 
received  signal,  or  approximately  the  number  of  range  bins 
on  the  target.  Since  the  time  duration  of  the  received 
signal  is  in  proportion  to  the  time  on  target  (TOT) 
it  follows  that  for  a  fixed  number  of  TBP,  the  required 
bandwidth  is  higher  in  the  case  of  a  faster  propagation 
speed  or  a  smaller  TOT. 

The  above  analysis  is  based  purely  on  mathematical 
comparison.  Because  of  the  large  characteristic  differences 
between  the  electromagnetic  waves  and  sonar  waves,  there  are 
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many  physical  considerations  to  be  made  before  actually 
applying  the  SAR  techniques  to  the  sonar  imaging  system. 

For  example,  in  radar  imaging  systems,  it  is  usually  assumed 
that  during  the  transmission  and  receiving  of  a  single  pulse, 
the  target  is  stationary.  This  assumption  makes  imaging 
equations  relatively  simple  and  is  justified  by  the  high 
speed  of  the  light.  In  the  sonar  system  the  signal  dwell 
time  on  the  target  could  be  so  long  (because  of  a  low  v)  that 
the  target  might  rotate  by  a  significant  amount  during  TOT 
and  violate  the  stationarity  assumption.  Compensating 
target  motion  within  one  pulse  time  could  be  very  difficult. 
Also,  the  PRF  requirement  is  not  as  easily  met  because  of 
the  slow  v.  From  this  point  of  view,  sonar  system  should 
be  better  suitable  to  imaging  slowly  moving  targets  like 
the  submarines.  Another  major  disadvantage  with  the  sonar 
system  is  the  high  attenuation  rate  of  the  sound  wave  which 
tends  to  impose  a  requirement  for  very  high  power  trans¬ 
mitters  even  in  the  modest  target  range.  Besides,  target 
reflectivity  under  sonar,  signal  coherence  and  overall 
system  designs  have  to  be  well  understood  and  many  experi¬ 
ments  performed  before  the  feasibility  of  sonar  imaging  can 
even  be  assumed. 
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